R Notebook
Step 1: revised catch for 2010-2020 Step 2: tune resilience of all the species (i.e., long-lived, slower growing species = vulnerable…etc…) Step 3: calibrated model with fishing mortality represented (time-averaged catch comparison, model vs observed) = base model Step 4: check base model without fishing to ensure co-existence is still present, observed biomass would not be specified Step 5: Experiments through time - effort functions - uses the base model with no fishing effort as a starting point Step 6: Expose the base model to fishing through time and compare to observed catch time series - start from 0 effort at first time-step to max effort approx around collapse - If the modelled versus observed don’t match, need to adjust effort/or initial biological params and experiment to see what improves - hand tune / optim options Adjust and recalibrate if needed (optim) Use rules to constrain - i.e., - time-averaged biomasses within calibration period must be +- xx % - can optim things like fishing-size selectivity, PPMRs etc
Can we reproduce the time series?
Load libraries
remotes::install_github("sizespectrum/mizerExperimental")Skipping install of 'mizerExperimental' from a github remote, the SHA1 (8279ac0d) has not changed since last install.
Use `force = TRUE` to force installationlibrary(mizerExperimental)
# remotes::install_github("sizespectrum/mizerMR")
# library(mizerMR)
# library(mizer)
library(tidyverse)
library(tictoc)
library(parallel)
# library(plotly)Catch data for all indivudals and summarised as totals per species each year
Maybe this isn’t expressing effort, rather fishing mortality
df_ind_CPUE <- readRDS("ind_catch_weight_BanzareBank_1930_2019_CPUE.rds")
df_CPUE_kg_day <- readRDS("catch_timeseries_BanzareBank_1930_2019_CPUE.rds")
glimpse(df_ind_CPUE)Rows: 107,032
Columns: 9
Groups: Exp.no, Species [104]
$ ID <int> 81681, 81687, 81689, 104868, 104869, 104879, 104881, 104886, 104907, 104908, 104909, 104910, 104916, 104921, 104959, 104960, 104963, 104981, 104984, 104985, 104986, 104987, …
$ Exp.no <int> 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781,…
$ Year <int> 1932, 1932, 1932, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933,…
$ Species <chr> "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "S…
$ Length <dbl> 1737.36, 1493.52, 1645.92, 1645.92, 1584.96, 1493.52, 1645.92, 1615.44, 1706.88, 1554.48, 1615.44, 1645.92, 1645.92, 1706.88, 1676.40, 1615.44, 1615.44, 1584.96, 1615.44, 15…
$ Length.unit <int> 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,…
$ Weight_kg <dbl> 57160.49, 36312.79, 48601.83, 48601.83, 43399.17, 36312.79, 48601.83, 45951.43, 54204.52, 40943.21, 45951.43, 48601.83, 48601.83, 54204.52, 51352.25, 45951.43, 45951.43, 433…
$ date <date> 1932-12-10, 1932-12-10, 1932-12-10, 1933-12-12, 1933-12-12, 1933-12-12, 1933-12-12, 1933-12-13, 1933-12-14, 1933-12-14, 1933-12-14, 1933-12-14, 1933-12-14, 1933-12-15, 1933…
$ Duration <int> 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 1…glimpse(df_CPUE_kg_day)Rows: 132
Columns: 5
Groups: Year [74]
$ Year <int> 1930, 1931, 1932, 1932, 1933, 1933, 1934, 1934, 1935, 1935, 1936, 1936, 1937, 1937, 1938, 1938, 1939, 1939, 1940, 1940, 1946, 1946, 1947, 1947, 1948, 1948, 1949, 1949, 19…
$ Species <chr> "Baleen", "Baleen", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm"…
$ total_catch_kg <dbl> 50155688.13, 335374080.08, 306611735.04, 142075.12, 496515326.66, 3108915.83, 279492543.54, 5993267.00, 220991003.48, 3701390.97, 184531268.07, 4019223.65, 198881583.46, …
$ effort_days <dbl> 79.5265942, 373.4554741, 260.8692295, 0.7931792, 471.8601991, 16.4294765, 309.9332648, 26.8628395, 241.8924453, 18.4724967, 211.9265770, 18.3120140, 218.3889041, 17.84867…
$ CPUE <dbl> 630678.186, 898029.627, 1175346.497, 179121.092, 1052250.916, 189227.930, 901782.981, 223106.236, 913592.002, 200373.076, 870732.074, 219485.615, 910676.228, 206482.577, …Yield for the period that matches Ecopath model, post 2000 annual average
df_yield_2000 <- df_CPUE_kg_day %>%
filter(Year > 2009) %>%
group_by(Species) %>%
summarise(yield = sum(total_catch_kg)/10)
df_yield_2000effort_days/max value plotted by species This might indicate the effort curve required, especially if there is a difference between different species
Relative change in catch pre-1930, could inform the relative change in effort Ask Cami, how could I reconstruct the effort to the initial point of catch when effort = 0…i.e., approx 1850
df_plot %>%
ggplot(aes(x = Year, y = effort_standard, colour = Species)) +
geom_smooth() +
geom_line()
# facet_wrap(~Species) 
Load catch length
catch_lengths <- readRDS("catch_lengths.rds")
glimpse(catch_lengths)Rows: 355
Columns: 5
Groups: Species [4]
$ Species <chr> "Antarctic Minke", "Antarctic Minke", "Antarctic Minke", "Antarctic Minke", "Antarctic Minke", "Antarctic Minke", "Antarctic Minke", "A…
$ species <chr> "minke whales", "minke whales", "minke whales", "minke whales", "minke whales", "minke whales", "minke whales", "minke whales", "minke …
$ catch <int> 1, 5, 5, 9, 16, 17, 27, 30, 41, 36, 46, 43, 26, 33, 26, 32, 34, 38, 55, 65, 95, 69, 112, 197, 149, 167, 390, 239, 305, 614, 503, 503, 1…
$ length <dbl> 480, 490, 500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 700, 710, 720, 730, 740, …
$ dl <dbl> 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,…The value for euphausiids obs_yield is from
calibration_catch_histsoc_1850_2004_regional_models.csv
found at http://portal.sf.utas.edu.au/thredds/catalog/gem/fishmip/ISIMIP3a/InputData/fishing/histsoc/catalog.html.
It is the annual average yield over a 22 year period for the Prydz Bay
Region.
The values of obs_yield for minke whales (2175476136),
orca (19923729), sperm whales (4062177445), and baleen whales
(50341696055) are the annual average yield from 1930 - 2019 from IWC
records of whaling in the Prydz Bay model domain (1,433,028 km2).
Adjust the catch to represent the 2010 - 2020 time period + incorporate toothfish catch if possible (i.e., Stacey’s paper values)
df_ind_CPUE <- readRDS("ind_catch_weight_BanzareBank_1930_2019_CPUE.rds")
df_CPUE_kg_day <- readRDS("catch_timeseries_BanzareBank_1930_2019_CPUE.rds")
glimpse(df_ind_CPUE)Rows: 107,032
Columns: 9
Groups: Exp.no, Species [104]
$ ID <int> 81681, 81687, 81689, 104868, 104869, 104879, 104881, 104886, 104907, 104908, 104909, 104910, 104916, 104921, 104959, 104960, 104963…
$ Exp.no <int> 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781, 5781,…
$ Year <int> 1932, 1932, 1932, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933, 1933,…
$ Species <chr> "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sperm", "Sper…
$ Length <dbl> 1737.36, 1493.52, 1645.92, 1645.92, 1584.96, 1493.52, 1645.92, 1615.44, 1706.88, 1554.48, 1615.44, 1645.92, 1645.92, 1706.88, 1676.…
$ Length.unit <int> 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,…
$ Weight_kg <dbl> 57160.49, 36312.79, 48601.83, 48601.83, 43399.17, 36312.79, 48601.83, 45951.43, 54204.52, 40943.21, 45951.43, 48601.83, 48601.83, 5…
$ date <date> 1932-12-10, 1932-12-10, 1932-12-10, 1933-12-12, 1933-12-12, 1933-12-12, 1933-12-12, 1933-12-13, 1933-12-14, 1933-12-14, 1933-12-14…
$ Duration <int> 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14,…glimpse(df_CPUE_kg_day)Rows: 132
Columns: 5
Groups: Year [74]
$ Year <int> 1930, 1931, 1932, 1932, 1933, 1933, 1934, 1934, 1935, 1935, 1936, 1936, 1937, 1937, 1938, 1938, 1939, 1939, 1940, 1940, 1946, 19…
$ Species <chr> "Baleen", "Baleen", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sperm", "Baleen", "Sp…
$ total_catch_kg <dbl> 50155688.1, 335374080.1, 306611735.0, 142075.1, 496515326.7, 3108915.8, 279492543.5, 5993267.0, 220991003.5, 3701391.0, 18453126…
$ effort_days <dbl> 79.5265942, 373.4554741, 260.8692295, 0.7931792, 471.8601991, 16.4294765, 309.9332648, 26.8628395, 241.8924453, 18.4724967, 211.…
$ CPUE <dbl> 630678.2, 898029.6, 1175346.5, 179121.1, 1052250.9, 189227.9, 901783.0, 223106.2, 913592.0, 200373.1, 870732.1, 219485.6, 910676…Yield for the period that matches Ecopath model, post 2000 annual average
df_yield_2000_2019 <- df_CPUE_kg_day %>%
filter(Year > 2000) %>%
group_by(Species) %>%
summarise(yield = sum(total_catch_kg)/19)
df_yield_2000_2019Add yield in tonnes per km2, same as g m2 and it is consistent with
biomass_observed
1.474341e+12 m^2 for model domain
1.474341e+12/1e+6 = 1474341 km^2
507511 kg of minke whale per year from 2000 - 2019
507.5/1474341 = 0.0003442216 tonnes Minke whale yield per km2 from 2000 - 2019
15619.24 kg of baleen whale per year from 2000 - 2019
15.6/1474341 = 1.0581e-05 tonnes of baleen whale per km2 from 2000 - 2019
obs_yield <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0003442216, 0, 0, 1.0581e-05) # to add a yield_observed column in species_params
biomass_cutoff <- c(0.1, 1, 1, 1, 40, 4500, 1, 1, 200000, 10000, 135000, 600000, 490000, 3650000, 2250000) # to add a biomass_cutoff column in species_params
Steady-state model with no fishing effort calibrated to observed biomass values from McCormack et al. 2020 that represent an average state of the food web from 2010-2020.
# params <- readRDS("stage1_steady_vXX.rds") # Incorrect orca and leopard seal w_max values
so_params <- readRDS("params_16_March_2023.rds") # Updated w_max for orca and leopard seals
# params <- setParams(setFishing(params, initial_effort = 0.2))
# species_params(params)$yield_observed <- obs_yield
# species_params(params)$biomass_cutoff <- biomass_cutoff
# species_params(params)$biomass_cutoffLow <- biomass_cutoff
# species_params(params)$biomass_cutoffHigh <- species_params(params)$w_max
species_params(so_params)$yield_observed <- obs_yield
species_params(so_params)$biomass_cutoff <- biomass_cutoff
species_params(so_params)$biomass_cutoffLow <- biomass_cutoff
species_params(so_params)$biomass_cutoffHigh <- species_params(so_params)$w_max
# species_params(params) |> dplyr::select(species, yield_observed, w_mat, w_max, R_max)
species_params(so_params) |> dplyr::select(species, yield_observed, w_mat, w_max, R_max)# groups |> dplyr::select(species, yield_observed, biomass_observed, biomass_cutoff, w_min)# yield_observed <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0003442216, 0, 0, 1.0581e-05) # to add a yield_observed column in species_params
#
# species_params(params)$yield_observed <- yield_observed# biomass_cutoff <- c(0.1, 1, 1, 1, 40, 4500, 1, 1, 200000, 10000, 135000, 600000, 490000, 3650000, 2250000) # to add a biomass_cutoff column in species_params
#
# species_params(params)$biomass_cutoff <- biomass_cutoffUpdate max and mat size for orca using sealifebase values rather than the previous w_max of 6t
# species_params(params)$species[13]
# species_params(params)$w_max[13]
# species_params(params)$w_mat[13]
#
# species_params(params)$w_max[13] <- 10628034
#
# species_params(params)$w_mat[13] <- 3198855
#
# species_params(params)$w_max[13]
# species_params(params)$w_mat[13]
# species_params(params)$w_mat25[13]Update max and mat size for leopard seal using sealifebase values. Max weight observed reported as 450kg. Using the sealifebase values for L-W conversion and max length results in an estimated max weight of 545875.2g, which is too high above the max weight observed.
# species_params(params)$species[9]
# species_params(params)$w_max[9]
# species_params(params)$w_mat[9]
#
# species_params(params)$w_max[9] <- 450000
#
# species_params(params)$w_max[9]
# species_params(params)$w_mat[9]
# species_params(params)$w_mat25[9]Adjust w_mat values that were changed by default in
newMultispeciesParams
# params_v1 <- params
#
# params_v1@species_params$w_mat[params_v1@species_params$species == "large divers"] <- params_v1@species_params$w_max[params_v1@species_params$species == "large divers"] * 0.9
# params_v1@species_params$w_mat[params_v1@species_params$species == "minke whales"] <- params_v1@species_params$w_max[params_v1@species_params$species == "minke whales"] * 0.9
# params_v1@species_params$w_mat[params_v1@species_params$species == "orca"] <- params_v1@species_params$w_max[params_v1@species_params$species == "orca"] * 0.9
# params_v1@species_params$w_mat[params_v1@species_params$species == "sperm whales"] <- params_v1@species_params$w_max[params_v1@species_params$species == "sperm whales"] * 0.9
#
# params_v1 <- setParams(params_v1)# params_v2 <- params_v1
#
# params_v2@species_params$w_min[params_v2@species_params$species == "small divers"] <- params_v2@species_params$w_mat[params_v2@species_params$species == "small divers"] * 0.85
# params_v2@species_params$w_min[params_v2@species_params$species == "leopard seals"] <- params_v2@species_params$w_mat[params_v2@species_params$species == "leopard seals"] * 0.85
#
# params_v2 <- setParams(params_v2)params <- steady(so_params)
Convergence was achieved in 1.5 years.
Warning: The following species require an unrealistic reproductive efficiency greater than 1: leopard sealsCheck gear params
gear_params(params)Adjust catchability to only select fishing on species with catch data This is a only starting point
gear_params(params)$catchability <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.49, 0, 0, 0.001)Catchability: 1 value Effort changes through time (voyage frequency and length)
gear_params(params)$gear <- c("Main", "Main", "Main", "Main", "Main", "Main",
"Main", "Main", "Main", "Main", "Main", "Main",
"Main", "Main", "Main")
gear_params(params)$l50 <- c(1, 5, 5, 5, 10, 10, 10, 20, 20, 20, 20, 850, 600, 1500, 2200) # values for minke, orca, sperm and baleen are estimated based off `catch_lengths`, while all others are rough guesses, purely as the param won't work with values missing.
gear_params(params)$l25 <- c(0.8, 4, 4, 4, 8, 8, 8, 16, 16, 16, 16, 800, 500, 1400, 1500)
gear_params(params)$sel_func <- c("sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length",
"sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length", "sigmoid_length",
"sigmoid_length", "sigmoid_length", "sigmoid_length")
gear_params(params)params_v02 <- setParams(setFishing(params, initial_effort = 0.2))
params_v02 <- steady(params_v02)Convergence was achieved in 52.5 years.
Warning: The following species require an unrealistic reproductive efficiency greater than 1: euphausiids, leopard seals, minke whalessim_v01 <- project(params_v02, t_max = 100)
[=======>------------------------------------] 19% ETA: 1s
[========>-----------------------------------] 20% ETA: 1s
[========>-----------------------------------] 21% ETA: 1s
[=========>----------------------------------] 22% ETA: 1s
[=========>----------------------------------] 23% ETA: 1s
[=========>----------------------------------] 24% ETA: 1s
[==========>---------------------------------] 25% ETA: 1s
[==========>---------------------------------] 26% ETA: 1s
[===========>--------------------------------] 27% ETA: 1s
[===========>--------------------------------] 28% ETA: 1s
[============>-------------------------------] 29% ETA: 1s
[============>-------------------------------] 30% ETA: 1s
[=============>------------------------------] 31% ETA: 1s
[=============>------------------------------] 32% ETA: 1s
[=============>------------------------------] 33% ETA: 1s
[==============>-----------------------------] 34% ETA: 1s
[==============>-----------------------------] 35% ETA: 1s
[===============>----------------------------] 36% ETA: 1s
[===============>----------------------------] 37% ETA: 1s
[================>---------------------------] 38% ETA: 1s
[================>---------------------------] 39% ETA: 1s
[================>---------------------------] 40% ETA: 1s
[=================>--------------------------] 41% ETA: 1s
[=================>--------------------------] 42% ETA: 1s
[==================>-------------------------] 43% ETA: 1s
[==================>-------------------------] 44% ETA: 1s
[===================>------------------------] 45% ETA: 1s
[===================>------------------------] 46% ETA: 1s
[===================>------------------------] 47% ETA: 1s
[====================>-----------------------] 48% ETA: 1s
[====================>-----------------------] 49% ETA: 1s
[=====================>----------------------] 50% ETA: 1s
[======================>---------------------] 51% ETA: 1s
[======================>---------------------] 52% ETA: 1s
[=======================>--------------------] 53% ETA: 1s
[=======================>--------------------] 54% ETA: 1s
[=======================>--------------------] 55% ETA: 1s
[========================>-------------------] 56% ETA: 1s
[========================>-------------------] 57% ETA: 1s
[=========================>------------------] 58% ETA: 1s
[=========================>------------------] 59% ETA: 1s
[==========================>-----------------] 60% ETA: 1s
[==========================>-----------------] 61% ETA: 1s
[==========================>-----------------] 62% ETA: 1s
[===========================>----------------] 63% ETA: 1s
[===========================>----------------] 64% ETA: 1s
[============================>---------------] 65% ETA: 1s
[============================>---------------] 66% ETA: 1s
[=============================>--------------] 67% ETA: 1s
[=============================>--------------] 68% ETA: 1s
[=============================>--------------] 69% ETA: 0s
[==============================>-------------] 70% ETA: 0s
[==============================>-------------] 71% ETA: 0s
[===============================>------------] 72% ETA: 0s
[===============================>------------] 73% ETA: 0s
[================================>-----------] 74% ETA: 0s
[================================>-----------] 75% ETA: 0s
[=================================>----------] 76% ETA: 0s
[=================================>----------] 77% ETA: 0s
[=================================>----------] 78% ETA: 0s
[==================================>---------] 79% ETA: 0s
[==================================>---------] 80% ETA: 0s
[===================================>--------] 81% ETA: 0s
[===================================>--------] 82% ETA: 0s
[====================================>-------] 83% ETA: 0s
[====================================>-------] 84% ETA: 0s
[====================================>-------] 85% ETA: 0s
[=====================================>------] 86% ETA: 0s
[=====================================>------] 87% ETA: 0s
[======================================>-----] 88% ETA: 0s
[======================================>-----] 89% ETA: 0s
[=======================================>----] 90% ETA: 0s
[=======================================>----] 91% ETA: 0s
[========================================>---] 92% ETA: 0s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[===========================================>] 99% ETA: 0splot(sim_v01)
plotlyBiomass(sim_v01)sim_v01@params@initial_effortMain
0.2 plotYieldVsSize(sim_v01, species = "baleen whales", catch = catch_lengths,
x_var = "Length")
# plotYieldVsSize(sim_v1, species = "sperm whales", catch = catch_lengths,
# x_var = "Length")
#
# plotYieldVsSize(sim_v1, species = "orca", catch = catch_lengths,
# x_var = "Length")
plotYieldVsSize(sim_v01, species = "minke whales", catch = catch_lengths,
x_var = "Length")
getErrorCustom <- function(vary, params, dat, tol = 0.001,
timetorun = 10)
{
params@species_params$R_max[1:15]<-10^vary[1:15] # R_max for 15 species
params@species_params$erepro[1:15]<-vary[16:30] # erepro for 15 species
params@species_params$interaction_resource[1:15] <- vary[31:45] # interaction_resource for 15 species
params <- setParams(params)
# interaction <- params@interaction
# interaction[] <- matrix(vary[28:108],nrow = 9) # stop at 54 if looking only at 3 biggest species
# params <- setInteraction(params,interaction)
params <- projectToSteady(params, distance_func = distanceSSLogN,
tol = tol, t_max = 200, return_sim = F)
sim <- project(params, t_max = timetorun, progress_bar = F)
sim_biomass = rep(0, length(params@species_params$species))
cutoffLow <- params@species_params$biomass_cutoffLow
if (is.null(cutoffLow))
cutoffLow <- rep(0, no_sp)
cutoffLow[is.na(cutoffLow)] <- 0
cutoffHigh <- params@species_params$biomass_cutoffHigh
if (is.null(cutoffHigh))
cutoffHigh <- rep(0, no_sp)
cutoffHigh[is.na(cutoffHigh)] <- 0
for (j in 1:length(sim_biomass)) {
sim_biomass[j] = sum((sim@n[dim(sim@n)[1],j,] * params@w *
params@dw)[params@w >= cutoffLow[j] & cutoffHigh[j] >= params@w])
}
pred <- log(sim_biomass)
dat <- log(dat)
discrep <- pred - dat
discrep <- (sum(discrep^2))
return(discrep)
}
# create set of params for the optimisation process
tic()
params_optim <- params_v02
vary <- c(log10(params_optim@species_params$R_max),
params_optim@species_params$erepro,
params_optim@species_params$interaction_resource)
params_optim<-setParams(params_optim)
# set up workers
noCores <- parallel::detectCores() - 1 # keep some spare core
cl <- parallel::makeCluster(noCores, setup_timeout = 0.5)
setDefaultCluster(cl = cl)
clusterExport(cl, varlist = "cl",envir=environment())
clusterEvalQ(cl, {
library(mizerExperimental)
library(optimParallel)
})
optim_result <- optimParallel::optimParallel(par=vary,getErrorCustom,params=params_optim,
dat = params_optim@species_params$biomass_observed, method ="L-BFGS-B",
lower=c(rep(-15,15),rep(1e-7,15),rep(.1,15)),
upper= c(rep(15,15),rep(1,15),rep(.99,15)),
parallel=list(loginfo=TRUE, forward=TRUE))
stopCluster(cl)
toc()
saveRDS(optim_result, file="optim_result_v00.RDS")
# optim_result <- readRDS("optim_result_v00.RDS")#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max<-10^optim_result$par[1:15]
species_params(params_optim)$erepro<-optim_result$par[16:30]
species_params(params_optim)$interaction_resource <-optim_result$par[31:45]
sim_optim <- project(params_optim, t_max = 2000)
[--------------------------------------------] 0% ETA: 40s
[--------------------------------------------] 1% ETA: 40s
[--------------------------------------------] 1% ETA: 41s
[--------------------------------------------] 1% ETA: 40s
[--------------------------------------------] 1% ETA: 41s
[--------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 2% ETA: 1m
[>-------------------------------------------] 3% ETA: 1m
[>-------------------------------------------] 3% ETA: 50s
[>-------------------------------------------] 3% ETA: 49s
[>-------------------------------------------] 3% ETA: 48s
[>-------------------------------------------] 3% ETA: 1m
[=>------------------------------------------] 3% ETA: 1m
[=>------------------------------------------] 4% ETA: 1m
[=>------------------------------------------] 4% ETA: 50s
[=>------------------------------------------] 4% ETA: 1m
[=>------------------------------------------] 4% ETA: 50s
[=>------------------------------------------] 4% ETA: 49s
[=>------------------------------------------] 4% ETA: 48s
[=>------------------------------------------] 4% ETA: 1m
[=>------------------------------------------] 5% ETA: 1m
[=>------------------------------------------] 5% ETA: 50s
[=>------------------------------------------] 5% ETA: 49s
[=>------------------------------------------] 5% ETA: 48s
[=>------------------------------------------] 6% ETA: 48s
[=>------------------------------------------] 6% ETA: 1m
[=>------------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 49s
[==>-----------------------------------------] 6% ETA: 48s
[==>-----------------------------------------] 6% ETA: 47s
[==>-----------------------------------------] 7% ETA: 47s
[==>-----------------------------------------] 7% ETA: 46s
[==>-----------------------------------------] 7% ETA: 48s
[==>-----------------------------------------] 7% ETA: 47s
[==>-----------------------------------------] 8% ETA: 47s
[==>-----------------------------------------] 8% ETA: 46s
[===>----------------------------------------] 8% ETA: 46s
[===>----------------------------------------] 8% ETA: 45s
[===>----------------------------------------] 9% ETA: 45s
[===>----------------------------------------] 9% ETA: 44s
[===>----------------------------------------] 9% ETA: 46s
[===>----------------------------------------] 9% ETA: 45s
[===>----------------------------------------] 9% ETA: 44s
[===>----------------------------------------] 10% ETA: 44s
[===>----------------------------------------] 10% ETA: 43s
[====>---------------------------------------] 10% ETA: 43s
[====>---------------------------------------] 10% ETA: 44s
[====>---------------------------------------] 11% ETA: 44s
[====>---------------------------------------] 11% ETA: 43s
[====>---------------------------------------] 11% ETA: 42s
[====>---------------------------------------] 12% ETA: 42s
[====>---------------------------------------] 12% ETA: 43s
[====>---------------------------------------] 12% ETA: 42s
[=====>--------------------------------------] 13% ETA: 42s
[=====>--------------------------------------] 13% ETA: 43s
[=====>--------------------------------------] 13% ETA: 42s
[=====>--------------------------------------] 13% ETA: 41s
[=====>--------------------------------------] 14% ETA: 41s
[=====>--------------------------------------] 14% ETA: 40s
[=====>--------------------------------------] 14% ETA: 41s
[=====>--------------------------------------] 15% ETA: 41s
[======>-------------------------------------] 15% ETA: 40s
[======>-------------------------------------] 16% ETA: 39s
[======>-------------------------------------] 16% ETA: 40s
[======>-------------------------------------] 16% ETA: 39s
[======>-------------------------------------] 17% ETA: 39s
[======>-------------------------------------] 17% ETA: 38s
[=======>------------------------------------] 17% ETA: 38s
[=======>------------------------------------] 17% ETA: 39s
[=======>------------------------------------] 18% ETA: 39s
[=======>------------------------------------] 18% ETA: 38s
[=======>------------------------------------] 19% ETA: 38s
[=======>------------------------------------] 19% ETA: 37s
[========>-----------------------------------] 19% ETA: 37s
[========>-----------------------------------] 20% ETA: 37s
[========>-----------------------------------] 21% ETA: 37s
[========>-----------------------------------] 21% ETA: 36s
[========>-----------------------------------] 22% ETA: 36s
[=========>----------------------------------] 22% ETA: 36s
[=========>----------------------------------] 22% ETA: 35s
[=========>----------------------------------] 23% ETA: 35s
[=========>----------------------------------] 24% ETA: 35s
[=========>----------------------------------] 24% ETA: 34s
[==========>---------------------------------] 24% ETA: 34s
[==========>---------------------------------] 25% ETA: 34s
[==========>---------------------------------] 26% ETA: 34s
[===========>--------------------------------] 26% ETA: 34s
[===========>--------------------------------] 26% ETA: 33s
[===========>--------------------------------] 27% ETA: 33s
[===========>--------------------------------] 28% ETA: 33s
[===========>--------------------------------] 28% ETA: 32s
[============>-------------------------------] 28% ETA: 32s
[============>-------------------------------] 29% ETA: 32s
[============>-------------------------------] 29% ETA: 33s
[============>-------------------------------] 29% ETA: 32s
[============>-------------------------------] 30% ETA: 32s
[============>-------------------------------] 31% ETA: 32s
[=============>------------------------------] 31% ETA: 32s
[=============>------------------------------] 31% ETA: 31s
[=============>------------------------------] 31% ETA: 32s
[=============>------------------------------] 31% ETA: 31s
[=============>------------------------------] 32% ETA: 31s
[=============>------------------------------] 33% ETA: 31s
[==============>-----------------------------] 33% ETA: 31s
[==============>-----------------------------] 33% ETA: 30s
[==============>-----------------------------] 34% ETA: 30s
[==============>-----------------------------] 35% ETA: 30s
[==============>-----------------------------] 35% ETA: 29s
[===============>----------------------------] 35% ETA: 29s
[===============>----------------------------] 36% ETA: 29s
[===============>----------------------------] 37% ETA: 29s
[===============>----------------------------] 37% ETA: 28s
[===============>----------------------------] 37% ETA: 29s
[================>---------------------------] 38% ETA: 29s
[================>---------------------------] 38% ETA: 28s
[================>---------------------------] 38% ETA: 29s
[================>---------------------------] 38% ETA: 28s
[================>---------------------------] 39% ETA: 28s
[================>---------------------------] 40% ETA: 28s
[=================>--------------------------] 40% ETA: 28s
[=================>--------------------------] 41% ETA: 28s
[=================>--------------------------] 41% ETA: 27s
[=================>--------------------------] 41% ETA: 28s
[=================>--------------------------] 42% ETA: 28s
[=================>--------------------------] 42% ETA: 27s
[==================>-------------------------] 42% ETA: 27s
[==================>-------------------------] 43% ETA: 27s
[==================>-------------------------] 44% ETA: 27s
[===================>------------------------] 44% ETA: 27s
[===================>------------------------] 44% ETA: 26s
[===================>------------------------] 45% ETA: 26s
[===================>------------------------] 46% ETA: 26s
[===================>------------------------] 47% ETA: 26s
[====================>-----------------------] 47% ETA: 26s
[====================>-----------------------] 47% ETA: 25s
[====================>-----------------------] 47% ETA: 26s
[====================>-----------------------] 47% ETA: 25s
[====================>-----------------------] 48% ETA: 25s
[====================>-----------------------] 49% ETA: 25s
[=====================>----------------------] 49% ETA: 25s
[=====================>----------------------] 49% ETA: 24s
[=====================>----------------------] 50% ETA: 24s
[=====================>----------------------] 51% ETA: 24s
[======================>---------------------] 51% ETA: 24s
[======================>---------------------] 52% ETA: 24s
[======================>---------------------] 52% ETA: 23s
[======================>---------------------] 53% ETA: 23s
[=======================>--------------------] 53% ETA: 23s
[=======================>--------------------] 54% ETA: 23s
[=======================>--------------------] 54% ETA: 22s
[=======================>--------------------] 54% ETA: 23s
[=======================>--------------------] 54% ETA: 22s
[=======================>--------------------] 55% ETA: 22s
[=======================>--------------------] 56% ETA: 22s
[========================>-------------------] 56% ETA: 22s
[========================>-------------------] 57% ETA: 22s
[========================>-------------------] 57% ETA: 21s
[========================>-------------------] 58% ETA: 21s
[=========================>------------------] 58% ETA: 21s
[=========================>------------------] 59% ETA: 21s
[=========================>------------------] 59% ETA: 20s
[=========================>------------------] 60% ETA: 20s
[==========================>-----------------] 60% ETA: 20s
[==========================>-----------------] 61% ETA: 20s
[==========================>-----------------] 61% ETA: 19s
[==========================>-----------------] 62% ETA: 19s
[===========================>----------------] 63% ETA: 19s
[===========================>----------------] 63% ETA: 18s
[===========================>----------------] 63% ETA: 19s
[===========================>----------------] 63% ETA: 18s
[===========================>----------------] 64% ETA: 18s
[===========================>----------------] 65% ETA: 18s
[============================>---------------] 65% ETA: 18s
[============================>---------------] 65% ETA: 17s
[============================>---------------] 66% ETA: 17s
[============================>---------------] 67% ETA: 17s
[=============================>--------------] 67% ETA: 16s
[=============================>--------------] 68% ETA: 16s
[=============================>--------------] 69% ETA: 16s
[=============================>--------------] 69% ETA: 15s
[=============================>--------------] 69% ETA: 16s
[=============================>--------------] 69% ETA: 15s
[==============================>-------------] 69% ETA: 15s
[==============================>-------------] 70% ETA: 15s
[==============================>-------------] 71% ETA: 15s
[==============================>-------------] 71% ETA: 14s
[==============================>-------------] 72% ETA: 14s
[===============================>------------] 72% ETA: 14s
[===============================>------------] 73% ETA: 14s
[===============================>------------] 73% ETA: 13s
[===============================>------------] 74% ETA: 13s
[================================>-----------] 74% ETA: 13s
[================================>-----------] 75% ETA: 13s
[================================>-----------] 75% ETA: 12s
[================================>-----------] 76% ETA: 12s
[=================================>----------] 76% ETA: 12s
[=================================>----------] 77% ETA: 12s
[=================================>----------] 77% ETA: 11s
[=================================>----------] 78% ETA: 11s
[==================================>---------] 78% ETA: 11s
[==================================>---------] 79% ETA: 11s
[==================================>---------] 79% ETA: 10s
[==================================>---------] 80% ETA: 10s
[==================================>---------] 81% ETA: 10s
[===================================>--------] 81% ETA: 10s
[===================================>--------] 81% ETA: 9s
[===================================>--------] 82% ETA: 9s
[===================================>--------] 83% ETA: 9s
[====================================>-------] 83% ETA: 9s
[====================================>-------] 83% ETA: 8s
[====================================>-------] 84% ETA: 8s
[====================================>-------] 85% ETA: 8s
[====================================>-------] 85% ETA: 7s
[=====================================>------] 85% ETA: 7s
[=====================================>------] 86% ETA: 7s
[=====================================>------] 87% ETA: 7s
[=====================================>------] 87% ETA: 6s
[======================================>-----] 88% ETA: 6s
[======================================>-----] 89% ETA: 6s
[======================================>-----] 89% ETA: 5s
[======================================>-----] 90% ETA: 5s
[=======================================>----] 90% ETA: 5s
[=======================================>----] 91% ETA: 5s
[=======================================>----] 91% ETA: 4s
[=======================================>----] 92% ETA: 4s
[========================================>---] 92% ETA: 4s
[========================================>---] 93% ETA: 4s
[========================================>---] 93% ETA: 3s
[========================================>---] 94% ETA: 3s
[=========================================>--] 94% ETA: 3s
[=========================================>--] 95% ETA: 3s
[=========================================>--] 95% ETA: 2s
[=========================================>--] 96% ETA: 2s
[=========================================>--] 97% ETA: 2s
[==========================================>-] 97% ETA: 2s
[==========================================>-] 97% ETA: 1s
[==========================================>-] 98% ETA: 1s
[==========================================>-] 99% ETA: 1s
[===========================================>] 99% ETA: 1s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_optim)
#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max [1] 4.816091e+00 2.694485e-03 2.079331e-03 1.087521e-04 1.486522e-07 1.727671e-06 6.342514e-07 5.682170e-07 3.435600e-10 1.145037e-08 1.296982e-09
[12] 9.091376e-10 7.838551e-12 3.411213e-11 1.439994e-10species_params(params_optim)$erepro [1] 1.0000000 0.9979727 0.9964291 0.9987793 0.9996089 0.7126939 0.9956963 0.9875731 0.9986068 0.8251681 0.7820067 0.9913709 0.7658567 0.7118458
[15] 0.7590756species_params(params_optim)$interaction_resource [1] 0.9900000 0.9900000 0.9900000 0.9900000 0.9408792 0.9338620 0.9900000 0.9900000 0.9655430 0.9863187 0.9056183 0.9847589 0.9851372 0.9199649
[15] 0.9309803# this function adds a lower boundary to selected size
plotBiomassObservedVsModelCustom <- function (object, species = NULL, ratio = FALSE, log_scale = TRUE,
return_data = FALSE, labels = TRUE, show_unobserved = FALSE)
{
if (is(object, "MizerSim")) {
params = object@params
n <- finalN(object)
}
else if (is(object, "MizerParams")) {
params = object
n <- initialN(params)
}
else {
stop("You have not provided a valid mizerSim or mizerParams object.")
}
sp_params <- params@species_params
species = valid_species_arg(object, species)
if (length(species) == 0)
stop("No species selected, please fix.")
row_select = match(species, sp_params$species)
if (!"biomass_observed" %in% names(sp_params)) {
stop("You have not provided values for the column 'biomass_observed' ",
"in the mizerParams/mizerSim object.")
}
else if (!is.numeric(sp_params$biomass_observed)) {
stop("The column 'biomass_observed' in the mizerParams/mizerSim object",
" is not numeric, please fix.")
}
else {
biomass_observed = sp_params$biomass_observed
}
cutoffLow <- sp_params$biomass_cutoffLow[row_select]
if (is.null(cutoffLow)) {
cutoffLow = rep(0, length(species))
}
else if (!is.numeric(cutoffLow)) {
stop("params@species_params$biomass_cutoffLow is not numeric, \",\n \"please fix.")
}
cutoffLow[is.na(cutoffLow)] <- 0
cutoffHigh <- sp_params$biomass_cutoffHigh[row_select]
if (is.null(cutoffHigh)) {
cutoffHigh = rep(0, length(species))
}
else if (!is.numeric(cutoffHigh)) {
stop("params@species_params$biomass_cutoffHigh is not numeric, \",\n \"please fix.")
}
cutoffHigh[is.na(cutoffHigh)] <- 0
sim_biomass = rep(0, length(species))
for (j in 1:length(species)) {
sim_biomass[j] = sum((n[row_select[j], ] * params@w *
params@dw)[params@w >= cutoffLow[j] & cutoffHigh[j] >= params@w])
}
dummy = data.frame(species = species, model = sim_biomass,
observed = biomass_observed[row_select]) %>% mutate(species = factor(species,
levels = species), is_observed = !is.na(observed) & observed >
0, observed = case_when(is_observed ~ observed, !is_observed ~
model), ratio = model/observed)
if (sum(dummy$is_observed) == 0) {
cat(paste("There are no observed biomasses to compare to model,",
"only plotting model biomasses.", sep = "\n"))
}
if (!show_unobserved) {
dummy <- filter(dummy, is_observed)
}
if (return_data == TRUE)
return(dummy)
tre <- round(sum(abs(1 - dummy$ratio)), digits = 3)
caption <- paste0("Total relative error = ", tre)
if (any(!dummy$is_observed)) {
caption <- paste(caption, "\n Open circles represent species without biomass observation.")
}
if (ratio == FALSE) {
gg <- ggplot(data = dummy, aes(x = observed, y = model,
colour = species, shape = is_observed)) + geom_abline(aes(intercept = 0,
slope = 1), colour = "purple", linetype = "dashed",
size = 1.3) + geom_point(size = 3) + labs(y = "model biomass [g]") +
coord_cartesian(ylim = range(dummy$model, dummy$observed))
}
else {
gg <- ggplot(data = dummy, aes(x = observed, y = ratio,
colour = species, shape = is_observed)) + geom_hline(aes(yintercept = 1),
linetype = "dashed", colour = "purple",
size = 1.3) + geom_point(size = 3) + labs(y = "model biomass / observed biomass") +
coord_cartesian(ylim = range(dummy$ratio))
}
gg <- gg + labs(x = "observed biomass [g]", caption = caption) +
scale_colour_manual(values = getColours(params)[dummy$species]) +
scale_shape_manual(values = c(`TRUE` = 19, `FALSE` = 1)) +
guides(shape = "none")
if (log_scale == TRUE & ratio == FALSE) {
gg = gg + scale_x_log10() + scale_y_log10()
}
if (log_scale == TRUE & ratio == TRUE) {
gg = gg + scale_x_log10()
}
if (labels == TRUE) {
gg = gg + ggrepel::geom_label_repel(aes(label = species),
box.padding = 0.35, point.padding = 0.5, segment.color = "grey50",
show.legend = FALSE, max.overlaps = Inf, seed = 42)
}
gg
}plotBiomassObservedVsModelCustom(sim_optim)
Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
Please use `linewidth` instead.
Use tuneParams() to investigate the species with erepro values too high
params_tuned_v01@species_params$erepro
[1] 1.0038283 0.9994687 0.9976619 0.9995862 1.0036193 0.7142939 0.9962189 0.9879108 0.9997113 0.8253208 0.7841239 0.9990174 0.7665787 0.7098369
[15] 0.7553690# create set of params for the optimisation process
tic()
params_optim <- params_tuned_v01
vary <- c(log10(params_optim@species_params$R_max),
params_optim@species_params$erepro,
params_optim@species_params$interaction_resource)
params_optim<-setParams(params_optim)
# set up workers
noCores <- parallel::detectCores() - 1 # keep some spare core
cl <- parallel::makeCluster(noCores, setup_timeout = 0.5)
setDefaultCluster(cl = cl)
clusterExport(cl, varlist = "cl",envir=environment())
clusterEvalQ(cl, {
library(mizerExperimental)
library(optimParallel)
})
optim_result <- optimParallel::optimParallel(par=vary,getErrorCustom,params=params_optim,
dat = params_optim@species_params$biomass_observed, method ="L-BFGS-B",
lower=c(rep(-15,15),rep(1e-7,15),rep(.1,15)),
upper= c(rep(15,15),rep(0.99,15),rep(.99,15)),
parallel=list(loginfo=TRUE, forward=TRUE))
stopCluster(cl)
toc()
saveRDS(optim_result, file="optim_result_v01.RDS")
# optim_result <- readRDS("optim_result_v00.RDS")#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max<-10^optim_result$par[1:15]
species_params(params_optim)$erepro<-optim_result$par[16:30]
species_params(params_optim)$interaction_resource <-optim_result$par[31:45]
species_params(params_optim)$R_max [1] 4.816945e+00 2.694890e-03 2.077221e-03 1.087269e-04 1.467918e-07 1.735034e-06 6.331959e-07 5.675777e-07 3.507060e-10 1.153278e-08 1.226957e-09
[12] 8.590879e-10 7.806507e-12 2.801063e-11 1.423327e-10species_params(params_optim)$erepro [1] 0.9900000 0.9900000 0.9900000 0.9900000 0.9898627 0.6810777 0.9900000 0.9879108 0.9875105 0.8256844 0.8030409 0.9867368 0.7664345 0.7315449
[15] 0.7560960species_params(params_optim)$interaction_resource [1] 0.9897483 0.9872286 0.9900000 0.9898699 0.9006595 0.8409188 0.9900000 0.9900000 0.9529289 0.9803054 0.7903730 0.9887297 0.9833562 0.6553079
[15] 0.8602099sim_optim_v02 <- project(params_optim, t_max = 1000)
[>-------------------------------------------] 1% ETA: 15s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 2% ETA: 17s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 3% ETA: 16s
[>-------------------------------------------] 3% ETA: 23s
[>-------------------------------------------] 3% ETA: 22s
[=>------------------------------------------] 3% ETA: 22s
[=>------------------------------------------] 4% ETA: 21s
[=>------------------------------------------] 4% ETA: 20s
[=>------------------------------------------] 5% ETA: 20s
[=>------------------------------------------] 5% ETA: 19s
[=>------------------------------------------] 6% ETA: 19s
[==>-----------------------------------------] 6% ETA: 19s
[==>-----------------------------------------] 7% ETA: 19s
[==>-----------------------------------------] 7% ETA: 18s
[==>-----------------------------------------] 7% ETA: 19s
[==>-----------------------------------------] 7% ETA: 18s
[==>-----------------------------------------] 7% ETA: 19s
[==>-----------------------------------------] 7% ETA: 18s
[==>-----------------------------------------] 8% ETA: 18s
[===>----------------------------------------] 8% ETA: 18s
[===>----------------------------------------] 9% ETA: 18s
[===>----------------------------------------] 9% ETA: 17s
[===>----------------------------------------] 10% ETA: 17s
[====>---------------------------------------] 10% ETA: 17s
[====>---------------------------------------] 10% ETA: 18s
[====>---------------------------------------] 11% ETA: 18s
[====>---------------------------------------] 11% ETA: 17s
[====>---------------------------------------] 12% ETA: 17s
[=====>--------------------------------------] 13% ETA: 17s
[=====>--------------------------------------] 13% ETA: 16s
[=====>--------------------------------------] 14% ETA: 16s
[=====>--------------------------------------] 15% ETA: 16s
[======>-------------------------------------] 15% ETA: 16s
[======>-------------------------------------] 16% ETA: 16s
[======>-------------------------------------] 16% ETA: 15s
[======>-------------------------------------] 16% ETA: 16s
[======>-------------------------------------] 16% ETA: 15s
[======>-------------------------------------] 17% ETA: 15s
[=======>------------------------------------] 17% ETA: 15s
[=======>------------------------------------] 18% ETA: 15s
[=======>------------------------------------] 19% ETA: 15s
[=======>------------------------------------] 19% ETA: 14s
[=======>------------------------------------] 19% ETA: 15s
[========>-----------------------------------] 19% ETA: 15s
[========>-----------------------------------] 20% ETA: 15s
[========>-----------------------------------] 21% ETA: 15s
[========>-----------------------------------] 21% ETA: 14s
[========>-----------------------------------] 22% ETA: 14s
[=========>----------------------------------] 22% ETA: 14s
[=========>----------------------------------] 23% ETA: 14s
[=========>----------------------------------] 24% ETA: 14s
[==========>---------------------------------] 24% ETA: 14s
[==========>---------------------------------] 24% ETA: 13s
[==========>---------------------------------] 25% ETA: 13s
[==========>---------------------------------] 26% ETA: 13s
[===========>--------------------------------] 26% ETA: 13s
[===========>--------------------------------] 27% ETA: 13s
[===========>--------------------------------] 28% ETA: 13s
[============>-------------------------------] 28% ETA: 13s
[============>-------------------------------] 29% ETA: 13s
[============>-------------------------------] 29% ETA: 12s
[============>-------------------------------] 30% ETA: 12s
[============>-------------------------------] 31% ETA: 12s
[=============>------------------------------] 31% ETA: 12s
[=============>------------------------------] 32% ETA: 12s
[=============>------------------------------] 33% ETA: 12s
[==============>-----------------------------] 33% ETA: 11s
[==============>-----------------------------] 34% ETA: 11s
[==============>-----------------------------] 35% ETA: 11s
[===============>----------------------------] 35% ETA: 11s
[===============>----------------------------] 36% ETA: 11s
[===============>----------------------------] 37% ETA: 11s
[================>---------------------------] 38% ETA: 11s
[================>---------------------------] 38% ETA: 10s
[================>---------------------------] 39% ETA: 10s
[================>---------------------------] 40% ETA: 10s
[=================>--------------------------] 40% ETA: 10s
[=================>--------------------------] 41% ETA: 10s
[=================>--------------------------] 42% ETA: 10s
[==================>-------------------------] 42% ETA: 10s
[==================>-------------------------] 43% ETA: 10s
[==================>-------------------------] 44% ETA: 9s
[===================>------------------------] 44% ETA: 9s
[===================>------------------------] 45% ETA: 9s
[===================>------------------------] 46% ETA: 9s
[===================>------------------------] 47% ETA: 9s
[====================>-----------------------] 47% ETA: 9s
[====================>-----------------------] 48% ETA: 9s
[====================>-----------------------] 49% ETA: 8s
[=====================>----------------------] 49% ETA: 8s
[=====================>----------------------] 50% ETA: 8s
[=====================>----------------------] 51% ETA: 8s
[======================>---------------------] 51% ETA: 8s
[======================>---------------------] 52% ETA: 8s
[======================>---------------------] 53% ETA: 8s
[=======================>--------------------] 53% ETA: 8s
[=======================>--------------------] 54% ETA: 8s
[=======================>--------------------] 55% ETA: 8s
[=======================>--------------------] 55% ETA: 7s
[=======================>--------------------] 56% ETA: 7s
[========================>-------------------] 56% ETA: 7s
[========================>-------------------] 57% ETA: 7s
[========================>-------------------] 58% ETA: 7s
[=========================>------------------] 58% ETA: 7s
[=========================>------------------] 59% ETA: 7s
[=========================>------------------] 60% ETA: 7s
[==========================>-----------------] 60% ETA: 7s
[==========================>-----------------] 61% ETA: 7s
[==========================>-----------------] 62% ETA: 6s
[===========================>----------------] 63% ETA: 6s
[===========================>----------------] 64% ETA: 6s
[===========================>----------------] 65% ETA: 6s
[============================>---------------] 65% ETA: 6s
[============================>---------------] 66% ETA: 6s
[============================>---------------] 67% ETA: 6s
[============================>---------------] 67% ETA: 5s
[=============================>--------------] 67% ETA: 5s
[=============================>--------------] 68% ETA: 5s
[=============================>--------------] 69% ETA: 5s
[==============================>-------------] 69% ETA: 5s
[==============================>-------------] 70% ETA: 5s
[==============================>-------------] 71% ETA: 5s
[==============================>-------------] 72% ETA: 5s
[===============================>------------] 72% ETA: 5s
[===============================>------------] 73% ETA: 5s
[===============================>------------] 73% ETA: 4s
[===============================>------------] 74% ETA: 4s
[================================>-----------] 74% ETA: 4s
[================================>-----------] 75% ETA: 4s
[================================>-----------] 76% ETA: 4s
[=================================>----------] 76% ETA: 4s
[=================================>----------] 77% ETA: 4s
[=================================>----------] 78% ETA: 4s
[==================================>---------] 78% ETA: 4s
[==================================>---------] 79% ETA: 4s
[==================================>---------] 79% ETA: 3s
[==================================>---------] 80% ETA: 3s
[==================================>---------] 81% ETA: 3s
[===================================>--------] 81% ETA: 3s
[===================================>--------] 82% ETA: 3s
[===================================>--------] 83% ETA: 3s
[====================================>-------] 83% ETA: 3s
[====================================>-------] 84% ETA: 3s
[====================================>-------] 85% ETA: 3s
[====================================>-------] 85% ETA: 2s
[=====================================>------] 85% ETA: 2s
[=====================================>------] 86% ETA: 2s
[=====================================>------] 87% ETA: 2s
[======================================>-----] 88% ETA: 2s
[======================================>-----] 89% ETA: 2s
[======================================>-----] 90% ETA: 2s
[=======================================>----] 90% ETA: 2s
[=======================================>----] 91% ETA: 2s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 1s
[========================================>---] 94% ETA: 1s
[=========================================>--] 94% ETA: 1s
[=========================================>--] 95% ETA: 1s
[=========================================>--] 96% ETA: 1s
[=========================================>--] 97% ETA: 1s
[==========================================>-] 97% ETA: 1s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_optim_v02)
params_v03 <- steady(params_optim)
Simulation run did not converge after 99 years. Value returned by the distance function was: 0.011566504266491params_loop <- params_v03 |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady() |>
matchBiomasses() |> steady() |> matchBiomasses() |> steady()Convergence was achieved in 10.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.
Convergence was achieved in 1.5 years.params_loop@species_params$erepro [1] 0.94874184 0.69740266 0.74789288 0.68164999 0.40986421 0.24796060 0.71872306 0.89462821 0.63623877 0.54381592 0.25947925 0.65810143 0.49211255
[14] 0.12865081 0.07863463sim_v02 <- project(params_loop, t_max = 2000)
Save progress
Try to ween the species of the resource spectra
params_loop@resource_params$kappa
[1] 19.76374params_v04@resource_params$w_pp_cutoff
[1] 10
plotBiomass(sim_v03)
plotBiomassObservedVsModelCustom(sim_v03)
plotBiomassObservedVsModel(sim_v03)
plotDiet(params_v04)
getErrorCustom_v02 <- function(vary, params, dat, tol = 0.001,
timetorun = 10)
{
params@species_params$R_max[1:15]<-10^vary[1:15] # R_max for 15 species
params@species_params$erepro[1:15]<-vary[16:30] # erepro for 15 species
params@species_params$interaction_resource[1:15] <- vary[31:45] # interaction_resource for 15 species
params@resource_params$w_pp_cutoff[1] <- vary[46] # interaction_resource for 15 species
params <- setParams(params)
# interaction <- params@interaction
# interaction[] <- matrix(vary[28:108],nrow = 9) # stop at 54 if looking only at 3 biggest species
# params <- setInteraction(params,interaction)
params <- projectToSteady(params, distance_func = distanceSSLogN,
tol = tol, t_max = 200, return_sim = F)
sim <- project(params, t_max = timetorun, progress_bar = F)
sim_biomass = rep(0, length(params@species_params$species))
cutoffLow <- params@species_params$biomass_cutoffLow
if (is.null(cutoffLow))
cutoffLow <- rep(0, no_sp)
cutoffLow[is.na(cutoffLow)] <- 0
cutoffHigh <- params@species_params$biomass_cutoffHigh
if (is.null(cutoffHigh))
cutoffHigh <- rep(0, no_sp)
cutoffHigh[is.na(cutoffHigh)] <- 0
for (j in 1:length(sim_biomass)) {
sim_biomass[j] = sum((sim@n[dim(sim@n)[1],j,] * params@w *
params@dw)[params@w >= cutoffLow[j] & cutoffHigh[j] >= params@w])
}
pred <- log(sim_biomass)
dat <- log(dat)
discrep <- pred - dat
discrep <- (sum(discrep^2))
return(discrep)
}
# create set of params for the optimisation process
tic()
params_optim <- params_v04
vary <- c(log10(params_optim@species_params$R_max),
params_optim@species_params$erepro,
params_optim@species_params$interaction_resource,
params_optim@resource_params$w_pp_cutoff)
params_optim<-setParams(params_optim)
# set up workers
noCores <- parallel::detectCores() - 1 # keep some spare core
cl <- parallel::makeCluster(noCores, setup_timeout = 0.5)
setDefaultCluster(cl = cl)
clusterExport(cl, varlist = "cl",envir=environment())
clusterEvalQ(cl, {
library(mizerExperimental)
library(optimParallel)
})
optim_result <- optimParallel::optimParallel(par=vary,getErrorCustom_v02,params=params_optim,
dat = params_optim@species_params$biomass_observed, method ="L-BFGS-B",
lower=c(rep(-15,15),rep(1e-7,15),rep(.1,15), 0.1),
upper= c(rep(15,15),rep(0.99,15),rep(.99,15),1000),
parallel=list(loginfo=TRUE, forward=TRUE))
stopCluster(cl)
toc()
saveRDS(optim_result, file="optim_result_v02.RDS")
# optim_result <- readRDS("optim_result_v00.RDS")#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max<-10^optim_result$par[1:15]
species_params(params_optim)$erepro<-optim_result$par[16:30]
species_params(params_optim)$interaction_resource <-optim_result$par[31:45]
resource_params(params_optim)$w_pp_cutoff <- optim_result$par[46]
species_params(params_optim)$R_max [1] 4.824845e+00 2.704878e-03 2.078290e-03 1.087379e-04 1.452868e-07 2.373280e-06 6.334551e-07 5.664204e-07
[9] 4.012486e-10 1.147254e-08 1.282038e-09 1.001733e-09 8.019699e-12 2.460401e-11 1.411257e-10species_params(params_optim)$erepro [1] 0.9487424 0.6974029 0.7478930 0.6816501 0.4110420 0.3877134 0.7187231 0.8946282 0.7424363 0.5446118 0.4133560
[12] 0.7679180 0.4922805 0.4796008 0.1651285species_params(params_optim)$interaction_resource [1] 0.9880931 0.9748302 0.9900000 0.9900000 0.8636127 0.8778679 0.9900000 0.9900000 0.9751581 0.9745590 0.7326587
[12] 0.9900000 0.9821797 0.4116474 0.7857996sim_v04 <- project(params_optim, t_max = 1000)
[>-------------------------------------------] 1% ETA: 15s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 3% ETA: 15s
[=>------------------------------------------] 3% ETA: 15s
[=>------------------------------------------] 4% ETA: 15s
[=>------------------------------------------] 4% ETA: 14s
[=>------------------------------------------] 4% ETA: 15s
[=>------------------------------------------] 4% ETA: 14s
[=>------------------------------------------] 5% ETA: 15s
[=>------------------------------------------] 5% ETA: 14s
[=>------------------------------------------] 6% ETA: 14s
[==>-----------------------------------------] 6% ETA: 14s
[==>-----------------------------------------] 7% ETA: 14s
[==>-----------------------------------------] 8% ETA: 14s
[===>----------------------------------------] 8% ETA: 13s
[===>----------------------------------------] 8% ETA: 14s
[===>----------------------------------------] 8% ETA: 13s
[===>----------------------------------------] 9% ETA: 13s
[===>----------------------------------------] 10% ETA: 15s
[====>---------------------------------------] 10% ETA: 15s
[====>---------------------------------------] 11% ETA: 15s
[====>---------------------------------------] 11% ETA: 14s
[====>---------------------------------------] 12% ETA: 14s
[=====>--------------------------------------] 13% ETA: 14s
[=====>--------------------------------------] 13% ETA: 13s
[=====>--------------------------------------] 14% ETA: 13s
[=====>--------------------------------------] 14% ETA: 14s
[=====>--------------------------------------] 14% ETA: 13s
[=====>--------------------------------------] 15% ETA: 13s
[======>-------------------------------------] 15% ETA: 13s
[======>-------------------------------------] 16% ETA: 13s
[======>-------------------------------------] 17% ETA: 13s
[=======>------------------------------------] 17% ETA: 13s
[=======>------------------------------------] 17% ETA: 12s
[=======>------------------------------------] 18% ETA: 12s
[=======>------------------------------------] 19% ETA: 12s
[========>-----------------------------------] 19% ETA: 12s
[========>-----------------------------------] 20% ETA: 12s
[========>-----------------------------------] 21% ETA: 12s
[========>-----------------------------------] 21% ETA: 11s
[========>-----------------------------------] 22% ETA: 11s
[=========>----------------------------------] 22% ETA: 11s
[=========>----------------------------------] 23% ETA: 12s
[=========>----------------------------------] 24% ETA: 12s
[==========>---------------------------------] 24% ETA: 12s
[==========>---------------------------------] 24% ETA: 11s
[==========>---------------------------------] 24% ETA: 12s
[==========>---------------------------------] 24% ETA: 11s
[==========>---------------------------------] 25% ETA: 11s
[==========>---------------------------------] 26% ETA: 11s
[===========>--------------------------------] 26% ETA: 11s
[===========>--------------------------------] 27% ETA: 11s
[===========>--------------------------------] 28% ETA: 11s
[============>-------------------------------] 28% ETA: 11s
[============>-------------------------------] 29% ETA: 11s
[============>-------------------------------] 29% ETA: 10s
[============>-------------------------------] 30% ETA: 10s
[============>-------------------------------] 31% ETA: 10s
[=============>------------------------------] 31% ETA: 10s
[=============>------------------------------] 32% ETA: 10s
[=============>------------------------------] 33% ETA: 10s
[==============>-----------------------------] 33% ETA: 10s
[==============>-----------------------------] 34% ETA: 10s
[==============>-----------------------------] 34% ETA: 9s
[==============>-----------------------------] 35% ETA: 9s
[===============>----------------------------] 35% ETA: 9s
[===============>----------------------------] 36% ETA: 9s
[===============>----------------------------] 37% ETA: 9s
[================>---------------------------] 38% ETA: 9s
[================>---------------------------] 39% ETA: 9s
[================>---------------------------] 40% ETA: 9s
[=================>--------------------------] 40% ETA: 9s
[=================>--------------------------] 41% ETA: 9s
[=================>--------------------------] 41% ETA: 8s
[=================>--------------------------] 42% ETA: 8s
[==================>-------------------------] 42% ETA: 8s
[==================>-------------------------] 43% ETA: 8s
[==================>-------------------------] 44% ETA: 8s
[===================>------------------------] 44% ETA: 8s
[===================>------------------------] 45% ETA: 8s
[===================>------------------------] 46% ETA: 8s
[===================>------------------------] 47% ETA: 8s
[====================>-----------------------] 47% ETA: 8s
[====================>-----------------------] 47% ETA: 7s
[====================>-----------------------] 48% ETA: 7s
[====================>-----------------------] 49% ETA: 7s
[=====================>----------------------] 49% ETA: 7s
[=====================>----------------------] 50% ETA: 7s
[=====================>----------------------] 51% ETA: 7s
[======================>---------------------] 51% ETA: 7s
[======================>---------------------] 52% ETA: 7s
[======================>---------------------] 53% ETA: 7s
[=======================>--------------------] 53% ETA: 7s
[=======================>--------------------] 54% ETA: 7s
[=======================>--------------------] 54% ETA: 6s
[=======================>--------------------] 55% ETA: 6s
[=======================>--------------------] 56% ETA: 6s
[========================>-------------------] 56% ETA: 6s
[========================>-------------------] 57% ETA: 6s
[========================>-------------------] 58% ETA: 6s
[=========================>------------------] 58% ETA: 6s
[=========================>------------------] 59% ETA: 6s
[=========================>------------------] 60% ETA: 6s
[==========================>-----------------] 60% ETA: 6s
[==========================>-----------------] 61% ETA: 5s
[==========================>-----------------] 62% ETA: 5s
[===========================>----------------] 63% ETA: 5s
[===========================>----------------] 64% ETA: 5s
[===========================>----------------] 65% ETA: 5s
[============================>---------------] 65% ETA: 5s
[============================>---------------] 66% ETA: 5s
[============================>---------------] 67% ETA: 5s
[=============================>--------------] 67% ETA: 5s
[=============================>--------------] 68% ETA: 5s
[=============================>--------------] 68% ETA: 4s
[=============================>--------------] 69% ETA: 4s
[==============================>-------------] 69% ETA: 4s
[==============================>-------------] 70% ETA: 4s
[==============================>-------------] 71% ETA: 4s
[==============================>-------------] 72% ETA: 4s
[===============================>------------] 72% ETA: 4s
[===============================>------------] 73% ETA: 4s
[===============================>------------] 74% ETA: 4s
[================================>-----------] 74% ETA: 4s
[================================>-----------] 75% ETA: 4s
[================================>-----------] 75% ETA: 3s
[================================>-----------] 76% ETA: 3s
[=================================>----------] 76% ETA: 3s
[=================================>----------] 77% ETA: 3s
[=================================>----------] 78% ETA: 3s
[==================================>---------] 78% ETA: 3s
[==================================>---------] 79% ETA: 3s
[==================================>---------] 80% ETA: 3s
[==================================>---------] 81% ETA: 3s
[===================================>--------] 81% ETA: 3s
[===================================>--------] 82% ETA: 3s
[===================================>--------] 82% ETA: 2s
[===================================>--------] 83% ETA: 2s
[====================================>-------] 83% ETA: 2s
[====================================>-------] 84% ETA: 2s
[====================================>-------] 85% ETA: 2s
[=====================================>------] 85% ETA: 2s
[=====================================>------] 86% ETA: 2s
[=====================================>------] 87% ETA: 2s
[======================================>-----] 88% ETA: 2s
[======================================>-----] 89% ETA: 2s
[======================================>-----] 89% ETA: 1s
[======================================>-----] 90% ETA: 1s
[=======================================>----] 90% ETA: 1s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 1s
[========================================>---] 94% ETA: 1s
[=========================================>--] 94% ETA: 1s
[=========================================>--] 95% ETA: 1s
[=========================================>--] 96% ETA: 1s
[=========================================>--] 96% ETA: 0s
[=========================================>--] 97% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_v04)
plotBiomassObservedVsModelCustom(sim_v04)
plotBiomassObservedVsModel(sim_v04)
plotDiet(params_optim)
params_tuned_v02@species_params$erepro
[1] 0.9485557 0.6971840 0.7477216 0.6816138 0.4110374 0.3877374 0.7186214 0.8940243 0.7432162 0.5455491 0.4207180
[12] 0.7685176 0.4934246 0.4891521 0.1642704species w_min w_mat w_max beta k_vb h min_depth max_depth water.column.use p_time_prydz pc_annual_offspring 5 salps 3.162278e-05 0.2511886 25.11886 1000 NA 33 0 500 non DVM 1 NA biomass_observed 5 0.652
params_v05@species_params$erepro
[1] 9.485557e-01 6.971840e-01 7.477216e-01 6.816138e-01 4.110374e-01 3.877374e-01 7.186214e-01 8.940243e-01
[9] 7.432162e-01 5.455491e-01 4.207180e-01 7.685176e-01 4.934246e-01 4.891521e-01 1.642704e-01 1.706412e-06
params_v05@interaction
prey
predator euphausiids mesopelagic fishes bathypelagic fishes shelf and coastal fishes flying birds
euphausiids 1.0000000 0.1250 0.1250 0.137362637 0.050000000
mesopelagic fishes 0.1250000 1.0000 1.0000 0.072800000 0.012500000
bathypelagic fishes 0.1250000 1.0000 1.0000 0.072800000 0.012500000
shelf and coastal fishes 0.1373626 0.0728 0.0728 1.000000000 0.006868132
flying birds 0.0500000 0.0125 0.0125 0.006868132 0.000000000
small divers 0.1000000 0.0250 0.0250 0.013736264 0.000000000
squids 0.2500000 0.5000 0.5000 0.072800000 0.012500000
toothfishes 0.2500000 0.5000 0.5000 0.072800000 0.012500000
leopard seals 0.1000000 0.0250 0.0250 0.013736264 0.000000000
medium divers 0.1666667 0.3750 0.3750 0.048533333 0.000000000
large divers 0.1666667 0.3750 0.3750 0.048533333 0.000000000
minke whales 0.0500000 0.0125 0.0125 0.006868132 0.000000000
orca 0.2000000 0.0500 0.0500 0.027472527 0.000000000
sperm whales 0.1250000 0.5000 0.5000 0.036400000 0.000000000
baleen whales 0.2000000 0.0500 0.0500 0.027472527 0.000000000
salps 1.0000000 1.0000 1.0000 1.000000000 1.000000000
prey
predator small divers squids toothfishes leopard seals medium divers large divers minke whales
euphausiids 0.10000000 0.2500 0.2500 0.10000000 0.16666667 0.16666667 0.050000000
mesopelagic fishes 0.02500000 0.5000 0.5000 0.02500000 0.37500000 0.37500000 0.012500000
bathypelagic fishes 0.02500000 0.5000 0.5000 0.02500000 0.37500000 0.37500000 0.012500000
shelf and coastal fishes 0.01373626 0.0728 0.0728 0.01373626 0.04853333 0.04853333 0.006868132
flying birds 0.00000000 0.0125 0.0125 0.00000000 0.00000000 0.00000000 0.000000000
small divers 0.00000000 0.0250 0.0250 1.00000000 0.00000000 0.00000000 0.000000000
squids 0.02500000 1.0000 1.0000 0.02500000 0.37500000 0.37500000 0.012500000
toothfishes 0.02500000 1.0000 1.0000 0.02500000 0.37500000 0.37500000 0.012500000
leopard seals 1.00000000 0.0250 0.0250 1.00000000 0.26666667 0.26666667 1.000000000
medium divers 0.00000000 0.3750 0.3750 0.26666667 0.00000000 0.00000000 0.000000000
large divers 0.00000000 0.3750 0.3750 0.26666667 0.00000000 0.00000000 0.000000000
minke whales 0.00000000 0.0125 0.0125 1.00000000 0.00000000 0.00000000 0.000000000
orca 1.00000000 0.0500 0.0500 1.00000000 0.53333333 0.53333333 1.000000000
sperm whales 0.00000000 0.5000 0.5000 0.20000000 0.00000000 0.00000000 0.000000000
baleen whales 0.00000000 0.0500 0.0500 1.00000000 0.00000000 0.00000000 0.000000000
salps 1.00000000 1.0000 1.0000 1.00000000 1.00000000 1.00000000 1.000000000
prey
predator orca sperm whales baleen whales salps
euphausiids 0.20000000 0.1250 0.20000000 1
mesopelagic fishes 0.05000000 0.5000 0.05000000 1
bathypelagic fishes 0.05000000 0.5000 0.05000000 1
shelf and coastal fishes 0.02747253 0.0364 0.02747253 1
flying birds 0.00000000 0.0000 0.00000000 1
small divers 1.00000000 0.0000 0.00000000 1
squids 0.05000000 0.5000 0.05000000 1
toothfishes 0.05000000 0.5000 0.05000000 1
leopard seals 1.00000000 0.2000 1.00000000 1
medium divers 0.53333333 0.0000 0.00000000 1
large divers 0.53333333 0.0000 0.00000000 1
minke whales 1.00000000 0.0000 0.00000000 1
orca 0.00000000 0.4000 1.00000000 1
sperm whales 0.40000000 0.0000 0.00000000 1
baleen whales 1.00000000 0.0000 0.00000000 1
salps 1.00000000 1.0000 1.00000000 1so_theta[16,] # automated salps interaction when added using addSpecies()
so_theta[16,] <- salp_row
so_theta[,16] # automated salps interaction when added using addSpecies()
so_theta[,16] <- salp_row
so_thetaparams_v06 <- setParams(params_v05, interaction = so_theta)
params_v06@interaction
initialN(params_v06) <- sim_v07@n[dim(sim_v07@n)[1],,]sim_v08 <- project(params_v06, t_max = 100)
[=========>----------------------------------] 23% ETA: 1s
[=========>----------------------------------] 24% ETA: 1s
[==========>---------------------------------] 25% ETA: 1s
[==========>---------------------------------] 26% ETA: 1s
[===========>--------------------------------] 27% ETA: 1s
[===========>--------------------------------] 28% ETA: 1s
[============>-------------------------------] 29% ETA: 1s
[============>-------------------------------] 30% ETA: 1s
[=============>------------------------------] 31% ETA: 1s
[=============>------------------------------] 32% ETA: 1s
[=============>------------------------------] 33% ETA: 1s
[==============>-----------------------------] 34% ETA: 1s
[==============>-----------------------------] 35% ETA: 1s
[===============>----------------------------] 36% ETA: 1s
[===============>----------------------------] 37% ETA: 1s
[================>---------------------------] 38% ETA: 1s
[================>---------------------------] 39% ETA: 1s
[================>---------------------------] 40% ETA: 1s
[=================>--------------------------] 41% ETA: 1s
[=================>--------------------------] 42% ETA: 1s
[==================>-------------------------] 43% ETA: 1s
[==================>-------------------------] 44% ETA: 1s
[===================>------------------------] 45% ETA: 1s
[===================>------------------------] 46% ETA: 1s
[===================>------------------------] 47% ETA: 1s
[====================>-----------------------] 48% ETA: 1s
[====================>-----------------------] 49% ETA: 1s
[=====================>----------------------] 50% ETA: 1s
[======================>---------------------] 51% ETA: 1s
[======================>---------------------] 52% ETA: 1s
[=======================>--------------------] 53% ETA: 1s
[=======================>--------------------] 54% ETA: 1s
[=======================>--------------------] 55% ETA: 0s
[========================>-------------------] 56% ETA: 0s
[========================>-------------------] 57% ETA: 0s
[=========================>------------------] 58% ETA: 0s
[=========================>------------------] 59% ETA: 0s
[==========================>-----------------] 60% ETA: 0s
[==========================>-----------------] 61% ETA: 0s
[==========================>-----------------] 62% ETA: 0s
[===========================>----------------] 63% ETA: 0s
[===========================>----------------] 64% ETA: 0s
[============================>---------------] 65% ETA: 0s
[============================>---------------] 66% ETA: 0s
[=============================>--------------] 67% ETA: 0s
[=============================>--------------] 68% ETA: 0s
[=============================>--------------] 69% ETA: 0s
[==============================>-------------] 70% ETA: 0s
[==============================>-------------] 71% ETA: 0s
[===============================>------------] 72% ETA: 0s
[===============================>------------] 73% ETA: 0s
[================================>-----------] 74% ETA: 0s
[================================>-----------] 75% ETA: 0s
[=================================>----------] 76% ETA: 0s
[=================================>----------] 77% ETA: 0s
[=================================>----------] 78% ETA: 0s
[==================================>---------] 79% ETA: 0s
[==================================>---------] 80% ETA: 0s
[===================================>--------] 81% ETA: 0s
[===================================>--------] 82% ETA: 0s
[====================================>-------] 83% ETA: 0s
[====================================>-------] 84% ETA: 0s
[====================================>-------] 85% ETA: 0s
[=====================================>------] 86% ETA: 0s
[=====================================>------] 87% ETA: 0s
[======================================>-----] 88% ETA: 0s
[======================================>-----] 89% ETA: 0s
[=======================================>----] 90% ETA: 0s
[=======================================>----] 91% ETA: 0s
[========================================>---] 92% ETA: 0s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[===========================================>] 99% ETA: 0splotBiomass(sim_v08)
plotDiet(params_v06)
Fill in missing and adjust beta values for zooplankton groups using Heneghan et al. 2020 (10.1016/j.ecolmodel.2020.109265)
Need a value for microzooplankton microzooplankton (in McCormack et al. 2020) is composed of Heterotrophic dinoflagellates, tintinnids, ciliates, copepod nauplii Heneghan et al. log10PPMR values for: Hetero.Flagellates = 0.2–0.72 -> 0.46 Hetero.Ciliates = 2.5–2.9 -> 2.7 Mean: (2.7+0.46)/2 = 1.58 10^1.58 = 38.01894
Need to adjust values for mesozoo, other macrozoo, euphausiids, salps
Heneghan et al. log10PPMR values (midpoints for range) for: salps = 6.8–11.7 -> 9.25 10^9.25 = 1778279410
euphausiids = 6.6–7.8 -> 7.2 10^7.2 = 15848932
mesozoo (copepods) Omni.Cop. = 3.6–4.6 -> 4.1 Carn.Cop. = 0.8–1.9 -> 1.35 Mean: (4.1+1.35)/2 = 2.725 10^2.725 = 530.8844
other macrozoo () Chaetognaths = 1.9–3.4 -> 2.65 10^2.65 = 446.6836
params_tuned_v03 <- tuneParams(params_v06)
params_tuned_v04 <- tuneParams(params_tuned_v03) # updated salps PPMR to match empirical estimates from Heneghan et al 2020getErrorCustom_v3 <- function(vary, params, dat, tol = 0.001,
timetorun = 10)
{
params@species_params$R_max[1:16]<-10^vary[1:16] # R_max for 15 species
params@species_params$erepro[1:16]<-vary[17:32] # erepro for 15 species
params@species_params$interaction_resource[1:16] <- vary[33:48] # interaction_resource for 15 species
params <- setParams(params)
# interaction <- params@interaction
# interaction[] <- matrix(vary[28:108],nrow = 9) # stop at 54 if looking only at 3 biggest species
# params <- setInteraction(params,interaction)
params <- projectToSteady(params, distance_func = distanceSSLogN,
tol = tol, t_max = 200, return_sim = F)
sim <- project(params, t_max = timetorun, progress_bar = F)
sim_biomass = rep(0, length(params@species_params$species))
cutoffLow <- params@species_params$biomass_cutoffLow
if (is.null(cutoffLow))
cutoffLow <- rep(0, no_sp)
cutoffLow[is.na(cutoffLow)] <- 0
cutoffHigh <- params@species_params$biomass_cutoffHigh
if (is.null(cutoffHigh))
cutoffHigh <- rep(0, no_sp)
cutoffHigh[is.na(cutoffHigh)] <- 0
for (j in 1:length(sim_biomass)) {
sim_biomass[j] = sum((sim@n[dim(sim@n)[1],j,] * params@w *
params@dw)[params@w >= cutoffLow[j] & cutoffHigh[j] >= params@w])
}
pred <- log(sim_biomass)
dat <- log(dat)
discrep <- pred - dat
discrep <- (sum(discrep^2))
return(discrep)
}# create set of params for the optimisation process
tic()
params_optim <- params_tuned_v04
vary <- c(log10(params_optim@species_params$R_max),
params_optim@species_params$erepro,
params_optim@species_params$interaction_resource)
params_optim<-setParams(params_optim)
# set up workers
noCores <- parallel::detectCores() - 1 # keep some spare core
cl <- parallel::makeCluster(noCores, setup_timeout = 0.5)
setDefaultCluster(cl = cl)
clusterExport(cl, varlist = "cl",envir=environment())
clusterEvalQ(cl, {
library(mizerExperimental)
library(optimParallel)
})
optim_result <- optimParallel::optimParallel(par=vary,getErrorCustom_v3,params=params_optim,
dat = params_optim@species_params$biomass_observed, method ="L-BFGS-B",
lower=c(rep(-15,16),rep(1e-7,16),rep(.1,16)),
upper= c(rep(15,16),rep(0.99,16),rep(.99,16)),
parallel=list(loginfo=TRUE, forward=TRUE))
stopCluster(cl)
toc()
saveRDS(optim_result, file="optim_result_v03.RDS")
# optim_result <- readRDS("optim_result_v00.RDS")#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max<-10^optim_result$par[1:16]
species_params(params_optim)$erepro<-optim_result$par[17:32]
species_params(params_optim)$interaction_resource <-optim_result$par[33:48]
species_params(params_optim)$R_max [1] 4.830654e+00 2.852640e-03 2.269715e-03 1.109938e-04 1.434663e-07 2.904891e-06 7.470823e-07 6.594487e-07
[9] 3.846612e-10 1.148052e-08 1.182556e-09 1.143978e-09 7.843004e-12 2.415528e-11 1.420689e-10 2.223152e-02species_params(params_optim)$erepro [1] 0.9485563 0.6971843 0.7477218 0.6816140 0.4121015 0.3240383 0.7186214 0.8940243 0.7998257 0.5462956 0.5012431
[12] 0.8250529 0.4934160 0.7479503 0.1740849 0.5808559species_params(params_optim)$interaction_resource [1] 0.9867241 0.9633324 0.9838021 0.9853247 0.8239974 0.8550857 0.9879850 0.9871328 0.9820318 0.9728534 0.6259341
[12] 0.9829386 0.9789605 0.3306710 0.7138859 0.9896762sim_v09 <- project(params_optim, t_max = 200)
[===>----------------------------------------] 9% ETA: 2s
[===>----------------------------------------] 10% ETA: 2s
[====>---------------------------------------] 10% ETA: 2s
[====>---------------------------------------] 11% ETA: 2s
[====>---------------------------------------] 12% ETA: 2s
[=====>--------------------------------------] 13% ETA: 2s
[=====>--------------------------------------] 14% ETA: 2s
[======>-------------------------------------] 15% ETA: 2s
[======>-------------------------------------] 16% ETA: 2s
[======>-------------------------------------] 17% ETA: 2s
[=======>------------------------------------] 17% ETA: 2s
[=======>------------------------------------] 18% ETA: 2s
[=======>------------------------------------] 19% ETA: 2s
[========>-----------------------------------] 19% ETA: 2s
[========>-----------------------------------] 20% ETA: 2s
[========>-----------------------------------] 21% ETA: 2s
[=========>----------------------------------] 22% ETA: 2s
[=========>----------------------------------] 23% ETA: 2s
[==========>---------------------------------] 24% ETA: 2s
[==========>---------------------------------] 25% ETA: 2sWarning in list(...) :
closing unused connection 32 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 31 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 30 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 29 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 28 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 27 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 26 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 25 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 24 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 23 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 22 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 21 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 20 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 19 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 18 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 17 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 16 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 15 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 14 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 13 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 12 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 11 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 10 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 9 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 8 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 7 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 6 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 5 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 4 (<-IMA-1KNH6J3.utas.ad.internal:11624)
Warning in list(...) :
closing unused connection 3 (<-IMA-1KNH6J3.utas.ad.internal:11624)
[==========>---------------------------------] 26% ETA: 2s
[===========>--------------------------------] 26% ETA: 2s
[===========>--------------------------------] 27% ETA: 2s
[===========>--------------------------------] 28% ETA: 2s
[============>-------------------------------] 29% ETA: 2s
[============>-------------------------------] 30% ETA: 2s
[=============>------------------------------] 31% ETA: 2s
[=============>------------------------------] 32% ETA: 2s
[=============>------------------------------] 33% ETA: 2s
[==============>-----------------------------] 33% ETA: 2s
[==============>-----------------------------] 34% ETA: 2s
[==============>-----------------------------] 35% ETA: 2s
[===============>----------------------------] 35% ETA: 2s
[===============>----------------------------] 36% ETA: 2s
[===============>----------------------------] 37% ETA: 2s
[================>---------------------------] 38% ETA: 2s
[================>---------------------------] 39% ETA: 2s
[=================>--------------------------] 40% ETA: 2s
[=================>--------------------------] 41% ETA: 2s
[=================>--------------------------] 42% ETA: 2s
[==================>-------------------------] 42% ETA: 2s
[==================>-------------------------] 43% ETA: 2s
[==================>-------------------------] 44% ETA: 2s
[===================>------------------------] 45% ETA: 2s
[===================>------------------------] 46% ETA: 2s
[====================>-----------------------] 47% ETA: 2s
[====================>-----------------------] 48% ETA: 2s
[====================>-----------------------] 48% ETA: 1s
[====================>-----------------------] 49% ETA: 1s
[=====================>----------------------] 49% ETA: 1s
[=====================>----------------------] 50% ETA: 1s
[=====================>----------------------] 51% ETA: 1s
[======================>---------------------] 51% ETA: 1s
[======================>---------------------] 52% ETA: 1s
[======================>---------------------] 53% ETA: 1s
[=======================>--------------------] 54% ETA: 1s
[=======================>--------------------] 55% ETA: 1s
[========================>-------------------] 56% ETA: 1s
[========================>-------------------] 57% ETA: 1s
[========================>-------------------] 58% ETA: 1s
[=========================>------------------] 58% ETA: 1s
[=========================>------------------] 59% ETA: 1s
[=========================>------------------] 60% ETA: 1s
[==========================>-----------------] 61% ETA: 1s
[==========================>-----------------] 62% ETA: 1s
[===========================>----------------] 63% ETA: 1s
[===========================>----------------] 64% ETA: 1s
[===========================>----------------] 65% ETA: 1s
[============================>---------------] 65% ETA: 1s
[============================>---------------] 66% ETA: 1s
[============================>---------------] 67% ETA: 1s
[=============================>--------------] 67% ETA: 1s
[=============================>--------------] 68% ETA: 1s
[=============================>--------------] 69% ETA: 1s
[==============================>-------------] 70% ETA: 1s
[==============================>-------------] 71% ETA: 1s
[===============================>------------] 72% ETA: 1s
[===============================>------------] 73% ETA: 1s
[===============================>------------] 74% ETA: 1s
[================================>-----------] 74% ETA: 1s
[================================>-----------] 75% ETA: 1s
[================================>-----------] 76% ETA: 1s
[=================================>----------] 77% ETA: 1s
[=================================>----------] 78% ETA: 1s
[==================================>---------] 79% ETA: 1s
[==================================>---------] 80% ETA: 1s
[==================================>---------] 81% ETA: 1s
[===================================>--------] 81% ETA: 1s
[===================================>--------] 82% ETA: 1s
[===================================>--------] 83% ETA: 1s
[====================================>-------] 83% ETA: 0s
[====================================>-------] 84% ETA: 0s
[====================================>-------] 85% ETA: 0s
[=====================================>------] 86% ETA: 0s
[=====================================>------] 87% ETA: 0s
[======================================>-----] 88% ETA: 0s
[======================================>-----] 89% ETA: 0s
[======================================>-----] 90% ETA: 0s
[=======================================>----] 90% ETA: 0s
[=======================================>----] 91% ETA: 0s
[=======================================>----] 92% ETA: 0s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[=========================================>--] 97% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_v09)
params_v07 <- params_optim
initialN(params_v07) <- sim_v09@n[dim(sim_v09@n)[1],,]sim_v10 <- project(params_v07, t_max = 200)
[====>---------------------------------------] 10% ETA: 4s
[====>---------------------------------------] 11% ETA: 4s
[====>---------------------------------------] 12% ETA: 3s
[=====>--------------------------------------] 13% ETA: 3s
[=====>--------------------------------------] 14% ETA: 3s
[======>-------------------------------------] 15% ETA: 3s
[======>-------------------------------------] 16% ETA: 3s
[======>-------------------------------------] 17% ETA: 3s
[=======>------------------------------------] 17% ETA: 3s
[=======>------------------------------------] 18% ETA: 3s
[=======>------------------------------------] 19% ETA: 3s
[========>-----------------------------------] 19% ETA: 3s
[========>-----------------------------------] 20% ETA: 3s
[========>-----------------------------------] 21% ETA: 3s
[========>-----------------------------------] 21% ETA: 2s
[=========>----------------------------------] 22% ETA: 2s
[=========>----------------------------------] 23% ETA: 2s
[==========>---------------------------------] 24% ETA: 2s
[==========>---------------------------------] 25% ETA: 2s
[==========>---------------------------------] 26% ETA: 2s
[===========>--------------------------------] 26% ETA: 2s
[===========>--------------------------------] 27% ETA: 2s
[===========>--------------------------------] 28% ETA: 2s
[============>-------------------------------] 29% ETA: 2s
[============>-------------------------------] 30% ETA: 2s
[=============>------------------------------] 31% ETA: 2s
[=============>------------------------------] 32% ETA: 2s
[=============>------------------------------] 33% ETA: 2s
[==============>-----------------------------] 33% ETA: 2s
[==============>-----------------------------] 34% ETA: 2s
[==============>-----------------------------] 35% ETA: 2s
[===============>----------------------------] 35% ETA: 2s
[===============>----------------------------] 36% ETA: 2s
[===============>----------------------------] 37% ETA: 2s
[================>---------------------------] 38% ETA: 2s
[================>---------------------------] 39% ETA: 2s
[=================>--------------------------] 40% ETA: 2s
[=================>--------------------------] 41% ETA: 2s
[=================>--------------------------] 42% ETA: 2s
[==================>-------------------------] 42% ETA: 2s
[==================>-------------------------] 43% ETA: 2s
[==================>-------------------------] 44% ETA: 2s
[===================>------------------------] 45% ETA: 2s
[===================>------------------------] 46% ETA: 2s
[====================>-----------------------] 47% ETA: 1s
[====================>-----------------------] 48% ETA: 1s
[====================>-----------------------] 49% ETA: 1s
[=====================>----------------------] 49% ETA: 1s
[=====================>----------------------] 50% ETA: 1s
[=====================>----------------------] 51% ETA: 1s
[======================>---------------------] 51% ETA: 1s
[======================>---------------------] 52% ETA: 1s
[======================>---------------------] 53% ETA: 1s
[=======================>--------------------] 54% ETA: 1s
[=======================>--------------------] 55% ETA: 1s
[========================>-------------------] 56% ETA: 1s
[========================>-------------------] 57% ETA: 1s
[========================>-------------------] 58% ETA: 1s
[=========================>------------------] 58% ETA: 1s
[=========================>------------------] 59% ETA: 1s
[=========================>------------------] 60% ETA: 1s
[==========================>-----------------] 61% ETA: 1s
[==========================>-----------------] 62% ETA: 1s
[===========================>----------------] 63% ETA: 1s
[===========================>----------------] 64% ETA: 1s
[===========================>----------------] 65% ETA: 1s
[============================>---------------] 65% ETA: 1s
[============================>---------------] 66% ETA: 1s
[============================>---------------] 67% ETA: 1s
[=============================>--------------] 67% ETA: 1s
[=============================>--------------] 68% ETA: 1s
[=============================>--------------] 69% ETA: 1s
[==============================>-------------] 70% ETA: 1s
[==============================>-------------] 71% ETA: 1s
[===============================>------------] 72% ETA: 1s
[===============================>------------] 73% ETA: 1s
[===============================>------------] 74% ETA: 1s
[================================>-----------] 74% ETA: 1s
[================================>-----------] 75% ETA: 1s
[================================>-----------] 76% ETA: 1s
[=================================>----------] 77% ETA: 1s
[=================================>----------] 78% ETA: 1s
[==================================>---------] 79% ETA: 1s
[==================================>---------] 80% ETA: 1s
[==================================>---------] 81% ETA: 1s
[===================================>--------] 81% ETA: 1s
[===================================>--------] 82% ETA: 1s
[===================================>--------] 83% ETA: 0s
[====================================>-------] 83% ETA: 0s
[====================================>-------] 84% ETA: 0s
[====================================>-------] 85% ETA: 0s
[=====================================>------] 86% ETA: 0s
[=====================================>------] 87% ETA: 0s
[======================================>-----] 88% ETA: 0s
[======================================>-----] 89% ETA: 0s
[======================================>-----] 90% ETA: 0s
[=======================================>----] 90% ETA: 0s
[=======================================>----] 91% ETA: 0s
[=======================================>----] 92% ETA: 0s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[=========================================>--] 97% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_v10)
plotDiet(params_v07)
params_tuned_v05@species_params$erepro
[1] 0.9482281 0.6939526 0.7420393 0.6801803 0.4120497 0.3239043 0.7167263 0.8926139 0.7990909 0.5454373 0.4995779
[12] 0.8216917 0.4943938 0.7415872 0.1725889 0.5806101# saveRDS(params_tuned_v05, "params_optim_v02.rds")
params_tuned_v05 <- readRDS("params_optim_v02.rds")sim_v12 <- project(params_tuned_v05, t_max = 500)
[>-------------------------------------------] 3% ETA: 7s
[=>------------------------------------------] 4% ETA: 7s
[=>------------------------------------------] 5% ETA: 7s
[=>------------------------------------------] 6% ETA: 7s
[==>-----------------------------------------] 6% ETA: 7s
[==>-----------------------------------------] 7% ETA: 7s
[==>-----------------------------------------] 8% ETA: 7s
[===>----------------------------------------] 8% ETA: 6s
[===>----------------------------------------] 8% ETA: 7s
[===>----------------------------------------] 9% ETA: 7s
[===>----------------------------------------] 10% ETA: 7s
[====>---------------------------------------] 10% ETA: 7s
[====>---------------------------------------] 11% ETA: 7s
[====>---------------------------------------] 12% ETA: 7s
[====>---------------------------------------] 12% ETA: 6s
[=====>--------------------------------------] 13% ETA: 6s
[=====>--------------------------------------] 13% ETA: 8s
[=====>--------------------------------------] 14% ETA: 8s
[=====>--------------------------------------] 15% ETA: 8s
[======>-------------------------------------] 15% ETA: 8s
[======>-------------------------------------] 16% ETA: 9s
[======>-------------------------------------] 17% ETA: 9s
[=======>------------------------------------] 17% ETA: 9s
[=======>------------------------------------] 18% ETA: 9s
[=======>------------------------------------] 18% ETA: 8s
[=======>------------------------------------] 19% ETA: 8s
[========>-----------------------------------] 19% ETA: 8s
[========>-----------------------------------] 20% ETA: 8s
[========>-----------------------------------] 21% ETA: 8s
[========>-----------------------------------] 22% ETA: 8s
[=========>----------------------------------] 22% ETA: 8s
[=========>----------------------------------] 23% ETA: 7s
[=========>----------------------------------] 23% ETA: 8s
[=========>----------------------------------] 23% ETA: 7s
[=========>----------------------------------] 24% ETA: 7s
[==========>---------------------------------] 24% ETA: 7s
[==========>---------------------------------] 25% ETA: 7s
[==========>---------------------------------] 26% ETA: 7s
[===========>--------------------------------] 26% ETA: 7s
[===========>--------------------------------] 27% ETA: 7s
[===========>--------------------------------] 28% ETA: 7s
[============>-------------------------------] 29% ETA: 7s
[============>-------------------------------] 30% ETA: 7s
[============>-------------------------------] 30% ETA: 6s
[============>-------------------------------] 31% ETA: 6s
[=============>------------------------------] 31% ETA: 6s
[=============>------------------------------] 32% ETA: 6s
[=============>------------------------------] 33% ETA: 6s
[==============>-----------------------------] 33% ETA: 6s
[==============>-----------------------------] 34% ETA: 6s
[==============>-----------------------------] 35% ETA: 6s
[===============>----------------------------] 35% ETA: 6s
[===============>----------------------------] 36% ETA: 6s
[===============>----------------------------] 37% ETA: 6s
[================>---------------------------] 38% ETA: 6s
[================>---------------------------] 38% ETA: 5s
[================>---------------------------] 39% ETA: 5s
[================>---------------------------] 40% ETA: 5s
[=================>--------------------------] 40% ETA: 5s
[=================>--------------------------] 41% ETA: 5s
[=================>--------------------------] 42% ETA: 5s
[==================>-------------------------] 42% ETA: 5s
[==================>-------------------------] 43% ETA: 5s
[==================>-------------------------] 44% ETA: 5s
[===================>------------------------] 45% ETA: 5s
[===================>------------------------] 46% ETA: 5s
[===================>------------------------] 47% ETA: 5s
[====================>-----------------------] 47% ETA: 5s
[====================>-----------------------] 48% ETA: 5s
[====================>-----------------------] 48% ETA: 4s
[====================>-----------------------] 49% ETA: 4s
[=====================>----------------------] 49% ETA: 4s
[=====================>----------------------] 50% ETA: 4s
[=====================>----------------------] 51% ETA: 4s
[======================>---------------------] 51% ETA: 4s
[======================>---------------------] 52% ETA: 4s
[======================>---------------------] 53% ETA: 4s
[=======================>--------------------] 53% ETA: 4s
[=======================>--------------------] 54% ETA: 4s
[=======================>--------------------] 55% ETA: 4s
[========================>-------------------] 56% ETA: 4s
[========================>-------------------] 57% ETA: 4s
[========================>-------------------] 58% ETA: 4s
[=========================>------------------] 58% ETA: 4s
[=========================>------------------] 58% ETA: 3s
[=========================>------------------] 59% ETA: 3s
[=========================>------------------] 60% ETA: 3s
[==========================>-----------------] 60% ETA: 3s
[==========================>-----------------] 61% ETA: 3s
[==========================>-----------------] 62% ETA: 3s
[===========================>----------------] 63% ETA: 3s
[===========================>----------------] 64% ETA: 3s
[===========================>----------------] 65% ETA: 3s
[============================>---------------] 65% ETA: 3s
[============================>---------------] 66% ETA: 3s
[============================>---------------] 67% ETA: 3s
[=============================>--------------] 67% ETA: 3s
[=============================>--------------] 68% ETA: 3s
[=============================>--------------] 69% ETA: 3s
[==============================>-------------] 69% ETA: 3s
[==============================>-------------] 70% ETA: 3s
[==============================>-------------] 71% ETA: 2s
[===============================>------------] 72% ETA: 2s
[===============================>------------] 73% ETA: 2s
[===============================>------------] 74% ETA: 2s
[================================>-----------] 74% ETA: 2s
[================================>-----------] 75% ETA: 2s
[================================>-----------] 76% ETA: 2s
[=================================>----------] 76% ETA: 2s
[=================================>----------] 77% ETA: 2s
[=================================>----------] 78% ETA: 2s
[==================================>---------] 78% ETA: 2s
[==================================>---------] 79% ETA: 2s
[==================================>---------] 80% ETA: 2s
[==================================>---------] 81% ETA: 2s
[===================================>--------] 81% ETA: 2s
[===================================>--------] 82% ETA: 2s
[===================================>--------] 82% ETA: 1s
[===================================>--------] 83% ETA: 1s
[====================================>-------] 83% ETA: 1s
[====================================>-------] 84% ETA: 1s
[====================================>-------] 85% ETA: 1s
[=====================================>------] 85% ETA: 1s
[=====================================>------] 86% ETA: 1s
[=====================================>------] 87% ETA: 1s
[======================================>-----] 88% ETA: 1s
[======================================>-----] 89% ETA: 1s
[======================================>-----] 90% ETA: 1s
[=======================================>----] 90% ETA: 1s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 1s
[========================================>---] 94% ETA: 1s
[========================================>---] 94% ETA: 0s
[=========================================>--] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0sinitialN(params_tuned_v05) <- sim_v12@n[dim(sim_v12@n)[1],,]
sim_v13 <- project(params_tuned_v05, t_max = 500)
[>-------------------------------------------] 3% ETA: 6s
[=>------------------------------------------] 4% ETA: 6s
[=>------------------------------------------] 5% ETA: 6s
[=>------------------------------------------] 6% ETA: 6s
[==>-----------------------------------------] 6% ETA: 6s
[==>-----------------------------------------] 7% ETA: 6s
[==>-----------------------------------------] 8% ETA: 6s
[===>----------------------------------------] 8% ETA: 6s
[===>----------------------------------------] 9% ETA: 6s
[===>----------------------------------------] 10% ETA: 6s
[====>---------------------------------------] 10% ETA: 6s
[====>---------------------------------------] 11% ETA: 6s
[====>---------------------------------------] 12% ETA: 6s
[=====>--------------------------------------] 13% ETA: 6s
[=====>--------------------------------------] 14% ETA: 6s
[=====>--------------------------------------] 15% ETA: 6s
[======>-------------------------------------] 15% ETA: 6s
[======>-------------------------------------] 16% ETA: 6s
[======>-------------------------------------] 17% ETA: 6s
[=======>------------------------------------] 17% ETA: 6s
[=======>------------------------------------] 18% ETA: 6s
[=======>------------------------------------] 19% ETA: 6s
[========>-----------------------------------] 19% ETA: 6s
[========>-----------------------------------] 20% ETA: 6s
[========>-----------------------------------] 20% ETA: 5s
[========>-----------------------------------] 20% ETA: 6s
[========>-----------------------------------] 21% ETA: 6s
[========>-----------------------------------] 21% ETA: 5s
[========>-----------------------------------] 22% ETA: 5s
[=========>----------------------------------] 22% ETA: 5s
[=========>----------------------------------] 23% ETA: 5s
[=========>----------------------------------] 24% ETA: 5s
[==========>---------------------------------] 24% ETA: 5s
[==========>---------------------------------] 25% ETA: 5s
[==========>---------------------------------] 26% ETA: 5s
[===========>--------------------------------] 26% ETA: 5s
[===========>--------------------------------] 27% ETA: 5s
[===========>--------------------------------] 28% ETA: 5s
[============>-------------------------------] 29% ETA: 5s
[============>-------------------------------] 30% ETA: 5s
[============>-------------------------------] 31% ETA: 5s
[=============>------------------------------] 31% ETA: 5s
[=============>------------------------------] 32% ETA: 5s
[=============>------------------------------] 33% ETA: 5s
[==============>-----------------------------] 33% ETA: 5s
[==============>-----------------------------] 34% ETA: 5s
[==============>-----------------------------] 35% ETA: 5s
[===============>----------------------------] 35% ETA: 4s
[===============>----------------------------] 36% ETA: 4s
[===============>----------------------------] 37% ETA: 4s
[================>---------------------------] 38% ETA: 4s
[================>---------------------------] 39% ETA: 4s
[================>---------------------------] 40% ETA: 4s
[=================>--------------------------] 40% ETA: 4s
[=================>--------------------------] 41% ETA: 4s
[=================>--------------------------] 42% ETA: 4s
[==================>-------------------------] 42% ETA: 4s
[==================>-------------------------] 43% ETA: 4s
[==================>-------------------------] 44% ETA: 4s
[===================>------------------------] 45% ETA: 4s
[===================>------------------------] 46% ETA: 4s
[===================>------------------------] 47% ETA: 4s
[====================>-----------------------] 47% ETA: 4s
[====================>-----------------------] 48% ETA: 4s
[====================>-----------------------] 49% ETA: 4s
[=====================>----------------------] 49% ETA: 4s
[=====================>----------------------] 50% ETA: 4s
[=====================>----------------------] 51% ETA: 4s
[======================>---------------------] 51% ETA: 4s
[======================>---------------------] 52% ETA: 4s
[======================>---------------------] 53% ETA: 4s
[======================>---------------------] 53% ETA: 3s
[=======================>--------------------] 53% ETA: 3s
[=======================>--------------------] 54% ETA: 3s
[=======================>--------------------] 55% ETA: 3s
[========================>-------------------] 56% ETA: 3s
[========================>-------------------] 57% ETA: 3s
[========================>-------------------] 58% ETA: 3s
[=========================>------------------] 58% ETA: 3s
[=========================>------------------] 59% ETA: 3s
[=========================>------------------] 60% ETA: 3s
[==========================>-----------------] 60% ETA: 3s
[==========================>-----------------] 61% ETA: 3s
[==========================>-----------------] 62% ETA: 3s
[===========================>----------------] 63% ETA: 3s
[===========================>----------------] 64% ETA: 3s
[===========================>----------------] 65% ETA: 3s
[============================>---------------] 65% ETA: 3s
[============================>---------------] 66% ETA: 3s
[============================>---------------] 66% ETA: 2s
[============================>---------------] 67% ETA: 2s
[=============================>--------------] 67% ETA: 2s
[=============================>--------------] 68% ETA: 2s
[=============================>--------------] 69% ETA: 2s
[==============================>-------------] 69% ETA: 2s
[==============================>-------------] 70% ETA: 2s
[==============================>-------------] 71% ETA: 2s
[===============================>------------] 72% ETA: 2s
[===============================>------------] 73% ETA: 2s
[===============================>------------] 74% ETA: 2s
[================================>-----------] 74% ETA: 2s
[================================>-----------] 75% ETA: 2s
[================================>-----------] 76% ETA: 2s
[=================================>----------] 76% ETA: 2s
[=================================>----------] 77% ETA: 2s
[=================================>----------] 78% ETA: 2s
[==================================>---------] 78% ETA: 2s
[==================================>---------] 79% ETA: 2s
[==================================>---------] 79% ETA: 1s
[==================================>---------] 80% ETA: 1s
[==================================>---------] 81% ETA: 1s
[===================================>--------] 81% ETA: 1s
[===================================>--------] 82% ETA: 1s
[===================================>--------] 83% ETA: 1s
[====================================>-------] 83% ETA: 1s
[====================================>-------] 84% ETA: 1s
[====================================>-------] 85% ETA: 1s
[=====================================>------] 85% ETA: 1s
[=====================================>------] 86% ETA: 1s
[=====================================>------] 87% ETA: 1s
[======================================>-----] 88% ETA: 1s
[======================================>-----] 89% ETA: 1s
[======================================>-----] 90% ETA: 1s
[=======================================>----] 90% ETA: 1s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 1s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotlyBiomass(sim_v13)plotBiomassObservedVsModel(sim_v13)
plotBiomassObservedVsModelCustom(sim_v13)
Gradually reducing resource maximum size in tuneParams(), as doing it directly in the setParams() route will tell you the w_pp_cutoff has changed, but it will still incorporate a background resource up to the w_pp_cutoff that was originally used in newMultispeciesParams()
params_tuned_v06@species_params$R_max
[1] 4.664858e+00 2.181613e-03 1.899593e-03 7.871131e-05 6.543008e-08 4.569620e-07 4.394343e-07
[8] 7.671631e-07 1.613946e-10 1.375584e-08 Inf 7.111924e-10 1.654097e-11 1.156828e-12
[15] 9.964411e-12 1.183440e-02params_tuned_v10@resource_params$w_pp_cutoff
[1] 100# create set of params for the optimisation process
tic()
params_optim <- params_tuned_v10
vary <- c(log10(params_optim@species_params$R_max),
params_optim@species_params$erepro,
params_optim@species_params$interaction_resource)
params_optim<-setParams(params_optim)
# set up workers
noCores <- parallel::detectCores() - 1 # keep some spare core
cl <- parallel::makeCluster(noCores, setup_timeout = 0.5)
setDefaultCluster(cl = cl)
clusterExport(cl, varlist = "cl",envir=environment())
clusterEvalQ(cl, {
library(mizerExperimental)
library(optimParallel)
})
optim_result <- optimParallel::optimParallel(par=vary,getErrorCustom_v3,params=params_optim,
dat = params_optim@species_params$biomass_observed, method ="L-BFGS-B",
lower=c(rep(-15,16),rep(1e-7,16),rep(.1,16)),
upper= c(rep(15,16),rep(0.99,16),rep(.99,16)),
parallel=list(loginfo=TRUE, forward=TRUE))
stopCluster(cl)
toc()
saveRDS(optim_result, file="optim_result_v04.RDS")
# optim_result <- readRDS("optim_result_v04.RDS")#put these new vals intospecies_params and go back to the top of this page to re-check the calibration
species_params(params_optim)$R_max<-10^optim_result$par[1:16]
species_params(params_optim)$erepro<-optim_result$par[17:32]
species_params(params_optim)$interaction_resource <-optim_result$par[33:48]
species_params(params_optim)$R_max [1] 4.455921e+00 2.984966e-03 2.445172e-03 1.104837e-04 1.274650e-07 2.220551e-06 7.643668e-07
[8] 9.188736e-07 3.643475e-10 1.265181e-08 1.011248e-09 1.083077e-09 8.519112e-12 1.730661e-11
[15] 1.425133e-10 2.233647e-02species_params(params_optim)$erepro [1] 0.9485960 0.6957507 0.7442478 0.6805366 0.4135304 0.4176807 0.7177938 0.8906210 0.8206870 0.5417358
[11] 0.8990635 0.8349824 0.4713801 0.7878068 0.1844293 0.5808197species_params(params_optim)$interaction_resource [1] 0.9900000 0.9757103 0.9900000 0.9898909 0.7710204 0.8767302 0.9900000 0.9900000 0.9532022 0.9867254
[11] 0.6868874 0.9900000 0.9787623 0.3054025 0.6562795 0.9897251sim_v14 <- project(params_optim, t_max = 2000)
[--------------------------------------------] 0% ETA: 1m
[--------------------------------------------] 0% ETA: 2m
[--------------------------------------------] 1% ETA: 2m
[--------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 2% ETA: 1m
[>-------------------------------------------] 3% ETA: 1m
[=>------------------------------------------] 3% ETA: 1m
[=>------------------------------------------] 4% ETA: 1m
[=>------------------------------------------] 5% ETA: 1m
[=>------------------------------------------] 6% ETA: 1m
[==>-----------------------------------------] 6% ETA: 1m
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 49s
[==>-----------------------------------------] 7% ETA: 49s
[==>-----------------------------------------] 7% ETA: 48s
[==>-----------------------------------------] 7% ETA: 47s
[==>-----------------------------------------] 8% ETA: 47s
[==>-----------------------------------------] 8% ETA: 46s
[==>-----------------------------------------] 8% ETA: 47s
[==>-----------------------------------------] 8% ETA: 46s
[==>-----------------------------------------] 8% ETA: 47s
[===>----------------------------------------] 8% ETA: 46s
[===>----------------------------------------] 8% ETA: 45s
[===>----------------------------------------] 9% ETA: 46s
[===>----------------------------------------] 9% ETA: 45s
[===>----------------------------------------] 9% ETA: 44s
[===>----------------------------------------] 9% ETA: 47s
[===>----------------------------------------] 9% ETA: 46s
[===>----------------------------------------] 10% ETA: 47s
[===>----------------------------------------] 10% ETA: 46s
[===>----------------------------------------] 10% ETA: 45s
[===>----------------------------------------] 10% ETA: 46s
[===>----------------------------------------] 10% ETA: 45s
[====>---------------------------------------] 10% ETA: 45s
[====>---------------------------------------] 11% ETA: 45s
[====>---------------------------------------] 11% ETA: 44s
[====>---------------------------------------] 12% ETA: 44s
[====>---------------------------------------] 12% ETA: 43s
[=====>--------------------------------------] 13% ETA: 43s
[=====>--------------------------------------] 13% ETA: 42s
[=====>--------------------------------------] 13% ETA: 43s
[=====>--------------------------------------] 13% ETA: 42s
[=====>--------------------------------------] 13% ETA: 43s
[=====>--------------------------------------] 13% ETA: 42s
[=====>--------------------------------------] 14% ETA: 42s
[=====>--------------------------------------] 14% ETA: 41s
[=====>--------------------------------------] 15% ETA: 41s
[======>-------------------------------------] 15% ETA: 41s
[======>-------------------------------------] 15% ETA: 40s
[======>-------------------------------------] 15% ETA: 41s
[======>-------------------------------------] 15% ETA: 40s
[======>-------------------------------------] 16% ETA: 40s
[======>-------------------------------------] 16% ETA: 39s
[======>-------------------------------------] 16% ETA: 40s
[======>-------------------------------------] 16% ETA: 39s
[======>-------------------------------------] 16% ETA: 40s
[======>-------------------------------------] 17% ETA: 40s
[=======>------------------------------------] 17% ETA: 40s
[=======>------------------------------------] 17% ETA: 39s
[=======>------------------------------------] 17% ETA: 41s
[=======>------------------------------------] 18% ETA: 41s
[=======>------------------------------------] 18% ETA: 40s
[=======>------------------------------------] 19% ETA: 40s
[=======>------------------------------------] 19% ETA: 39s
[========>-----------------------------------] 19% ETA: 39s
[========>-----------------------------------] 20% ETA: 39s
[========>-----------------------------------] 20% ETA: 38s
[========>-----------------------------------] 21% ETA: 38s
[========>-----------------------------------] 21% ETA: 37s
[========>-----------------------------------] 21% ETA: 38s
[========>-----------------------------------] 21% ETA: 37s
[========>-----------------------------------] 22% ETA: 37s
[=========>----------------------------------] 22% ETA: 37s
[=========>----------------------------------] 23% ETA: 37s
[=========>----------------------------------] 23% ETA: 36s
[=========>----------------------------------] 23% ETA: 37s
[=========>----------------------------------] 23% ETA: 36s
[=========>----------------------------------] 24% ETA: 36s
[==========>---------------------------------] 24% ETA: 36s
[==========>---------------------------------] 24% ETA: 35s
[==========>---------------------------------] 25% ETA: 35s
[==========>---------------------------------] 25% ETA: 36s
[==========>---------------------------------] 25% ETA: 35s
[==========>---------------------------------] 26% ETA: 35s
[===========>--------------------------------] 26% ETA: 35s
[===========>--------------------------------] 26% ETA: 34s
[===========>--------------------------------] 27% ETA: 34s
[===========>--------------------------------] 28% ETA: 34s
[===========>--------------------------------] 28% ETA: 33s
[============>-------------------------------] 28% ETA: 33s
[============>-------------------------------] 29% ETA: 33s
[============>-------------------------------] 29% ETA: 32s
[============>-------------------------------] 30% ETA: 32s
[============>-------------------------------] 31% ETA: 32s
[=============>------------------------------] 31% ETA: 32s
[=============>------------------------------] 31% ETA: 31s
[=============>------------------------------] 32% ETA: 31s
[=============>------------------------------] 32% ETA: 30s
[=============>------------------------------] 33% ETA: 30s
[=============>------------------------------] 33% ETA: 31s
[==============>-----------------------------] 33% ETA: 31s
[==============>-----------------------------] 33% ETA: 30s
[==============>-----------------------------] 34% ETA: 30s
[==============>-----------------------------] 35% ETA: 30s
[==============>-----------------------------] 35% ETA: 29s
[===============>----------------------------] 35% ETA: 29s
[===============>----------------------------] 36% ETA: 29s
[===============>----------------------------] 37% ETA: 29s
[===============>----------------------------] 37% ETA: 28s
[===============>----------------------------] 37% ETA: 29s
[===============>----------------------------] 37% ETA: 28s
[================>---------------------------] 38% ETA: 28s
[================>---------------------------] 38% ETA: 27s
[================>---------------------------] 39% ETA: 27s
[================>---------------------------] 40% ETA: 27s
[=================>--------------------------] 40% ETA: 27s
[=================>--------------------------] 40% ETA: 26s
[=================>--------------------------] 40% ETA: 27s
[=================>--------------------------] 41% ETA: 27s
[=================>--------------------------] 41% ETA: 26s
[=================>--------------------------] 42% ETA: 26s
[==================>-------------------------] 42% ETA: 26s
[==================>-------------------------] 43% ETA: 26s
[==================>-------------------------] 43% ETA: 25s
[==================>-------------------------] 43% ETA: 26s
[==================>-------------------------] 43% ETA: 25s
[==================>-------------------------] 44% ETA: 25s
[===================>------------------------] 44% ETA: 25s
[===================>------------------------] 45% ETA: 25s
[===================>------------------------] 45% ETA: 24s
[===================>------------------------] 46% ETA: 24s
[===================>------------------------] 47% ETA: 24s
[====================>-----------------------] 47% ETA: 24s
[====================>-----------------------] 47% ETA: 23s
[====================>-----------------------] 48% ETA: 23s
[====================>-----------------------] 49% ETA: 23s
[=====================>----------------------] 49% ETA: 23s
[=====================>----------------------] 49% ETA: 22s
[=====================>----------------------] 50% ETA: 22s
[=====================>----------------------] 51% ETA: 22s
[======================>---------------------] 51% ETA: 22s
[======================>---------------------] 52% ETA: 21s
[======================>---------------------] 53% ETA: 21s
[=======================>--------------------] 53% ETA: 21s
[=======================>--------------------] 54% ETA: 21s
[=======================>--------------------] 54% ETA: 20s
[=======================>--------------------] 55% ETA: 20s
[=======================>--------------------] 56% ETA: 20s
[========================>-------------------] 56% ETA: 20s
[========================>-------------------] 56% ETA: 19s
[========================>-------------------] 57% ETA: 19s
[========================>-------------------] 58% ETA: 19s
[=========================>------------------] 58% ETA: 19s
[=========================>------------------] 59% ETA: 19s
[=========================>------------------] 59% ETA: 18s
[=========================>------------------] 60% ETA: 18s
[==========================>-----------------] 60% ETA: 18s
[==========================>-----------------] 61% ETA: 18s
[==========================>-----------------] 61% ETA: 17s
[==========================>-----------------] 62% ETA: 17s
[===========================>----------------] 63% ETA: 17s
[===========================>----------------] 64% ETA: 17s
[===========================>----------------] 64% ETA: 16s
[===========================>----------------] 65% ETA: 16s
[============================>---------------] 65% ETA: 16s
[============================>---------------] 66% ETA: 16s
[============================>---------------] 66% ETA: 15s
[============================>---------------] 67% ETA: 15s
[=============================>--------------] 67% ETA: 15s
[=============================>--------------] 68% ETA: 15s
[=============================>--------------] 68% ETA: 14s
[=============================>--------------] 69% ETA: 14s
[==============================>-------------] 69% ETA: 14s
[==============================>-------------] 70% ETA: 14s
[==============================>-------------] 71% ETA: 14s
[==============================>-------------] 71% ETA: 13s
[==============================>-------------] 72% ETA: 13s
[===============================>------------] 72% ETA: 13s
[===============================>------------] 73% ETA: 13s
[===============================>------------] 73% ETA: 12s
[===============================>------------] 74% ETA: 12s
[================================>-----------] 74% ETA: 12s
[================================>-----------] 75% ETA: 12s
[================================>-----------] 75% ETA: 11s
[================================>-----------] 76% ETA: 11s
[=================================>----------] 76% ETA: 11s
[=================================>----------] 77% ETA: 11s
[=================================>----------] 77% ETA: 10s
[=================================>----------] 78% ETA: 10s
[==================================>---------] 78% ETA: 10s
[==================================>---------] 79% ETA: 10s
[==================================>---------] 79% ETA: 9s
[==================================>---------] 80% ETA: 9s
[==================================>---------] 81% ETA: 9s
[===================================>--------] 81% ETA: 9s
[===================================>--------] 82% ETA: 9s
[===================================>--------] 82% ETA: 8s
[===================================>--------] 83% ETA: 8s
[====================================>-------] 83% ETA: 8s
[====================================>-------] 84% ETA: 8s
[====================================>-------] 84% ETA: 7s
[====================================>-------] 85% ETA: 7s
[=====================================>------] 85% ETA: 7s
[=====================================>------] 86% ETA: 7s
[=====================================>------] 86% ETA: 6s
[=====================================>------] 87% ETA: 6s
[======================================>-----] 88% ETA: 6s
[======================================>-----] 88% ETA: 5s
[======================================>-----] 89% ETA: 5s
[======================================>-----] 90% ETA: 5s
[=======================================>----] 90% ETA: 5s
[=======================================>----] 90% ETA: 4s
[=======================================>----] 91% ETA: 4s
[=======================================>----] 92% ETA: 4s
[========================================>---] 92% ETA: 4s
[========================================>---] 93% ETA: 4s
[========================================>---] 93% ETA: 3s
[========================================>---] 94% ETA: 3s
[=========================================>--] 94% ETA: 3s
[=========================================>--] 95% ETA: 3s
[=========================================>--] 95% ETA: 2s
[=========================================>--] 96% ETA: 2s
[=========================================>--] 97% ETA: 2s
[==========================================>-] 97% ETA: 2s
[==========================================>-] 97% ETA: 1s
[==========================================>-] 98% ETA: 1s
[==========================================>-] 99% ETA: 1s
[===========================================>] 99% ETA: 1s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_v14)
params_v08 <- params_optim
initialN(params_v08) <- sim_v14@n[dim(sim_v14@n)[1],,]params_tuned_v07@species_params$erepro
[1] 0.9485821 0.6956672 0.7441494 0.6805253 0.4135289 0.4176799 0.7177760 0.8905317 0.8214728 0.5421544
[11] 0.9005723 0.8378706 0.4727959 0.7956823 0.1849946 0.5808061params_tuned_v07@species_params$R_max [1] 4.289003e+00 2.077178e-03 1.798165e-03 7.486774e-05 6.487583e-08 3.340159e-07 4.677235e-07
[8] 5.228974e-07 1.661754e-10 7.504264e-09 3.456007e-10 6.683396e-10 7.048537e-12 3.203121e-13
[15] 5.664588e-12 1.181176e-02sim_v15 <- project(params_tuned_v07, t_max = 2000)
[--------------------------------------------] 0% ETA: 46s
[--------------------------------------------] 0% ETA: 47s
[--------------------------------------------] 1% ETA: 47s
[--------------------------------------------] 1% ETA: 1m
[--------------------------------------------] 1% ETA: 50s
[--------------------------------------------] 1% ETA: 1m
[--------------------------------------------] 1% ETA: 50s
[>-------------------------------------------] 1% ETA: 49s
[>-------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 1% ETA: 50s
[>-------------------------------------------] 1% ETA: 1m
[>-------------------------------------------] 2% ETA: 50s
[>-------------------------------------------] 2% ETA: 49s
[>-------------------------------------------] 2% ETA: 1m
[>-------------------------------------------] 2% ETA: 50s
[>-------------------------------------------] 2% ETA: 1m
[>-------------------------------------------] 2% ETA: 50s
[>-------------------------------------------] 2% ETA: 1m
[>-------------------------------------------] 2% ETA: 50s
[>-------------------------------------------] 2% ETA: 49s
[>-------------------------------------------] 3% ETA: 49s
[>-------------------------------------------] 3% ETA: 1m
[>-------------------------------------------] 3% ETA: 50s
[>-------------------------------------------] 3% ETA: 49s
[>-------------------------------------------] 3% ETA: 1m
[>-------------------------------------------] 3% ETA: 50s
[>-------------------------------------------] 3% ETA: 1m
[>-------------------------------------------] 3% ETA: 50s
[>-------------------------------------------] 3% ETA: 49s
[>-------------------------------------------] 3% ETA: 1m
[>-------------------------------------------] 3% ETA: 50s
[=>------------------------------------------] 3% ETA: 50s
[=>------------------------------------------] 4% ETA: 49s
[=>------------------------------------------] 4% ETA: 50s
[=>------------------------------------------] 4% ETA: 49s
[=>------------------------------------------] 4% ETA: 48s
[=>------------------------------------------] 4% ETA: 49s
[=>------------------------------------------] 5% ETA: 49s
[=>------------------------------------------] 5% ETA: 48s
[=>------------------------------------------] 5% ETA: 49s
[=>------------------------------------------] 5% ETA: 1m
[=>------------------------------------------] 5% ETA: 50s
[=>------------------------------------------] 5% ETA: 1m
[=>------------------------------------------] 5% ETA: 50s
[=>------------------------------------------] 5% ETA: 1m
[=>------------------------------------------] 5% ETA: 50s
[=>------------------------------------------] 5% ETA: 49s
[=>------------------------------------------] 5% ETA: 50s
[=>------------------------------------------] 6% ETA: 49s
[=>------------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 49s
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 49s
[==>-----------------------------------------] 6% ETA: 1m
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 6% ETA: 1m
[==>-----------------------------------------] 6% ETA: 50s
[==>-----------------------------------------] 7% ETA: 50s
[==>-----------------------------------------] 7% ETA: 1m
[==>-----------------------------------------] 7% ETA: 50s
[==>-----------------------------------------] 7% ETA: 1m
[==>-----------------------------------------] 7% ETA: 50s
[==>-----------------------------------------] 7% ETA: 49s
[==>-----------------------------------------] 7% ETA: 50s
[==>-----------------------------------------] 7% ETA: 49s
[==>-----------------------------------------] 8% ETA: 49s
[==>-----------------------------------------] 8% ETA: 48s
[==>-----------------------------------------] 8% ETA: 49s
[==>-----------------------------------------] 8% ETA: 48s
[==>-----------------------------------------] 8% ETA: 1m
[===>----------------------------------------] 8% ETA: 1m
[===>----------------------------------------] 9% ETA: 1m
[===>----------------------------------------] 9% ETA: 50s
[===>----------------------------------------] 10% ETA: 50s
[===>----------------------------------------] 10% ETA: 49s
[===>----------------------------------------] 10% ETA: 50s
[===>----------------------------------------] 10% ETA: 49s
[===>----------------------------------------] 10% ETA: 50s
[===>----------------------------------------] 10% ETA: 49s
[====>---------------------------------------] 10% ETA: 49s
[====>---------------------------------------] 11% ETA: 49s
[====>---------------------------------------] 11% ETA: 48s
[====>---------------------------------------] 12% ETA: 48s
[====>---------------------------------------] 12% ETA: 47s
[====>---------------------------------------] 12% ETA: 48s
[====>---------------------------------------] 12% ETA: 47s
[=====>--------------------------------------] 13% ETA: 47s
[=====>--------------------------------------] 13% ETA: 46s
[=====>--------------------------------------] 14% ETA: 46s
[=====>--------------------------------------] 14% ETA: 45s
[=====>--------------------------------------] 15% ETA: 45s
[======>-------------------------------------] 15% ETA: 45s
[======>-------------------------------------] 15% ETA: 46s
[======>-------------------------------------] 16% ETA: 46s
[======>-------------------------------------] 16% ETA: 45s
[======>-------------------------------------] 17% ETA: 45s
[=======>------------------------------------] 17% ETA: 45s
[=======>------------------------------------] 18% ETA: 45s
[=======>------------------------------------] 18% ETA: 44s
[=======>------------------------------------] 19% ETA: 44s
[=======>------------------------------------] 19% ETA: 43s
[========>-----------------------------------] 19% ETA: 43s
[========>-----------------------------------] 20% ETA: 43s
[========>-----------------------------------] 20% ETA: 42s
[========>-----------------------------------] 21% ETA: 42s
[========>-----------------------------------] 21% ETA: 41s
[========>-----------------------------------] 22% ETA: 41s
[=========>----------------------------------] 22% ETA: 41s
[=========>----------------------------------] 22% ETA: 42s
[=========>----------------------------------] 23% ETA: 42s
[=========>----------------------------------] 23% ETA: 41s
[=========>----------------------------------] 23% ETA: 42s
[=========>----------------------------------] 24% ETA: 42s
[==========>---------------------------------] 24% ETA: 41s
[==========>---------------------------------] 24% ETA: 42s
[==========>---------------------------------] 24% ETA: 41s
[==========>---------------------------------] 24% ETA: 42s
[==========>---------------------------------] 25% ETA: 42s
[==========>---------------------------------] 25% ETA: 41s
[==========>---------------------------------] 26% ETA: 41s
[===========>--------------------------------] 26% ETA: 41s
[===========>--------------------------------] 27% ETA: 40s
[===========>--------------------------------] 28% ETA: 40s
[===========>--------------------------------] 28% ETA: 39s
[============>-------------------------------] 28% ETA: 39s
[============>-------------------------------] 29% ETA: 39s
[============>-------------------------------] 30% ETA: 39s
[============>-------------------------------] 30% ETA: 38s
[============>-------------------------------] 30% ETA: 40s
[============>-------------------------------] 31% ETA: 39s
[=============>------------------------------] 31% ETA: 39s
[=============>------------------------------] 32% ETA: 39s
[=============>------------------------------] 32% ETA: 38s
[=============>------------------------------] 33% ETA: 38s
[==============>-----------------------------] 33% ETA: 38s
[==============>-----------------------------] 34% ETA: 38s
[==============>-----------------------------] 34% ETA: 37s
[==============>-----------------------------] 35% ETA: 37s
[==============>-----------------------------] 35% ETA: 36s
[===============>----------------------------] 35% ETA: 36s
[===============>----------------------------] 36% ETA: 36s
[===============>----------------------------] 36% ETA: 35s
[===============>----------------------------] 37% ETA: 35s
[================>---------------------------] 38% ETA: 35s
[================>---------------------------] 38% ETA: 34s
[================>---------------------------] 38% ETA: 35s
[================>---------------------------] 38% ETA: 34s
[================>---------------------------] 38% ETA: 35s
[================>---------------------------] 39% ETA: 35s
[================>---------------------------] 39% ETA: 34s
[================>---------------------------] 40% ETA: 34s
[=================>--------------------------] 40% ETA: 34s
[=================>--------------------------] 41% ETA: 33s
[=================>--------------------------] 42% ETA: 33s
[=================>--------------------------] 42% ETA: 32s
[==================>-------------------------] 42% ETA: 32s
[==================>-------------------------] 43% ETA: 32s
[==================>-------------------------] 43% ETA: 31s
[==================>-------------------------] 44% ETA: 31s
[===================>------------------------] 44% ETA: 31s
[===================>------------------------] 45% ETA: 31s
[===================>------------------------] 45% ETA: 30s
[===================>------------------------] 46% ETA: 30s
[===================>------------------------] 47% ETA: 30s
[====================>-----------------------] 47% ETA: 30s
[====================>-----------------------] 47% ETA: 29s
[====================>-----------------------] 48% ETA: 29s
[====================>-----------------------] 48% ETA: 28s
[====================>-----------------------] 49% ETA: 28s
[=====================>----------------------] 49% ETA: 28s
[=====================>----------------------] 50% ETA: 28s
[=====================>----------------------] 50% ETA: 27s
[=====================>----------------------] 51% ETA: 27s
[======================>---------------------] 51% ETA: 27s
[======================>---------------------] 52% ETA: 27s
[======================>---------------------] 52% ETA: 26s
[======================>---------------------] 53% ETA: 26s
[======================>---------------------] 53% ETA: 25s
[=======================>--------------------] 53% ETA: 25s
[=======================>--------------------] 54% ETA: 26s
[=======================>--------------------] 54% ETA: 25s
[=======================>--------------------] 55% ETA: 25s
[=======================>--------------------] 55% ETA: 24s
[=======================>--------------------] 56% ETA: 24s
[========================>-------------------] 56% ETA: 24s
[========================>-------------------] 57% ETA: 24s
[========================>-------------------] 57% ETA: 23s
[========================>-------------------] 58% ETA: 23s
[=========================>------------------] 58% ETA: 23s
[=========================>------------------] 59% ETA: 23s
[=========================>------------------] 59% ETA: 22s
[=========================>------------------] 60% ETA: 22s
[==========================>-----------------] 60% ETA: 22s
[==========================>-----------------] 60% ETA: 21s
[==========================>-----------------] 61% ETA: 21s
[==========================>-----------------] 62% ETA: 21s
[==========================>-----------------] 62% ETA: 20s
[===========================>----------------] 63% ETA: 20s
[===========================>----------------] 64% ETA: 20s
[===========================>----------------] 64% ETA: 19s
[===========================>----------------] 65% ETA: 19s
[============================>---------------] 65% ETA: 19s
[============================>---------------] 66% ETA: 19s
[============================>---------------] 66% ETA: 18s
[============================>---------------] 67% ETA: 18s
[=============================>--------------] 67% ETA: 18s
[=============================>--------------] 68% ETA: 18s
[=============================>--------------] 68% ETA: 17s
[=============================>--------------] 68% ETA: 18s
[=============================>--------------] 68% ETA: 17s
[=============================>--------------] 69% ETA: 17s
[==============================>-------------] 69% ETA: 17s
[==============================>-------------] 70% ETA: 17s
[==============================>-------------] 70% ETA: 16s
[==============================>-------------] 71% ETA: 16s
[==============================>-------------] 71% ETA: 15s
[==============================>-------------] 72% ETA: 15s
[===============================>------------] 72% ETA: 15s
[===============================>------------] 73% ETA: 15s
[===============================>------------] 73% ETA: 14s
[===============================>------------] 74% ETA: 14s
[================================>-----------] 74% ETA: 14s
[================================>-----------] 75% ETA: 14s
[================================>-----------] 75% ETA: 13s
[================================>-----------] 76% ETA: 13s
[=================================>----------] 76% ETA: 13s
[=================================>----------] 77% ETA: 13s
[=================================>----------] 77% ETA: 12s
[=================================>----------] 78% ETA: 12s
[==================================>---------] 78% ETA: 12s
[==================================>---------] 79% ETA: 12s
[==================================>---------] 79% ETA: 11s
[==================================>---------] 80% ETA: 11s
[==================================>---------] 81% ETA: 11s
[==================================>---------] 81% ETA: 10s
[===================================>--------] 81% ETA: 10s
[===================================>--------] 82% ETA: 10s
[===================================>--------] 82% ETA: 9s
[===================================>--------] 83% ETA: 9s
[====================================>-------] 83% ETA: 9s
[====================================>-------] 84% ETA: 9s
[====================================>-------] 84% ETA: 8s
[====================================>-------] 85% ETA: 8s
[=====================================>------] 85% ETA: 8s
[=====================================>------] 86% ETA: 8s
[=====================================>------] 86% ETA: 7s
[=====================================>------] 87% ETA: 7s
[======================================>-----] 88% ETA: 7s
[======================================>-----] 88% ETA: 6s
[======================================>-----] 89% ETA: 6s
[======================================>-----] 90% ETA: 6s
[=======================================>----] 90% ETA: 6s
[=======================================>----] 90% ETA: 5s
[=======================================>----] 91% ETA: 5s
[=======================================>----] 92% ETA: 5s
[=======================================>----] 92% ETA: 4s
[========================================>---] 92% ETA: 4s
[========================================>---] 93% ETA: 4s
[========================================>---] 94% ETA: 4s
[========================================>---] 94% ETA: 3s
[=========================================>--] 94% ETA: 3s
[=========================================>--] 95% ETA: 3s
[=========================================>--] 95% ETA: 2s
[=========================================>--] 96% ETA: 2s
[=========================================>--] 97% ETA: 2s
[==========================================>-] 97% ETA: 2s
[==========================================>-] 97% ETA: 1s
[==========================================>-] 98% ETA: 1s
[==========================================>-] 99% ETA: 1s
[===========================================>] 99% ETA: 1s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotBiomass(sim_v15)
plotDiet(params_tuned_v07)
sim_v17 <- project(params_loop, t_max = 500)
[>-------------------------------------------] 3% ETA: 6s
[>-------------------------------------------] 3% ETA: 7s
[=>------------------------------------------] 4% ETA: 6s
[=>------------------------------------------] 5% ETA: 6s
[=>------------------------------------------] 6% ETA: 6s
[==>-----------------------------------------] 6% ETA: 6s
[==>-----------------------------------------] 7% ETA: 6s
[==>-----------------------------------------] 8% ETA: 6s
[===>----------------------------------------] 8% ETA: 6s
[===>----------------------------------------] 9% ETA: 6s
[===>----------------------------------------] 10% ETA: 6s
[====>---------------------------------------] 10% ETA: 6s
[====>---------------------------------------] 11% ETA: 6s
[====>---------------------------------------] 12% ETA: 6s
[=====>--------------------------------------] 13% ETA: 6s
[=====>--------------------------------------] 14% ETA: 6s
[=====>--------------------------------------] 15% ETA: 6s
[======>-------------------------------------] 15% ETA: 5s
[======>-------------------------------------] 16% ETA: 6s
[======>-------------------------------------] 16% ETA: 5s
[======>-------------------------------------] 17% ETA: 5s
[=======>------------------------------------] 17% ETA: 5s
[=======>------------------------------------] 18% ETA: 5s
[=======>------------------------------------] 19% ETA: 6s
[========>-----------------------------------] 19% ETA: 6s
[========>-----------------------------------] 20% ETA: 6s
[========>-----------------------------------] 21% ETA: 6s
[========>-----------------------------------] 22% ETA: 6s
[=========>----------------------------------] 22% ETA: 6s
[=========>----------------------------------] 23% ETA: 6s
[=========>----------------------------------] 24% ETA: 6s
[==========>---------------------------------] 24% ETA: 6s
[==========>---------------------------------] 25% ETA: 6s
[==========>---------------------------------] 26% ETA: 6s
[===========>--------------------------------] 26% ETA: 6s
[===========>--------------------------------] 27% ETA: 5s
[===========>--------------------------------] 28% ETA: 5s
[============>-------------------------------] 29% ETA: 5s
[============>-------------------------------] 30% ETA: 5s
[============>-------------------------------] 31% ETA: 5s
[=============>------------------------------] 31% ETA: 5s
[=============>------------------------------] 32% ETA: 5s
[=============>------------------------------] 33% ETA: 5s
[==============>-----------------------------] 33% ETA: 5s
[==============>-----------------------------] 34% ETA: 5s
[==============>-----------------------------] 35% ETA: 5s
[===============>----------------------------] 35% ETA: 5s
[===============>----------------------------] 36% ETA: 5s
[===============>----------------------------] 37% ETA: 5s
[================>---------------------------] 38% ETA: 4s
[================>---------------------------] 39% ETA: 4s
[================>---------------------------] 40% ETA: 4s
[=================>--------------------------] 40% ETA: 4s
[=================>--------------------------] 41% ETA: 4s
[=================>--------------------------] 42% ETA: 4s
[==================>-------------------------] 42% ETA: 4s
[==================>-------------------------] 43% ETA: 4s
[==================>-------------------------] 44% ETA: 4s
[===================>------------------------] 45% ETA: 4s
[===================>------------------------] 46% ETA: 4s
[===================>------------------------] 47% ETA: 4s
[====================>-----------------------] 47% ETA: 4s
[====================>-----------------------] 48% ETA: 4s
[====================>-----------------------] 49% ETA: 4s
[=====================>----------------------] 49% ETA: 4s
[=====================>----------------------] 50% ETA: 4s
[=====================>----------------------] 50% ETA: 3s
[=====================>----------------------] 51% ETA: 3s
[======================>---------------------] 51% ETA: 3s
[======================>---------------------] 52% ETA: 3s
[======================>---------------------] 53% ETA: 3s
[=======================>--------------------] 53% ETA: 3s
[=======================>--------------------] 54% ETA: 3s
[=======================>--------------------] 55% ETA: 3s
[========================>-------------------] 56% ETA: 3s
[========================>-------------------] 57% ETA: 3s
[========================>-------------------] 58% ETA: 3s
[=========================>------------------] 58% ETA: 3s
[=========================>------------------] 59% ETA: 3s
[=========================>------------------] 60% ETA: 3s
[==========================>-----------------] 60% ETA: 3s
[==========================>-----------------] 61% ETA: 3s
[==========================>-----------------] 62% ETA: 3s
[===========================>----------------] 63% ETA: 3s
[===========================>----------------] 63% ETA: 2s
[===========================>----------------] 64% ETA: 2s
[===========================>----------------] 65% ETA: 2s
[============================>---------------] 65% ETA: 2s
[============================>---------------] 66% ETA: 2s
[============================>---------------] 67% ETA: 2s
[=============================>--------------] 67% ETA: 2s
[=============================>--------------] 68% ETA: 2s
[=============================>--------------] 69% ETA: 2s
[==============================>-------------] 69% ETA: 2s
[==============================>-------------] 70% ETA: 2s
[==============================>-------------] 71% ETA: 2s
[===============================>------------] 72% ETA: 2s
[===============================>------------] 73% ETA: 2s
[===============================>------------] 74% ETA: 2s
[================================>-----------] 74% ETA: 2s
[================================>-----------] 75% ETA: 2s
[================================>-----------] 76% ETA: 2s
[=================================>----------] 76% ETA: 2s
[=================================>----------] 77% ETA: 2s
[=================================>----------] 78% ETA: 2s
[=================================>----------] 78% ETA: 1s
[==================================>---------] 78% ETA: 1s
[==================================>---------] 79% ETA: 1s
[==================================>---------] 80% ETA: 1s
[==================================>---------] 81% ETA: 1s
[===================================>--------] 81% ETA: 1s
[===================================>--------] 82% ETA: 1s
[===================================>--------] 83% ETA: 1s
[====================================>-------] 83% ETA: 1s
[====================================>-------] 84% ETA: 1s
[====================================>-------] 85% ETA: 1s
[=====================================>------] 85% ETA: 1s
[=====================================>------] 86% ETA: 1s
[=====================================>------] 87% ETA: 1s
[======================================>-----] 88% ETA: 1s
[======================================>-----] 89% ETA: 1s
[======================================>-----] 90% ETA: 1s
[=======================================>----] 90% ETA: 1s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 0s
[========================================>---] 94% ETA: 0s
[=========================================>--] 94% ETA: 0s
[=========================================>--] 95% ETA: 0s
[=========================================>--] 96% ETA: 0s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotlyBiomass(sim_v17)plotBiomassObservedVsModel(sim_v17)
plotBiomassObservedVsModelCustom(sim_v17)
plotBiomassRelative(sim_v17)
plotDiet(params_loop)
sim_v18 <- project(params_loop, t_max = 1000)
[>-------------------------------------------] 1% ETA: 14s
[>-------------------------------------------] 2% ETA: 14s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 2% ETA: 16s
[>-------------------------------------------] 2% ETA: 15s
[>-------------------------------------------] 3% ETA: 15s
[=>------------------------------------------] 3% ETA: 15s
[=>------------------------------------------] 4% ETA: 15s
[=>------------------------------------------] 4% ETA: 14s
[=>------------------------------------------] 4% ETA: 15s
[=>------------------------------------------] 5% ETA: 15s
[=>------------------------------------------] 5% ETA: 14s
[=>------------------------------------------] 6% ETA: 15s
[==>-----------------------------------------] 6% ETA: 15s
[==>-----------------------------------------] 6% ETA: 14s
[==>-----------------------------------------] 6% ETA: 15s
[==>-----------------------------------------] 7% ETA: 15s
[==>-----------------------------------------] 7% ETA: 14s
[==>-----------------------------------------] 7% ETA: 15s
[==>-----------------------------------------] 7% ETA: 14s
[==>-----------------------------------------] 8% ETA: 14s
[==>-----------------------------------------] 8% ETA: 15s
[==>-----------------------------------------] 8% ETA: 14s
[===>----------------------------------------] 8% ETA: 14s
[===>----------------------------------------] 8% ETA: 15s
[===>----------------------------------------] 9% ETA: 14s
[===>----------------------------------------] 9% ETA: 15s
[===>----------------------------------------] 9% ETA: 14s
[===>----------------------------------------] 9% ETA: 18s
[===>----------------------------------------] 9% ETA: 17s
[===>----------------------------------------] 10% ETA: 17s
[====>---------------------------------------] 10% ETA: 17s
[====>---------------------------------------] 11% ETA: 17s
[====>---------------------------------------] 11% ETA: 16s
[====>---------------------------------------] 12% ETA: 16s
[====>---------------------------------------] 12% ETA: 15s
[=====>--------------------------------------] 13% ETA: 16s
[=====>--------------------------------------] 13% ETA: 15s
[=====>--------------------------------------] 14% ETA: 15s
[=====>--------------------------------------] 15% ETA: 15s
[======>-------------------------------------] 15% ETA: 15s
[======>-------------------------------------] 16% ETA: 15s
[======>-------------------------------------] 16% ETA: 14s
[======>-------------------------------------] 17% ETA: 14s
[=======>------------------------------------] 17% ETA: 14s
[=======>------------------------------------] 18% ETA: 14s
[=======>------------------------------------] 19% ETA: 14s
[=======>------------------------------------] 19% ETA: 13s
[========>-----------------------------------] 19% ETA: 13s
[========>-----------------------------------] 20% ETA: 13s
[========>-----------------------------------] 21% ETA: 13s
[========>-----------------------------------] 22% ETA: 13s
[=========>----------------------------------] 22% ETA: 13s
[=========>----------------------------------] 23% ETA: 13s
[=========>----------------------------------] 23% ETA: 12s
[=========>----------------------------------] 24% ETA: 12s
[==========>---------------------------------] 24% ETA: 12s
[==========>---------------------------------] 25% ETA: 12s
[==========>---------------------------------] 26% ETA: 12s
[===========>--------------------------------] 26% ETA: 12s
[===========>--------------------------------] 27% ETA: 12s
[===========>--------------------------------] 28% ETA: 12s
[===========>--------------------------------] 28% ETA: 11s
[============>-------------------------------] 28% ETA: 11s
[============>-------------------------------] 29% ETA: 11s
[============>-------------------------------] 30% ETA: 11s
[============>-------------------------------] 31% ETA: 11s
[=============>------------------------------] 31% ETA: 11s
[=============>------------------------------] 32% ETA: 11s
[=============>------------------------------] 33% ETA: 11s
[==============>-----------------------------] 33% ETA: 10s
[==============>-----------------------------] 34% ETA: 10s
[==============>-----------------------------] 35% ETA: 10s
[===============>----------------------------] 35% ETA: 10s
[===============>----------------------------] 36% ETA: 10s
[===============>----------------------------] 37% ETA: 10s
[================>---------------------------] 38% ETA: 10s
[================>---------------------------] 38% ETA: 9s
[================>---------------------------] 39% ETA: 9s
[================>---------------------------] 40% ETA: 9s
[=================>--------------------------] 40% ETA: 9s
[=================>--------------------------] 41% ETA: 9s
[=================>--------------------------] 42% ETA: 9s
[==================>-------------------------] 42% ETA: 9s
[==================>-------------------------] 43% ETA: 9s
[==================>-------------------------] 44% ETA: 9s
[==================>-------------------------] 44% ETA: 8s
[===================>------------------------] 44% ETA: 8s
[===================>------------------------] 45% ETA: 8s
[===================>------------------------] 46% ETA: 8s
[===================>------------------------] 47% ETA: 8s
[====================>-----------------------] 47% ETA: 8s
[====================>-----------------------] 48% ETA: 8s
[====================>-----------------------] 49% ETA: 8s
[=====================>----------------------] 49% ETA: 8s
[=====================>----------------------] 50% ETA: 8s
[=====================>----------------------] 50% ETA: 7s
[=====================>----------------------] 51% ETA: 7s
[======================>---------------------] 51% ETA: 7s
[======================>---------------------] 52% ETA: 7s
[======================>---------------------] 53% ETA: 7s
[=======================>--------------------] 53% ETA: 7s
[=======================>--------------------] 54% ETA: 7s
[=======================>--------------------] 55% ETA: 7s
[=======================>--------------------] 56% ETA: 7s
[========================>-------------------] 56% ETA: 7s
[========================>-------------------] 56% ETA: 6s
[========================>-------------------] 57% ETA: 6s
[========================>-------------------] 58% ETA: 6s
[=========================>------------------] 58% ETA: 6s
[=========================>------------------] 59% ETA: 6s
[=========================>------------------] 60% ETA: 6s
[==========================>-----------------] 60% ETA: 6s
[==========================>-----------------] 61% ETA: 6s
[==========================>-----------------] 62% ETA: 6s
[===========================>----------------] 63% ETA: 6s
[===========================>----------------] 63% ETA: 5s
[===========================>----------------] 64% ETA: 5s
[===========================>----------------] 65% ETA: 5s
[============================>---------------] 65% ETA: 5s
[============================>---------------] 66% ETA: 5s
[============================>---------------] 67% ETA: 5s
[=============================>--------------] 67% ETA: 5s
[=============================>--------------] 68% ETA: 5s
[=============================>--------------] 69% ETA: 5s
[==============================>-------------] 69% ETA: 5s
[==============================>-------------] 70% ETA: 5s
[==============================>-------------] 70% ETA: 4s
[==============================>-------------] 71% ETA: 4s
[==============================>-------------] 72% ETA: 4s
[===============================>------------] 72% ETA: 4s
[===============================>------------] 73% ETA: 4s
[===============================>------------] 74% ETA: 4s
[================================>-----------] 74% ETA: 4s
[================================>-----------] 75% ETA: 4s
[================================>-----------] 76% ETA: 4s
[=================================>----------] 76% ETA: 4s
[=================================>----------] 77% ETA: 4s
[=================================>----------] 77% ETA: 3s
[=================================>----------] 78% ETA: 3s
[==================================>---------] 78% ETA: 3s
[==================================>---------] 79% ETA: 3s
[==================================>---------] 80% ETA: 3s
[==================================>---------] 81% ETA: 3s
[===================================>--------] 81% ETA: 3s
[===================================>--------] 82% ETA: 3s
[===================================>--------] 83% ETA: 3s
[====================================>-------] 83% ETA: 3s
[====================================>-------] 84% ETA: 2s
[====================================>-------] 85% ETA: 2s
[=====================================>------] 85% ETA: 2s
[=====================================>------] 86% ETA: 2s
[=====================================>------] 87% ETA: 2s
[======================================>-----] 88% ETA: 2s
[======================================>-----] 89% ETA: 2s
[======================================>-----] 90% ETA: 2s
[=======================================>----] 90% ETA: 2s
[=======================================>----] 90% ETA: 1s
[=======================================>----] 91% ETA: 1s
[=======================================>----] 92% ETA: 1s
[========================================>---] 92% ETA: 1s
[========================================>---] 93% ETA: 1s
[========================================>---] 94% ETA: 1s
[=========================================>--] 94% ETA: 1s
[=========================================>--] 95% ETA: 1s
[=========================================>--] 96% ETA: 1s
[=========================================>--] 97% ETA: 1s
[==========================================>-] 97% ETA: 1s
[==========================================>-] 97% ETA: 0s
[==========================================>-] 98% ETA: 0s
[==========================================>-] 99% ETA: 0s
[===========================================>] 99% ETA: 0s
[===========================================>] 100% ETA: 0splotlyBiomass(sim_v18)plotBiomassObservedVsModel(sim_v18)
plotBiomassObservedVsModelCustom(sim_v18)
plotBiomassRelative(sim_v18)
plotDiet(params_loop)
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpTdGVwIDE6IHJldmlzZWQgY2F0Y2ggZm9yIDIwMTAtMjAyMA0KU3RlcCAyOiB0dW5lIHJlc2lsaWVuY2Ugb2YgYWxsIHRoZSBzcGVjaWVzIChpLmUuLCBsb25nLWxpdmVkLCBzbG93ZXIgZ3Jvd2luZyBzcGVjaWVzID0gdnVsbmVyYWJsZS4uLmV0Yy4uLikNClN0ZXAgMzogY2FsaWJyYXRlZCBtb2RlbCB3aXRoIGZpc2hpbmcgbW9ydGFsaXR5IHJlcHJlc2VudGVkICh0aW1lLWF2ZXJhZ2VkIGNhdGNoIGNvbXBhcmlzb24sIG1vZGVsIHZzIG9ic2VydmVkKSA9IGJhc2UgbW9kZWwNClN0ZXAgNDogY2hlY2sgYmFzZSBtb2RlbCB3aXRob3V0IGZpc2hpbmcgdG8gZW5zdXJlIGNvLWV4aXN0ZW5jZSBpcyBzdGlsbCBwcmVzZW50LCBvYnNlcnZlZCBiaW9tYXNzIHdvdWxkIG5vdCBiZSBzcGVjaWZpZWQNClN0ZXAgNTogRXhwZXJpbWVudHMgdGhyb3VnaCB0aW1lIC0gZWZmb3J0IGZ1bmN0aW9ucyAtIHVzZXMgdGhlIGJhc2UgbW9kZWwgd2l0aCBubyBmaXNoaW5nIGVmZm9ydCBhcyBhIHN0YXJ0aW5nIHBvaW50DQpTdGVwIDY6IEV4cG9zZSB0aGUgYmFzZSBtb2RlbCB0byBmaXNoaW5nIHRocm91Z2ggdGltZSBhbmQgY29tcGFyZSB0byBvYnNlcnZlZCBjYXRjaCB0aW1lIHNlcmllcyANCiAtIHN0YXJ0IGZyb20gMCBlZmZvcnQgYXQgZmlyc3QgdGltZS1zdGVwIHRvIG1heCBlZmZvcnQgYXBwcm94IGFyb3VuZCBjb2xsYXBzZQ0KIC0gSWYgdGhlIG1vZGVsbGVkIHZlcnN1cyBvYnNlcnZlZCBkb24ndCBtYXRjaCwgbmVlZCB0byBhZGp1c3QgZWZmb3J0L29yIGluaXRpYWwgYmlvbG9naWNhbCBwYXJhbXMgYW5kIGV4cGVyaW1lbnQgdG8gc2VlIHdoYXQgaW1wcm92ZXMNCiAtIGhhbmQgdHVuZSAvIG9wdGltIG9wdGlvbnMNCkFkanVzdCBhbmQgcmVjYWxpYnJhdGUgaWYgbmVlZGVkIChvcHRpbSkNClVzZSBydWxlcyB0byBjb25zdHJhaW4gLSBpLmUuLCANCi0gdGltZS1hdmVyYWdlZCBiaW9tYXNzZXMgd2l0aGluIGNhbGlicmF0aW9uIHBlcmlvZCBtdXN0IGJlICstIHh4ICUNCi0gY2FuIG9wdGltIHRoaW5ncyBsaWtlIGZpc2hpbmctc2l6ZSBzZWxlY3Rpdml0eSwgUFBNUnMgZXRjDQoNCg0KQ2FuIHdlIHJlcHJvZHVjZSB0aGUgdGltZSBzZXJpZXM/DQoNCg0KTG9hZCBsaWJyYXJpZXMgDQpgYGB7cn0NCnJlbW90ZXM6Omluc3RhbGxfZ2l0aHViKCJzaXplc3BlY3RydW0vbWl6ZXJFeHBlcmltZW50YWwiKQ0KbGlicmFyeShtaXplckV4cGVyaW1lbnRhbCkNCiMgcmVtb3Rlczo6aW5zdGFsbF9naXRodWIoInNpemVzcGVjdHJ1bS9taXplck1SIikNCiMgbGlicmFyeShtaXplck1SKQ0KIyBsaWJyYXJ5KG1pemVyKQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KHRpY3RvYykNCmxpYnJhcnkocGFyYWxsZWwpDQojIGxpYnJhcnkocGxvdGx5KQ0KYGBgDQoNCkNhdGNoIGRhdGEgZm9yIGFsbCBpbmRpdnVkYWxzIGFuZCBzdW1tYXJpc2VkIGFzIHRvdGFscyBwZXIgc3BlY2llcyBlYWNoIHllYXINCg0KTWF5YmUgdGhpcyBpc24ndCBleHByZXNzaW5nIGVmZm9ydCwgcmF0aGVyIGZpc2hpbmcgbW9ydGFsaXR5DQoNCg0KYGBge3J9DQojIGRmX2luZF9DUFVFIDwtIHJlYWRSRFMoImluZF9jYXRjaF93ZWlnaHRfQmFuemFyZUJhbmtfMTkzMF8yMDE5X0NQVUUucmRzIikNCiMgZGZfQ1BVRV9rZ19kYXkgPC0gcmVhZFJEUygiY2F0Y2hfdGltZXNlcmllc19CYW56YXJlQmFua18xOTMwXzIwMTlfQ1BVRS5yZHMiKQ0KIyANCiMgZ2xpbXBzZShkZl9pbmRfQ1BVRSkNCiMgZ2xpbXBzZShkZl9DUFVFX2tnX2RheSkNCmBgYA0KDQpZaWVsZCBmb3IgdGhlIHBlcmlvZCB0aGF0IG1hdGNoZXMgRWNvcGF0aCBtb2RlbCwgcG9zdCAyMDAwIGFubnVhbCBhdmVyYWdlDQpgYGB7cn0NCiMgZGZfeWllbGRfMjAwMCA8LSBkZl9DUFVFX2tnX2RheSAlPiUgDQojICAgZmlsdGVyKFllYXIgPiAyMDA5KSAlPiUgDQojICAgZ3JvdXBfYnkoU3BlY2llcykgJT4lIA0KIyAgIHN1bW1hcmlzZSh5aWVsZCA9IHN1bSh0b3RhbF9jYXRjaF9rZykvMTApDQojIA0KIyBkZl95aWVsZF8yMDAwDQpgYGANCg0KDQoNCmVmZm9ydF9kYXlzL21heCB2YWx1ZSBwbG90dGVkIGJ5IHNwZWNpZXMNClRoaXMgbWlnaHQgaW5kaWNhdGUgdGhlIGVmZm9ydCBjdXJ2ZSByZXF1aXJlZCwgZXNwZWNpYWxseSBpZiB0aGVyZSBpcyBhIGRpZmZlcmVuY2UgYmV0d2VlbiBkaWZmZXJlbnQgc3BlY2llcw0KYGBge3J9DQojIGRmX3Bsb3QgPC0gZGZfQ1BVRV9rZ19kYXkgJT4lIA0KIyAgIGZpbHRlcihZZWFyID4gMjAwMCkgJT4lIA0KIyAgIGdyb3VwX2J5KFNwZWNpZXMpICU+JSANCiMgICBtdXRhdGUobWF4X2VmZm9ydCA9IG1heChlZmZvcnRfZGF5cykpICU+JSANCiMgICBncm91cF9ieShTcGVjaWVzLFllYXIpICU+JSANCiMgICBtdXRhdGUoZWZmb3J0X3N0YW5kYXJkID0gZWZmb3J0X2RheXMvbWF4X2VmZm9ydCkNCmBgYA0KDQpSZWxhdGl2ZSBjaGFuZ2UgaW4gY2F0Y2ggcHJlLTE5MzAsIGNvdWxkIGluZm9ybSB0aGUgcmVsYXRpdmUgY2hhbmdlIGluIGVmZm9ydA0KQXNrIENhbWksIGhvdyBjb3VsZCBJIHJlY29uc3RydWN0IHRoZSBlZmZvcnQgdG8gdGhlIGluaXRpYWwgcG9pbnQgb2YgY2F0Y2ggd2hlbiBlZmZvcnQgPSAwLi4uaS5lLiwgYXBwcm94IDE4NTAgDQpgYGB7cn0NCiMgZGZfcGxvdCAlPiUgDQojICAgZ2dwbG90KGFlcyh4ID0gWWVhciwgeSA9IGVmZm9ydF9zdGFuZGFyZCwgY29sb3VyID0gU3BlY2llcykpICsNCiMgICBnZW9tX3Ntb290aCgpICsNCiMgICBnZW9tX2xpbmUoKQ0KIyAgICMgZmFjZXRfd3JhcCh+U3BlY2llcykgDQpgYGANCg0KYGBge3J9DQojIGRmX3Bsb3QgJT4lIA0KIyAgIGdncGxvdChhZXMoeCA9IFllYXIsIHkgPSBDUFVFLCBjb2xvdXIgPSBTcGVjaWVzKSkgKw0KIyAgICMgZ2VvbV9zbW9vdGgoKSArDQojICAgZ2VvbV9saW5lKCkNCiMgICAjIGZhY2V0X3dyYXAoflNwZWNpZXMpIA0KYGBgDQoNCkxvYWQgY2F0Y2ggbGVuZ3RoDQpgYGB7cn0NCmNhdGNoX2xlbmd0aHMgPC0gcmVhZFJEUygiY2F0Y2hfbGVuZ3Rocy5yZHMiKQ0KDQpnbGltcHNlKGNhdGNoX2xlbmd0aHMpDQpgYGANClRoZSB2YWx1ZSBmb3IgZXVwaGF1c2lpZHMgYG9ic195aWVsZGAgaXMgZnJvbSAgYGNhbGlicmF0aW9uX2NhdGNoX2hpc3Rzb2NfMTg1MF8yMDA0X3JlZ2lvbmFsX21vZGVscy5jc3ZgIGZvdW5kIGF0IGh0dHA6Ly9wb3J0YWwuc2YudXRhcy5lZHUuYXUvdGhyZWRkcy9jYXRhbG9nL2dlbS9maXNobWlwL0lTSU1JUDNhL0lucHV0RGF0YS9maXNoaW5nL2hpc3Rzb2MvY2F0YWxvZy5odG1sLiBJdCBpcyB0aGUgYW5udWFsIGF2ZXJhZ2UgeWllbGQgb3ZlciBhIDIyIHllYXIgcGVyaW9kIGZvciB0aGUgUHJ5ZHogQmF5IFJlZ2lvbi4NCg0KVGhlIHZhbHVlcyBvZiBgb2JzX3lpZWxkYCBmb3IgbWlua2Ugd2hhbGVzICgyMTc1NDc2MTM2KSwgb3JjYSAoMTk5MjM3MjkpLCBzcGVybSB3aGFsZXMgKDQwNjIxNzc0NDUpLCBhbmQgYmFsZWVuIHdoYWxlcyAoNTAzNDE2OTYwNTUpIGFyZSB0aGUgYW5udWFsIGF2ZXJhZ2UgeWllbGQgZnJvbSAxOTMwIC0gMjAxOSBmcm9tIElXQyByZWNvcmRzIG9mIHdoYWxpbmcgaW4gdGhlIFByeWR6IEJheSBtb2RlbCBkb21haW4gKDEsNDMzLDAyOCBrbTIpLg0KDQpBZGp1c3QgdGhlIGNhdGNoIHRvIHJlcHJlc2VudCB0aGUgMjAxMCAtIDIwMjAgdGltZSBwZXJpb2QgKyBpbmNvcnBvcmF0ZSB0b290aGZpc2ggY2F0Y2ggaWYgcG9zc2libGUgKGkuZS4sIFN0YWNleSdzIHBhcGVyIHZhbHVlcykNCg0KDQoNCg0KYGBge3J9DQpkZl9pbmRfQ1BVRSA8LSByZWFkUkRTKCJpbmRfY2F0Y2hfd2VpZ2h0X0JhbnphcmVCYW5rXzE5MzBfMjAxOV9DUFVFLnJkcyIpDQpkZl9DUFVFX2tnX2RheSA8LSByZWFkUkRTKCJjYXRjaF90aW1lc2VyaWVzX0JhbnphcmVCYW5rXzE5MzBfMjAxOV9DUFVFLnJkcyIpDQoNCmdsaW1wc2UoZGZfaW5kX0NQVUUpDQpnbGltcHNlKGRmX0NQVUVfa2dfZGF5KQ0KYGBgDQoNCllpZWxkIGZvciB0aGUgcGVyaW9kIHRoYXQgbWF0Y2hlcyBFY29wYXRoIG1vZGVsLCBwb3N0IDIwMDAgYW5udWFsIGF2ZXJhZ2UNCmBgYHtyfQ0KZGZfeWllbGRfMjAwMF8yMDE5IDwtIGRmX0NQVUVfa2dfZGF5ICU+JSANCiAgZmlsdGVyKFllYXIgPiAyMDAwKSAlPiUgDQogIGdyb3VwX2J5KFNwZWNpZXMpICU+JSANCiAgc3VtbWFyaXNlKHlpZWxkID0gc3VtKHRvdGFsX2NhdGNoX2tnKS8xOSkNCg0KZGZfeWllbGRfMjAwMF8yMDE5DQpgYGANCg0KQWRkIHlpZWxkIGluIHRvbm5lcyBwZXIga20yLCBzYW1lIGFzIGcgbTIgYW5kIGl0IGlzIGNvbnNpc3RlbnQgd2l0aCBgYmlvbWFzc19vYnNlcnZlZGAgDQoNCjEuNDc0MzQxZSsxMiBtXjIgZm9yIG1vZGVsIGRvbWFpbg0KDQoxLjQ3NDM0MWUrMTIvMWUrNiA9IDE0NzQzNDEga21eMg0KDQo1MDc1MTEga2cgb2YgbWlua2Ugd2hhbGUgcGVyIHllYXIgZnJvbSAyMDAwIC0gMjAxOQ0KDQo1MDcuNS8xNDc0MzQxID0gMC4wMDAzNDQyMjE2IHRvbm5lcyBNaW5rZSB3aGFsZSB5aWVsZCBwZXIga20yIGZyb20gMjAwMCAtIDIwMTkNCg0KMTU2MTkuMjQga2cgb2YgYmFsZWVuIHdoYWxlIHBlciB5ZWFyIGZyb20gMjAwMCAtIDIwMTkNCg0KMTUuNi8xNDc0MzQxID0gMS4wNTgxZS0wNSB0b25uZXMgb2YgYmFsZWVuIHdoYWxlIHBlciBrbTIgZnJvbSAyMDAwIC0gMjAxOQ0KDQpgYGB7cn0NCm9ic195aWVsZCA8LSBjKDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAuMDAwMzQ0MjIxNiwgMCwgMCwgMS4wNTgxZS0wNSkgIyB0byBhZGQgYSB5aWVsZF9vYnNlcnZlZCBjb2x1bW4gaW4gc3BlY2llc19wYXJhbXMNCg0KYGBgDQoNCg0KYGBge3J9DQpiaW9tYXNzX2N1dG9mZiA8LSBjKDAuMSwgMSwgMSwgMSwgNDAsIDQ1MDAsIDEsIDEsIDIwMDAwMCwgMTAwMDAsIDEzNTAwMCwgNjAwMDAwLCA0OTAwMDAsIDM2NTAwMDAsIDIyNTAwMDApICMgdG8gYWRkIGEgYmlvbWFzc19jdXRvZmYgY29sdW1uIGluIHNwZWNpZXNfcGFyYW1zDQoNCmBgYA0KDQoNClN0ZWFkeS1zdGF0ZSBtb2RlbCB3aXRoIG5vIGZpc2hpbmcgZWZmb3J0IGNhbGlicmF0ZWQgdG8gb2JzZXJ2ZWQgYmlvbWFzcyB2YWx1ZXMgZnJvbSBNY0Nvcm1hY2sgZXQgYWwuIDIwMjAgdGhhdCByZXByZXNlbnQgYW4gYXZlcmFnZSBzdGF0ZSBvZiB0aGUgZm9vZCB3ZWIgZnJvbSAyMDEwLTIwMjAuDQoNCmBgYHtyfQ0KIyBwYXJhbXMgPC0gIHJlYWRSRFMoInN0YWdlMV9zdGVhZHlfdlhYLnJkcyIpICMgSW5jb3JyZWN0IG9yY2EgYW5kIGxlb3BhcmQgc2VhbCB3X21heCB2YWx1ZXMNCnNvX3BhcmFtcyA8LSByZWFkUkRTKCJwYXJhbXNfMTZfTWFyY2hfMjAyMy5yZHMiKSAjIFVwZGF0ZWQgd19tYXggZm9yIG9yY2EgYW5kIGxlb3BhcmQgc2VhbHMNCg0KIyBwYXJhbXMgPC0gc2V0UGFyYW1zKHNldEZpc2hpbmcocGFyYW1zLCBpbml0aWFsX2VmZm9ydCA9IDAuMikpDQoNCmBgYA0KDQoNCmBgYHtyfQ0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJHlpZWxkX29ic2VydmVkIDwtIG9ic195aWVsZA0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJGJpb21hc3NfY3V0b2ZmIDwtIGJpb21hc3NfY3V0b2ZmDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkYmlvbWFzc19jdXRvZmZMb3cgPC0gYmlvbWFzc19jdXRvZmYNCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSRiaW9tYXNzX2N1dG9mZkhpZ2ggPC0gc3BlY2llc19wYXJhbXMocGFyYW1zKSR3X21heA0KDQpzcGVjaWVzX3BhcmFtcyhzb19wYXJhbXMpJHlpZWxkX29ic2VydmVkIDwtIG9ic195aWVsZA0Kc3BlY2llc19wYXJhbXMoc29fcGFyYW1zKSRiaW9tYXNzX2N1dG9mZiA8LSBiaW9tYXNzX2N1dG9mZg0Kc3BlY2llc19wYXJhbXMoc29fcGFyYW1zKSRiaW9tYXNzX2N1dG9mZkxvdyA8LSBiaW9tYXNzX2N1dG9mZg0Kc3BlY2llc19wYXJhbXMoc29fcGFyYW1zKSRiaW9tYXNzX2N1dG9mZkhpZ2ggPC0gc3BlY2llc19wYXJhbXMoc29fcGFyYW1zKSR3X21heA0KDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykgfD4gZHBseXI6OnNlbGVjdChzcGVjaWVzLCB5aWVsZF9vYnNlcnZlZCwgd19tYXQsIHdfbWF4LCBSX21heCkNCnNwZWNpZXNfcGFyYW1zKHNvX3BhcmFtcykgfD4gZHBseXI6OnNlbGVjdChzcGVjaWVzLCB5aWVsZF9vYnNlcnZlZCwgd19tYXQsIHdfbWF4LCBSX21heCkNCmBgYA0KDQoNCg0KYGBge3J9DQojIGdyb3VwcyB8PiBkcGx5cjo6c2VsZWN0KHNwZWNpZXMsIHlpZWxkX29ic2VydmVkLCBiaW9tYXNzX29ic2VydmVkLCBiaW9tYXNzX2N1dG9mZiwgd19taW4pDQpgYGANCg0KDQpgYGB7cn0NCiMgeWllbGRfb2JzZXJ2ZWQgPC0gYygwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLCAwLjAwMDM0NDIyMTYsIDAsIDAsIDEuMDU4MWUtMDUpICMgdG8gYWRkIGEgeWllbGRfb2JzZXJ2ZWQgY29sdW1uIGluIHNwZWNpZXNfcGFyYW1zDQojIA0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJHlpZWxkX29ic2VydmVkIDwtIHlpZWxkX29ic2VydmVkDQpgYGANCg0KYGBge3J9DQojIGJpb21hc3NfY3V0b2ZmIDwtIGMoMC4xLCAxLCAxLCAxLCA0MCwgNDUwMCwgMSwgMSwgMjAwMDAwLCAxMDAwMCwgMTM1MDAwLCA2MDAwMDAsIDQ5MDAwMCwgMzY1MDAwMCwgMjI1MDAwMCkgIyB0byBhZGQgYSBiaW9tYXNzX2N1dG9mZiBjb2x1bW4gaW4gc3BlY2llc19wYXJhbXMNCiMgDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkYmlvbWFzc19jdXRvZmYgPC0gYmlvbWFzc19jdXRvZmYNCmBgYA0KDQpVcGRhdGUgbWF4IGFuZCBtYXQgc2l6ZSBmb3Igb3JjYSB1c2luZyBzZWFsaWZlYmFzZSB2YWx1ZXMgcmF0aGVyIHRoYW4gdGhlIHByZXZpb3VzIHdfbWF4IG9mIDZ0DQpgYGB7cn0NCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSRzcGVjaWVzWzEzXQ0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJHdfbWF4WzEzXQ0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJHdfbWF0WzEzXQ0KIyANCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSR3X21heFsxM10gPC0gMTA2MjgwMzQNCiMgDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXRbMTNdIDwtIDMxOTg4NTUNCiMgDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXhbMTNdDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXRbMTNdDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXQyNVsxM10NCmBgYA0KDQoNClVwZGF0ZSBtYXggYW5kIG1hdCBzaXplIGZvciBsZW9wYXJkIHNlYWwgdXNpbmcgc2VhbGlmZWJhc2UgdmFsdWVzLiBNYXggd2VpZ2h0IG9ic2VydmVkIHJlcG9ydGVkIGFzIDQ1MGtnLiBVc2luZyB0aGUgc2VhbGlmZWJhc2UgdmFsdWVzIGZvciBMLVcgY29udmVyc2lvbiBhbmQgbWF4IGxlbmd0aCByZXN1bHRzIGluIGFuIGVzdGltYXRlZCBtYXggd2VpZ2h0IG9mIDU0NTg3NS4yZywgd2hpY2ggaXMgdG9vIGhpZ2ggYWJvdmUgdGhlIG1heCB3ZWlnaHQgb2JzZXJ2ZWQuDQpgYGB7cn0NCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSRzcGVjaWVzWzldDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXhbOV0NCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSR3X21hdFs5XQ0KIyANCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSR3X21heFs5XSA8LSA0NTAwMDANCiMgIA0KIyBzcGVjaWVzX3BhcmFtcyhwYXJhbXMpJHdfbWF4WzldDQojIHNwZWNpZXNfcGFyYW1zKHBhcmFtcykkd19tYXRbOV0NCiMgc3BlY2llc19wYXJhbXMocGFyYW1zKSR3X21hdDI1WzldDQpgYGANCg0KQWRqdXN0IGB3X21hdGAgdmFsdWVzIHRoYXQgd2VyZSBjaGFuZ2VkIGJ5IGRlZmF1bHQgaW4gYG5ld011bHRpc3BlY2llc1BhcmFtc2ANCmBgYHtyfQ0KIyBwYXJhbXNfdjEgPC0gcGFyYW1zDQojIA0KIyBwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkd19tYXRbcGFyYW1zX3YxQHNwZWNpZXNfcGFyYW1zJHNwZWNpZXMgPT0gImxhcmdlIGRpdmVycyJdICA8LSBwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkd19tYXhbcGFyYW1zX3YxQHNwZWNpZXNfcGFyYW1zJHNwZWNpZXMgPT0gImxhcmdlIGRpdmVycyJdICogMC45DQojIHBhcmFtc192MUBzcGVjaWVzX3BhcmFtcyR3X21hdFtwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAibWlua2Ugd2hhbGVzIl0gIDwtIHBhcmFtc192MUBzcGVjaWVzX3BhcmFtcyR3X21heFtwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAibWlua2Ugd2hhbGVzIl0gKiAwLjkNCiMgcGFyYW1zX3YxQHNwZWNpZXNfcGFyYW1zJHdfbWF0W3BhcmFtc192MUBzcGVjaWVzX3BhcmFtcyRzcGVjaWVzID09ICJvcmNhIl0gIDwtIHBhcmFtc192MUBzcGVjaWVzX3BhcmFtcyR3X21heFtwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAib3JjYSJdICogMC45DQojIHBhcmFtc192MUBzcGVjaWVzX3BhcmFtcyR3X21hdFtwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAic3Blcm0gd2hhbGVzIl0gIDwtIHBhcmFtc192MUBzcGVjaWVzX3BhcmFtcyR3X21heFtwYXJhbXNfdjFAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAic3Blcm0gd2hhbGVzIl0gKiAwLjkNCiMgDQojIHBhcmFtc192MSA8LSBzZXRQYXJhbXMocGFyYW1zX3YxKQ0KYGBgDQoNCg0KYGBge3J9DQojIHBhcmFtc192MiA8LSBwYXJhbXNfdjENCiMgDQojIHBhcmFtc192MkBzcGVjaWVzX3BhcmFtcyR3X21pbltwYXJhbXNfdjJAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAic21hbGwgZGl2ZXJzIl0gIDwtIHBhcmFtc192MkBzcGVjaWVzX3BhcmFtcyR3X21hdFtwYXJhbXNfdjJAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAic21hbGwgZGl2ZXJzIl0gKiAwLjg1DQojIHBhcmFtc192MkBzcGVjaWVzX3BhcmFtcyR3X21pbltwYXJhbXNfdjJAc3BlY2llc19wYXJhbXMkc3BlY2llcyA9PSAibGVvcGFyZCBzZWFscyJdICA8LSBwYXJhbXNfdjJAc3BlY2llc19wYXJhbXMkd19tYXRbcGFyYW1zX3YyQHNwZWNpZXNfcGFyYW1zJHNwZWNpZXMgPT0gImxlb3BhcmQgc2VhbHMiXSAqIDAuODUNCiMgDQojIHBhcmFtc192MiA8LSBzZXRQYXJhbXMocGFyYW1zX3YyKQ0KYGBgDQoNCg0KYGBge3J9DQpwYXJhbXMgPC0gc3RlYWR5KHNvX3BhcmFtcykNCmBgYA0KDQpDaGVjayBnZWFyIHBhcmFtcw0KDQpgYGB7cn0NCmdlYXJfcGFyYW1zKHBhcmFtcykNCmBgYA0KDQpBZGp1c3QgY2F0Y2hhYmlsaXR5IHRvIG9ubHkgc2VsZWN0IGZpc2hpbmcgb24gc3BlY2llcyB3aXRoIGNhdGNoIGRhdGENClRoaXMgaXMgYSBvbmx5IHN0YXJ0aW5nIHBvaW50DQpgYGB7cn0NCmdlYXJfcGFyYW1zKHBhcmFtcykkY2F0Y2hhYmlsaXR5IDwtICBjKDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAsIDAuNDksIDAsIDAsIDAuMDAxKQ0KYGBgDQoNCg0KQ2F0Y2hhYmlsaXR5OiAxIHZhbHVlIA0KRWZmb3J0IGNoYW5nZXMgdGhyb3VnaCB0aW1lICh2b3lhZ2UgZnJlcXVlbmN5IGFuZCBsZW5ndGgpDQoNCmBgYHtyfQ0KZ2Vhcl9wYXJhbXMocGFyYW1zKSRnZWFyIDwtICBjKCJNYWluIiwgIk1haW4iLCAiTWFpbiIsICJNYWluIiwgIk1haW4iLCAiTWFpbiIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJNYWluIiwgIk1haW4iLCAiTWFpbiIsICJNYWluIiwgIk1haW4iLCAiTWFpbiIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJNYWluIiwgIk1haW4iLCAiTWFpbiIpDQoNCmdlYXJfcGFyYW1zKHBhcmFtcykkbDUwIDwtICBjKDEsIDUsIDUsIDUsIDEwLCAxMCwgMTAsIDIwLCAyMCwgMjAsIDIwLCA4NTAsIDYwMCwgMTUwMCwgMjIwMCkgIyB2YWx1ZXMgZm9yIG1pbmtlLCBvcmNhLCBzcGVybSBhbmQgYmFsZWVuIGFyZSBlc3RpbWF0ZWQgYmFzZWQgb2ZmIGBjYXRjaF9sZW5ndGhzYCwgd2hpbGUgYWxsIG90aGVycyBhcmUgcm91Z2ggZ3Vlc3NlcywgcHVyZWx5IGFzIHRoZSBwYXJhbSB3b24ndCB3b3JrIHdpdGggdmFsdWVzIG1pc3NpbmcuDQoNCmdlYXJfcGFyYW1zKHBhcmFtcykkbDI1IDwtICBjKDAuOCwgNCwgNCwgNCwgOCwgOCwgOCwgMTYsIDE2LCAxNiwgMTYsIDgwMCwgNTAwLCAxNDAwLCAxNTAwKQ0KDQpnZWFyX3BhcmFtcyhwYXJhbXMpJHNlbF9mdW5jIDwtICBjKCJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIsICJzaWdtb2lkX2xlbmd0aCIpDQoNCmdlYXJfcGFyYW1zKHBhcmFtcykNCmBgYA0KDQoNCmBgYHtyfQ0KcGFyYW1zX3YwMiA8LSBzZXRQYXJhbXMoc2V0RmlzaGluZyhwYXJhbXMsIGluaXRpYWxfZWZmb3J0ID0gMC4yKSkNCg0KcGFyYW1zX3YwMiA8LSBzdGVhZHkocGFyYW1zX3YwMikNCg0Kc2ltX3YwMSA8LSBwcm9qZWN0KHBhcmFtc192MDIsIHRfbWF4ID0gMTAwKQ0KcGxvdChzaW1fdjAxKQ0KcGxvdGx5QmlvbWFzcyhzaW1fdjAxKQ0Kc2ltX3YwMUBwYXJhbXNAaW5pdGlhbF9lZmZvcnQNCmBgYA0KDQoNCmBgYHtyfQ0KcGxvdFlpZWxkVnNTaXplKHNpbV92MDEsIHNwZWNpZXMgPSAiYmFsZWVuIHdoYWxlcyIsIGNhdGNoID0gY2F0Y2hfbGVuZ3RocywgDQogICAgICAgICAgICAgICAgeF92YXIgPSAiTGVuZ3RoIikNCg0KIyBwbG90WWllbGRWc1NpemUoc2ltX3YxLCBzcGVjaWVzID0gInNwZXJtIHdoYWxlcyIsIGNhdGNoID0gY2F0Y2hfbGVuZ3RocywgDQojICAgICAgICAgICAgICAgICB4X3ZhciA9ICJMZW5ndGgiKQ0KIyANCiMgcGxvdFlpZWxkVnNTaXplKHNpbV92MSwgc3BlY2llcyA9ICJvcmNhIiwgY2F0Y2ggPSBjYXRjaF9sZW5ndGhzLCANCiMgICAgICAgICAgICAgICAgIHhfdmFyID0gIkxlbmd0aCIpDQoNCnBsb3RZaWVsZFZzU2l6ZShzaW1fdjAxLCBzcGVjaWVzID0gIm1pbmtlIHdoYWxlcyIsIGNhdGNoID0gY2F0Y2hfbGVuZ3RocywgDQogICAgICAgICAgICAgICAgeF92YXIgPSAiTGVuZ3RoIikNCmBgYA0KDQpgYGB7cn0NCmdldEVycm9yQ3VzdG9tIDwtIGZ1bmN0aW9uKHZhcnksIHBhcmFtcywgZGF0LCB0b2wgPSAwLjAwMSwgDQogICAgdGltZXRvcnVuID0gMTApDQp7DQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRSX21heFsxOjE1XTwtMTBedmFyeVsxOjE1XSAjIFJfbWF4IGZvciAxNSBzcGVjaWVzDQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRlcmVwcm9bMToxNV08LXZhcnlbMTY6MzBdICMgZXJlcHJvIGZvciAxNSBzcGVjaWVzDQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRpbnRlcmFjdGlvbl9yZXNvdXJjZVsxOjE1XSA8LSB2YXJ5WzMxOjQ1XSAjIGludGVyYWN0aW9uX3Jlc291cmNlIGZvciAxNSBzcGVjaWVzDQogIHBhcmFtcyA8LSBzZXRQYXJhbXMocGFyYW1zKQ0KICAjIGludGVyYWN0aW9uIDwtIHBhcmFtc0BpbnRlcmFjdGlvbg0KICAjIGludGVyYWN0aW9uW10gPC0gbWF0cml4KHZhcnlbMjg6MTA4XSxucm93ID0gOSkgIyBzdG9wIGF0IDU0IGlmIGxvb2tpbmcgb25seSBhdCAzIGJpZ2dlc3Qgc3BlY2llcw0KICANCiAgIyBwYXJhbXMgPC0gc2V0SW50ZXJhY3Rpb24ocGFyYW1zLGludGVyYWN0aW9uKQ0KICAgIHBhcmFtcyA8LSBwcm9qZWN0VG9TdGVhZHkocGFyYW1zLCBkaXN0YW5jZV9mdW5jID0gZGlzdGFuY2VTU0xvZ04sIA0KICAgICAgICB0b2wgPSB0b2wsIHRfbWF4ID0gMjAwLCByZXR1cm5fc2ltID0gRikNCiAgICANCiAgICBzaW0gPC0gcHJvamVjdChwYXJhbXMsIHRfbWF4ID0gdGltZXRvcnVuLCBwcm9ncmVzc19iYXIgPSBGKQ0KICAgIA0KICAgIHNpbV9iaW9tYXNzID0gcmVwKDAsIGxlbmd0aChwYXJhbXNAc3BlY2llc19wYXJhbXMkc3BlY2llcykpDQogICAgDQogICAgICAgIGN1dG9mZkxvdyA8LSBwYXJhbXNAc3BlY2llc19wYXJhbXMkYmlvbWFzc19jdXRvZmZMb3cNCiAgICBpZiAoaXMubnVsbChjdXRvZmZMb3cpKSANCiAgICAgICAgY3V0b2ZmTG93IDwtIHJlcCgwLCBub19zcCkNCiAgICBjdXRvZmZMb3dbaXMubmEoY3V0b2ZmTG93KV0gPC0gMA0KICAgIA0KICAgICAgICBjdXRvZmZIaWdoIDwtIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRiaW9tYXNzX2N1dG9mZkhpZ2gNCiAgICBpZiAoaXMubnVsbChjdXRvZmZIaWdoKSkgDQogICAgICAgIGN1dG9mZkhpZ2ggPC0gcmVwKDAsIG5vX3NwKQ0KICAgIGN1dG9mZkhpZ2hbaXMubmEoY3V0b2ZmSGlnaCldIDwtIDANCiAgICAgICAgDQogICAgZm9yIChqIGluIDE6bGVuZ3RoKHNpbV9iaW9tYXNzKSkgew0KICAgICAgICBzaW1fYmlvbWFzc1tqXSA9IHN1bSgoc2ltQG5bZGltKHNpbUBuKVsxXSxqLF0gKiBwYXJhbXNAdyAqIA0KICAgICAgICAgICAgcGFyYW1zQGR3KVtwYXJhbXNAdyA+PSBjdXRvZmZMb3dbal0gJiBjdXRvZmZIaWdoW2pdID49IHBhcmFtc0B3XSkNCiAgICB9DQogICAgDQogICAgIHByZWQgPC0gbG9nKHNpbV9iaW9tYXNzKQ0KICAgIGRhdCA8LSBsb2coZGF0KQ0KICAgIGRpc2NyZXAgPC0gcHJlZCAtIGRhdA0KICAgIGRpc2NyZXAgPC0gKHN1bShkaXNjcmVwXjIpKQ0KICAgIHJldHVybihkaXNjcmVwKQ0KfQ0KICAgIA0KYGBgDQoNCg0KYGBge3J9DQojIGNyZWF0ZSBzZXQgb2YgcGFyYW1zIGZvciB0aGUgb3B0aW1pc2F0aW9uIHByb2Nlc3MNCnRpYygpDQoNCnBhcmFtc19vcHRpbSA8LSBwYXJhbXNfdjAyDQoNCnZhcnkgPC0gYyhsb2cxMChwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkUl9tYXgpLA0KICAgICAgICAgIHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRlcmVwcm8sDQogICAgICAgICAgcGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJGludGVyYWN0aW9uX3Jlc291cmNlKQ0KDQpwYXJhbXNfb3B0aW08LXNldFBhcmFtcyhwYXJhbXNfb3B0aW0pDQoNCiMgc2V0IHVwIHdvcmtlcnMNCm5vQ29yZXMgPC0gcGFyYWxsZWw6OmRldGVjdENvcmVzKCkgLSAxICMga2VlcCBzb21lIHNwYXJlIGNvcmUNCmNsIDwtIHBhcmFsbGVsOjptYWtlQ2x1c3Rlcihub0NvcmVzLCBzZXR1cF90aW1lb3V0ID0gMC41KQ0Kc2V0RGVmYXVsdENsdXN0ZXIoY2wgPSBjbCkNCmNsdXN0ZXJFeHBvcnQoY2wsIHZhcmxpc3QgPSAiY2wiLGVudmlyPWVudmlyb25tZW50KCkpDQpjbHVzdGVyRXZhbFEoY2wsIHsNCiAgbGlicmFyeShtaXplckV4cGVyaW1lbnRhbCkNCiAgbGlicmFyeShvcHRpbVBhcmFsbGVsKQ0KfSkNCg0Kb3B0aW1fcmVzdWx0IDwtIG9wdGltUGFyYWxsZWw6Om9wdGltUGFyYWxsZWwocGFyPXZhcnksZ2V0RXJyb3JDdXN0b20scGFyYW1zPXBhcmFtc19vcHRpbSwgDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkYXQgPSBwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkYmlvbWFzc19vYnNlcnZlZCwgbWV0aG9kID0iTC1CRkdTLUIiLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxvd2VyPWMocmVwKC0xNSwxNSkscmVwKDFlLTcsMTUpLHJlcCguMSwxNSkpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdXBwZXI9IGMocmVwKDE1LDE1KSxyZXAoMSwxNSkscmVwKC45OSwxNSkpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFyYWxsZWw9bGlzdChsb2dpbmZvPVRSVUUsIGZvcndhcmQ9VFJVRSkpDQpzdG9wQ2x1c3RlcihjbCkNCg0KdG9jKCkNCg0Kc2F2ZVJEUyhvcHRpbV9yZXN1bHQsIGZpbGU9Im9wdGltX3Jlc3VsdF92MDAuUkRTIikNCiMgb3B0aW1fcmVzdWx0IDwtIHJlYWRSRFMoIm9wdGltX3Jlc3VsdF92MDAuUkRTIikNCmBgYA0KDQoNCmBgYHtyfQ0KI3B1dCB0aGVzZSBuZXcgdmFscyBpbnRvc3BlY2llc19wYXJhbXMgYW5kIGdvIGJhY2sgdG8gdGhlIHRvcCBvZiB0aGlzIHBhZ2UgdG8gcmUtY2hlY2sgdGhlIGNhbGlicmF0aW9uIA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heDwtMTBeb3B0aW1fcmVzdWx0JHBhclsxOjE1XQ0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm88LW9wdGltX3Jlc3VsdCRwYXJbMTY6MzBdDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGludGVyYWN0aW9uX3Jlc291cmNlIDwtb3B0aW1fcmVzdWx0JHBhclszMTo0NV0NCiANCnNpbV9vcHRpbSA8LSBwcm9qZWN0KHBhcmFtc19vcHRpbSwgdF9tYXggPSAyMDAwKQ0KcGxvdEJpb21hc3Moc2ltX29wdGltKQ0KYGBgDQoNCmBgYHtyfQ0KI3B1dCB0aGVzZSBuZXcgdmFscyBpbnRvc3BlY2llc19wYXJhbXMgYW5kIGdvIGJhY2sgdG8gdGhlIHRvcCBvZiB0aGlzIHBhZ2UgdG8gcmUtY2hlY2sgdGhlIGNhbGlicmF0aW9uIA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm8NCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkaW50ZXJhY3Rpb25fcmVzb3VyY2UNCmBgYA0KDQoNCmBgYHtyfQ0KIyB0aGlzIGZ1bmN0aW9uIGFkZHMgYSBsb3dlciBib3VuZGFyeSB0byBzZWxlY3RlZCBzaXplDQpwbG90QmlvbWFzc09ic2VydmVkVnNNb2RlbEN1c3RvbSA8LSBmdW5jdGlvbiAob2JqZWN0LCBzcGVjaWVzID0gTlVMTCwgcmF0aW8gPSBGQUxTRSwgbG9nX3NjYWxlID0gVFJVRSwgDQogICAgcmV0dXJuX2RhdGEgPSBGQUxTRSwgbGFiZWxzID0gVFJVRSwgc2hvd191bm9ic2VydmVkID0gRkFMU0UpIA0Kew0KICAgIGlmIChpcyhvYmplY3QsICJNaXplclNpbSIpKSB7DQogICAgICAgIHBhcmFtcyA9IG9iamVjdEBwYXJhbXMNCiAgICAgICAgbiA8LSBmaW5hbE4ob2JqZWN0KQ0KICAgIH0NCiAgICBlbHNlIGlmIChpcyhvYmplY3QsICJNaXplclBhcmFtcyIpKSB7DQogICAgICAgIHBhcmFtcyA9IG9iamVjdA0KICAgICAgICBuIDwtIGluaXRpYWxOKHBhcmFtcykNCiAgICB9DQogICAgZWxzZSB7DQogICAgICAgIHN0b3AoIllvdSBoYXZlIG5vdCBwcm92aWRlZCBhIHZhbGlkIG1pemVyU2ltIG9yIG1pemVyUGFyYW1zIG9iamVjdC4iKQ0KICAgIH0NCiAgICBzcF9wYXJhbXMgPC0gcGFyYW1zQHNwZWNpZXNfcGFyYW1zDQogICAgc3BlY2llcyA9IHZhbGlkX3NwZWNpZXNfYXJnKG9iamVjdCwgc3BlY2llcykNCiAgICBpZiAobGVuZ3RoKHNwZWNpZXMpID09IDApIA0KICAgICAgICBzdG9wKCJObyBzcGVjaWVzIHNlbGVjdGVkLCBwbGVhc2UgZml4LiIpDQogICAgcm93X3NlbGVjdCA9IG1hdGNoKHNwZWNpZXMsIHNwX3BhcmFtcyRzcGVjaWVzKQ0KICAgIGlmICghImJpb21hc3Nfb2JzZXJ2ZWQiICVpbiUgbmFtZXMoc3BfcGFyYW1zKSkgew0KICAgICAgICBzdG9wKCJZb3UgaGF2ZSBub3QgcHJvdmlkZWQgdmFsdWVzIGZvciB0aGUgY29sdW1uICdiaW9tYXNzX29ic2VydmVkJyAiLCANCiAgICAgICAgICAgICJpbiB0aGUgbWl6ZXJQYXJhbXMvbWl6ZXJTaW0gb2JqZWN0LiIpDQogICAgfQ0KICAgIGVsc2UgaWYgKCFpcy5udW1lcmljKHNwX3BhcmFtcyRiaW9tYXNzX29ic2VydmVkKSkgew0KICAgICAgICBzdG9wKCJUaGUgY29sdW1uICdiaW9tYXNzX29ic2VydmVkJyBpbiB0aGUgbWl6ZXJQYXJhbXMvbWl6ZXJTaW0gb2JqZWN0IiwgDQogICAgICAgICAgICAiIGlzIG5vdCBudW1lcmljLCBwbGVhc2UgZml4LiIpDQogICAgfQ0KICAgIGVsc2Ugew0KICAgICAgICBiaW9tYXNzX29ic2VydmVkID0gc3BfcGFyYW1zJGJpb21hc3Nfb2JzZXJ2ZWQNCiAgICB9DQogICAgDQogICAgY3V0b2ZmTG93IDwtIHNwX3BhcmFtcyRiaW9tYXNzX2N1dG9mZkxvd1tyb3dfc2VsZWN0XQ0KICAgIGlmIChpcy5udWxsKGN1dG9mZkxvdykpIHsNCiAgICAgICAgY3V0b2ZmTG93ID0gcmVwKDAsIGxlbmd0aChzcGVjaWVzKSkNCiAgICB9DQogICAgZWxzZSBpZiAoIWlzLm51bWVyaWMoY3V0b2ZmTG93KSkgew0KICAgICAgICBzdG9wKCJwYXJhbXNAc3BlY2llc19wYXJhbXMkYmlvbWFzc19jdXRvZmZMb3cgaXMgbm90IG51bWVyaWMsIFwiLFxuICAgICAgICAgICAgICAgICBcInBsZWFzZSBmaXguIikNCiAgICB9DQogICAgY3V0b2ZmTG93W2lzLm5hKGN1dG9mZkxvdyldIDwtIDANCiAgICANCiAgICBjdXRvZmZIaWdoIDwtIHNwX3BhcmFtcyRiaW9tYXNzX2N1dG9mZkhpZ2hbcm93X3NlbGVjdF0NCiAgICBpZiAoaXMubnVsbChjdXRvZmZIaWdoKSkgew0KICAgICAgICBjdXRvZmZIaWdoID0gcmVwKDAsIGxlbmd0aChzcGVjaWVzKSkNCiAgICB9DQogICAgZWxzZSBpZiAoIWlzLm51bWVyaWMoY3V0b2ZmSGlnaCkpIHsNCiAgICAgICAgc3RvcCgicGFyYW1zQHNwZWNpZXNfcGFyYW1zJGJpb21hc3NfY3V0b2ZmSGlnaCBpcyBub3QgbnVtZXJpYywgXCIsXG4gICAgICAgICAgICAgICAgIFwicGxlYXNlIGZpeC4iKQ0KICAgIH0NCiAgICBjdXRvZmZIaWdoW2lzLm5hKGN1dG9mZkhpZ2gpXSA8LSAwDQogICAgDQogICAgc2ltX2Jpb21hc3MgPSByZXAoMCwgbGVuZ3RoKHNwZWNpZXMpKQ0KICAgIGZvciAoaiBpbiAxOmxlbmd0aChzcGVjaWVzKSkgew0KICAgICAgICBzaW1fYmlvbWFzc1tqXSA9IHN1bSgobltyb3dfc2VsZWN0W2pdLCBdICogcGFyYW1zQHcgKiANCiAgICAgICAgICAgIHBhcmFtc0BkdylbcGFyYW1zQHcgPj0gY3V0b2ZmTG93W2pdICYgY3V0b2ZmSGlnaFtqXSA+PSBwYXJhbXNAd10pDQogICAgfQ0KICAgIGR1bW15ID0gZGF0YS5mcmFtZShzcGVjaWVzID0gc3BlY2llcywgbW9kZWwgPSBzaW1fYmlvbWFzcywgDQogICAgICAgIG9ic2VydmVkID0gYmlvbWFzc19vYnNlcnZlZFtyb3dfc2VsZWN0XSkgJT4lIG11dGF0ZShzcGVjaWVzID0gZmFjdG9yKHNwZWNpZXMsIA0KICAgICAgICBsZXZlbHMgPSBzcGVjaWVzKSwgaXNfb2JzZXJ2ZWQgPSAhaXMubmEob2JzZXJ2ZWQpICYgb2JzZXJ2ZWQgPiANCiAgICAgICAgMCwgb2JzZXJ2ZWQgPSBjYXNlX3doZW4oaXNfb2JzZXJ2ZWQgfiBvYnNlcnZlZCwgIWlzX29ic2VydmVkIH4gDQogICAgICAgIG1vZGVsKSwgcmF0aW8gPSBtb2RlbC9vYnNlcnZlZCkNCiAgICBpZiAoc3VtKGR1bW15JGlzX29ic2VydmVkKSA9PSAwKSB7DQogICAgICAgIGNhdChwYXN0ZSgiVGhlcmUgYXJlIG5vIG9ic2VydmVkIGJpb21hc3NlcyB0byBjb21wYXJlIHRvIG1vZGVsLCIsIA0KICAgICAgICAgICAgIm9ubHkgcGxvdHRpbmcgbW9kZWwgYmlvbWFzc2VzLiIsIHNlcCA9ICJcbiIpKQ0KICAgIH0NCiAgICBpZiAoIXNob3dfdW5vYnNlcnZlZCkgew0KICAgICAgICBkdW1teSA8LSBmaWx0ZXIoZHVtbXksIGlzX29ic2VydmVkKQ0KICAgIH0NCiAgICBpZiAocmV0dXJuX2RhdGEgPT0gVFJVRSkgDQogICAgICAgIHJldHVybihkdW1teSkNCiAgICB0cmUgPC0gcm91bmQoc3VtKGFicygxIC0gZHVtbXkkcmF0aW8pKSwgZGlnaXRzID0gMykNCiAgICBjYXB0aW9uIDwtIHBhc3RlMCgiVG90YWwgcmVsYXRpdmUgZXJyb3IgPSAiLCB0cmUpDQogICAgaWYgKGFueSghZHVtbXkkaXNfb2JzZXJ2ZWQpKSB7DQogICAgICAgIGNhcHRpb24gPC0gcGFzdGUoY2FwdGlvbiwgIlxuIE9wZW4gY2lyY2xlcyByZXByZXNlbnQgc3BlY2llcyB3aXRob3V0IGJpb21hc3Mgb2JzZXJ2YXRpb24uIikNCiAgICB9DQogICAgaWYgKHJhdGlvID09IEZBTFNFKSB7DQogICAgICAgIGdnIDwtIGdncGxvdChkYXRhID0gZHVtbXksIGFlcyh4ID0gb2JzZXJ2ZWQsIHkgPSBtb2RlbCwgDQogICAgICAgICAgICBjb2xvdXIgPSBzcGVjaWVzLCBzaGFwZSA9IGlzX29ic2VydmVkKSkgKyBnZW9tX2FibGluZShhZXMoaW50ZXJjZXB0ID0gMCwgDQogICAgICAgICAgICBzbG9wZSA9IDEpLCBjb2xvdXIgPSAicHVycGxlIiwgbGluZXR5cGUgPSAiZGFzaGVkIiwgDQogICAgICAgICAgICBzaXplID0gMS4zKSArIGdlb21fcG9pbnQoc2l6ZSA9IDMpICsgbGFicyh5ID0gIm1vZGVsIGJpb21hc3MgW2ddIikgKyANCiAgICAgICAgICAgIGNvb3JkX2NhcnRlc2lhbih5bGltID0gcmFuZ2UoZHVtbXkkbW9kZWwsIGR1bW15JG9ic2VydmVkKSkNCiAgICB9DQogICAgZWxzZSB7DQogICAgICAgIGdnIDwtIGdncGxvdChkYXRhID0gZHVtbXksIGFlcyh4ID0gb2JzZXJ2ZWQsIHkgPSByYXRpbywgDQogICAgICAgICAgICBjb2xvdXIgPSBzcGVjaWVzLCBzaGFwZSA9IGlzX29ic2VydmVkKSkgKyBnZW9tX2hsaW5lKGFlcyh5aW50ZXJjZXB0ID0gMSksIA0KICAgICAgICAgICAgbGluZXR5cGUgPSAiZGFzaGVkIiwgY29sb3VyID0gInB1cnBsZSIsIA0KICAgICAgICAgICAgc2l6ZSA9IDEuMykgKyBnZW9tX3BvaW50KHNpemUgPSAzKSArIGxhYnMoeSA9ICJtb2RlbCBiaW9tYXNzIC8gb2JzZXJ2ZWQgYmlvbWFzcyIpICsgDQogICAgICAgICAgICBjb29yZF9jYXJ0ZXNpYW4oeWxpbSA9IHJhbmdlKGR1bW15JHJhdGlvKSkNCiAgICB9DQogICAgZ2cgPC0gZ2cgKyBsYWJzKHggPSAib2JzZXJ2ZWQgYmlvbWFzcyBbZ10iLCBjYXB0aW9uID0gY2FwdGlvbikgKyANCiAgICAgICAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXMgPSBnZXRDb2xvdXJzKHBhcmFtcylbZHVtbXkkc3BlY2llc10pICsgDQogICAgICAgIHNjYWxlX3NoYXBlX21hbnVhbCh2YWx1ZXMgPSBjKGBUUlVFYCA9IDE5LCBgRkFMU0VgID0gMSkpICsgDQogICAgICAgIGd1aWRlcyhzaGFwZSA9ICJub25lIikNCiAgICBpZiAobG9nX3NjYWxlID09IFRSVUUgJiByYXRpbyA9PSBGQUxTRSkgew0KICAgICAgICBnZyA9IGdnICsgc2NhbGVfeF9sb2cxMCgpICsgc2NhbGVfeV9sb2cxMCgpDQogICAgfQ0KICAgIGlmIChsb2dfc2NhbGUgPT0gVFJVRSAmIHJhdGlvID09IFRSVUUpIHsNCiAgICAgICAgZ2cgPSBnZyArIHNjYWxlX3hfbG9nMTAoKQ0KICAgIH0NCiAgICBpZiAobGFiZWxzID09IFRSVUUpIHsNCiAgICAgICAgZ2cgPSBnZyArIGdncmVwZWw6Omdlb21fbGFiZWxfcmVwZWwoYWVzKGxhYmVsID0gc3BlY2llcyksIA0KICAgICAgICAgICAgYm94LnBhZGRpbmcgPSAwLjM1LCBwb2ludC5wYWRkaW5nID0gMC41LCBzZWdtZW50LmNvbG9yID0gImdyZXk1MCIsIA0KICAgICAgICAgICAgc2hvdy5sZWdlbmQgPSBGQUxTRSwgbWF4Lm92ZXJsYXBzID0gSW5mLCBzZWVkID0gNDIpDQogICAgfQ0KICAgIGdnDQp9DQpgYGANCg0KDQpgYGB7cn0NCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsQ3VzdG9tKHNpbV9vcHRpbSkNCmBgYA0KDQoNClVzZSB0dW5lUGFyYW1zKCkgdG8gaW52ZXN0aWdhdGUgdGhlIHNwZWNpZXMgd2l0aCBlcmVwcm8gdmFsdWVzIHRvbyBoaWdoDQpgYGB7cn0NCnBhcmFtc190dW5lZF92MDEgPC0gdHVuZVBhcmFtcyhwYXJhbXNfb3B0aW0pDQoNCnBhcmFtc190dW5lZF92MDFAc3BlY2llc19wYXJhbXMkZXJlcHJvDQoNCnBhcmFtc190dW5lZF92MDEgPC0gc3RlYWR5KHBhcmFtc190dW5lZF92MDEpDQpgYGANCg0KYGBge3J9DQojIGNyZWF0ZSBzZXQgb2YgcGFyYW1zIGZvciB0aGUgb3B0aW1pc2F0aW9uIHByb2Nlc3MNCnRpYygpDQoNCnBhcmFtc19vcHRpbSA8LSBwYXJhbXNfdHVuZWRfdjAxDQoNCnZhcnkgPC0gYyhsb2cxMChwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkUl9tYXgpLA0KICAgICAgICAgIHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRlcmVwcm8sDQogICAgICAgICAgcGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJGludGVyYWN0aW9uX3Jlc291cmNlKQ0KDQpwYXJhbXNfb3B0aW08LXNldFBhcmFtcyhwYXJhbXNfb3B0aW0pDQoNCiMgc2V0IHVwIHdvcmtlcnMNCm5vQ29yZXMgPC0gcGFyYWxsZWw6OmRldGVjdENvcmVzKCkgLSAxICMga2VlcCBzb21lIHNwYXJlIGNvcmUNCmNsIDwtIHBhcmFsbGVsOjptYWtlQ2x1c3Rlcihub0NvcmVzLCBzZXR1cF90aW1lb3V0ID0gMC41KQ0Kc2V0RGVmYXVsdENsdXN0ZXIoY2wgPSBjbCkNCmNsdXN0ZXJFeHBvcnQoY2wsIHZhcmxpc3QgPSAiY2wiLGVudmlyPWVudmlyb25tZW50KCkpDQpjbHVzdGVyRXZhbFEoY2wsIHsNCiAgbGlicmFyeShtaXplckV4cGVyaW1lbnRhbCkNCiAgbGlicmFyeShvcHRpbVBhcmFsbGVsKQ0KfSkNCg0Kb3B0aW1fcmVzdWx0IDwtIG9wdGltUGFyYWxsZWw6Om9wdGltUGFyYWxsZWwocGFyPXZhcnksZ2V0RXJyb3JDdXN0b20scGFyYW1zPXBhcmFtc19vcHRpbSwgDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkYXQgPSBwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkYmlvbWFzc19vYnNlcnZlZCwgbWV0aG9kID0iTC1CRkdTLUIiLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxvd2VyPWMocmVwKC0xNSwxNSkscmVwKDFlLTcsMTUpLHJlcCguMSwxNSkpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdXBwZXI9IGMocmVwKDE1LDE1KSxyZXAoMC45OSwxNSkscmVwKC45OSwxNSkpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFyYWxsZWw9bGlzdChsb2dpbmZvPVRSVUUsIGZvcndhcmQ9VFJVRSkpDQpzdG9wQ2x1c3RlcihjbCkNCg0KdG9jKCkNCg0Kc2F2ZVJEUyhvcHRpbV9yZXN1bHQsIGZpbGU9Im9wdGltX3Jlc3VsdF92MDEuUkRTIikNCiMgb3B0aW1fcmVzdWx0IDwtIHJlYWRSRFMoIm9wdGltX3Jlc3VsdF92MDAuUkRTIikNCmBgYA0KDQoNCmBgYHtyfQ0KI3B1dCB0aGVzZSBuZXcgdmFscyBpbnRvc3BlY2llc19wYXJhbXMgYW5kIGdvIGJhY2sgdG8gdGhlIHRvcCBvZiB0aGlzIHBhZ2UgdG8gcmUtY2hlY2sgdGhlIGNhbGlicmF0aW9uIA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heDwtMTBeb3B0aW1fcmVzdWx0JHBhclsxOjE1XQ0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm88LW9wdGltX3Jlc3VsdCRwYXJbMTY6MzBdDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGludGVyYWN0aW9uX3Jlc291cmNlIDwtb3B0aW1fcmVzdWx0JHBhclszMTo0NV0NCg0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm8NCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkaW50ZXJhY3Rpb25fcmVzb3VyY2UNCmBgYA0KDQoNCmBgYHtyfQ0Kc2ltX29wdGltX3YwMiA8LSBwcm9qZWN0KHBhcmFtc19vcHRpbSwgdF9tYXggPSAxMDAwKQ0KcGxvdEJpb21hc3Moc2ltX29wdGltX3YwMikNCmBgYA0KDQpgYGB7cn0NCnBhcmFtc192MDMgPC0gc3RlYWR5KHBhcmFtc19vcHRpbSkNCmBgYA0KDQpgYGB7cn0NCnBhcmFtc19sb29wIDwtIHBhcmFtc192MDMgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkgfD4NCiAgICBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+IG1hdGNoQmlvbWFzc2VzKCkgfD4gc3RlYWR5KCkNCg0KcGFyYW1zX2xvb3BAc3BlY2llc19wYXJhbXMkZXJlcHJvDQpgYGANCg0KYGBge3J9DQpzaW1fdjAyIDwtIHByb2plY3QocGFyYW1zX2xvb3AsIHRfbWF4ID0gMjAwMCkNCmBgYA0KDQoNCg0KYGBge3J9DQpwbG90QmlvbWFzcyhzaW1fdjAyKQ0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWxDdXN0b20oc2ltX3YwMikNCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsKHNpbV92MDIpDQpwbG90RGlldChwYXJhbXNfbG9vcCkNCmBgYA0KU2F2ZSBwcm9ncmVzcw0KYGBge3J9DQojIHNhdmVSRFMocGFyYW1zX2xvb3AsICJwYXJhbXNfb3B0aW1fdjAxLnJkcyIpDQpgYGANCg0KDQpUcnkgdG8gd2VlbiB0aGUgc3BlY2llcyBvZiB0aGUgcmVzb3VyY2Ugc3BlY3RyYQ0KYGBge3J9DQpwYXJhbXNfbG9vcEByZXNvdXJjZV9wYXJhbXMkd19wcF9jdXRvZmYNCnBhcmFtc19sb29wQHJlc291cmNlX3BhcmFtcyRrYXBwYQ0KYGBgDQoNCmBgYHtyfQ0KcmVzb3VyY2VfcGFyYW1zKHBhcmFtc19sb29wKVtbIndfcHBfY3V0b2ZmIl1dIDwtIDEwDQoNCnBhcmFtc192MDQgPC0gc2V0UGFyYW1zKHBhcmFtc19sb29wKQ0KDQpwYXJhbXNfdjA0QHJlc291cmNlX3BhcmFtcyR3X3BwX2N1dG9mZg0KYGBgDQoNCmBgYHtyfQ0Kc2ltX3YwMyA8LSBwcm9qZWN0KHBhcmFtc192MDQsIHRfbWF4ID0gMjAwKQ0KDQpwbG90QmlvbWFzcyhzaW1fdjAzKQ0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWxDdXN0b20oc2ltX3YwMykNCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsKHNpbV92MDMpDQpwbG90RGlldChwYXJhbXNfdjA0KQ0KYGBgDQoNCmBgYHtyfQ0KZ2V0RXJyb3JDdXN0b21fdjAyIDwtIGZ1bmN0aW9uKHZhcnksIHBhcmFtcywgZGF0LCB0b2wgPSAwLjAwMSwgDQogICAgdGltZXRvcnVuID0gMTApDQp7DQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRSX21heFsxOjE1XTwtMTBedmFyeVsxOjE1XSAjIFJfbWF4IGZvciAxNSBzcGVjaWVzDQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRlcmVwcm9bMToxNV08LXZhcnlbMTY6MzBdICMgZXJlcHJvIGZvciAxNSBzcGVjaWVzDQogIHBhcmFtc0BzcGVjaWVzX3BhcmFtcyRpbnRlcmFjdGlvbl9yZXNvdXJjZVsxOjE1XSA8LSB2YXJ5WzMxOjQ1XSAjIGludGVyYWN0aW9uX3Jlc291cmNlIGZvciAxNSBzcGVjaWVzDQogIHBhcmFtc0ByZXNvdXJjZV9wYXJhbXMkd19wcF9jdXRvZmZbMV0gPC0gdmFyeVs0Nl0gIyBpbnRlcmFjdGlvbl9yZXNvdXJjZSBmb3IgMTUgc3BlY2llcw0KICBwYXJhbXMgPC0gc2V0UGFyYW1zKHBhcmFtcykNCiAgIyBpbnRlcmFjdGlvbiA8LSBwYXJhbXNAaW50ZXJhY3Rpb24NCiAgIyBpbnRlcmFjdGlvbltdIDwtIG1hdHJpeCh2YXJ5WzI4OjEwOF0sbnJvdyA9IDkpICMgc3RvcCBhdCA1NCBpZiBsb29raW5nIG9ubHkgYXQgMyBiaWdnZXN0IHNwZWNpZXMNCiAgDQogICMgcGFyYW1zIDwtIHNldEludGVyYWN0aW9uKHBhcmFtcyxpbnRlcmFjdGlvbikNCiAgICBwYXJhbXMgPC0gcHJvamVjdFRvU3RlYWR5KHBhcmFtcywgZGlzdGFuY2VfZnVuYyA9IGRpc3RhbmNlU1NMb2dOLCANCiAgICAgICAgdG9sID0gdG9sLCB0X21heCA9IDIwMCwgcmV0dXJuX3NpbSA9IEYpDQogICAgDQogICAgc2ltIDwtIHByb2plY3QocGFyYW1zLCB0X21heCA9IHRpbWV0b3J1biwgcHJvZ3Jlc3NfYmFyID0gRikNCiAgICANCiAgICBzaW1fYmlvbWFzcyA9IHJlcCgwLCBsZW5ndGgocGFyYW1zQHNwZWNpZXNfcGFyYW1zJHNwZWNpZXMpKQ0KICAgIA0KICAgICAgICBjdXRvZmZMb3cgPC0gcGFyYW1zQHNwZWNpZXNfcGFyYW1zJGJpb21hc3NfY3V0b2ZmTG93DQogICAgaWYgKGlzLm51bGwoY3V0b2ZmTG93KSkgDQogICAgICAgIGN1dG9mZkxvdyA8LSByZXAoMCwgbm9fc3ApDQogICAgY3V0b2ZmTG93W2lzLm5hKGN1dG9mZkxvdyldIDwtIDANCiAgICANCiAgICAgICAgY3V0b2ZmSGlnaCA8LSBwYXJhbXNAc3BlY2llc19wYXJhbXMkYmlvbWFzc19jdXRvZmZIaWdoDQogICAgaWYgKGlzLm51bGwoY3V0b2ZmSGlnaCkpIA0KICAgICAgICBjdXRvZmZIaWdoIDwtIHJlcCgwLCBub19zcCkNCiAgICBjdXRvZmZIaWdoW2lzLm5hKGN1dG9mZkhpZ2gpXSA8LSAwDQogICAgICAgIA0KICAgIGZvciAoaiBpbiAxOmxlbmd0aChzaW1fYmlvbWFzcykpIHsNCiAgICAgICAgc2ltX2Jpb21hc3Nbal0gPSBzdW0oKHNpbUBuW2RpbShzaW1AbilbMV0saixdICogcGFyYW1zQHcgKiANCiAgICAgICAgICAgIHBhcmFtc0BkdylbcGFyYW1zQHcgPj0gY3V0b2ZmTG93W2pdICYgY3V0b2ZmSGlnaFtqXSA+PSBwYXJhbXNAd10pDQogICAgfQ0KICAgIA0KICAgICBwcmVkIDwtIGxvZyhzaW1fYmlvbWFzcykNCiAgICBkYXQgPC0gbG9nKGRhdCkNCiAgICBkaXNjcmVwIDwtIHByZWQgLSBkYXQNCiAgICBkaXNjcmVwIDwtIChzdW0oZGlzY3JlcF4yKSkNCiAgICByZXR1cm4oZGlzY3JlcCkNCn0NCiAgICANCmBgYA0KDQoNCmBgYHtyfQ0KIyBjcmVhdGUgc2V0IG9mIHBhcmFtcyBmb3IgdGhlIG9wdGltaXNhdGlvbiBwcm9jZXNzDQp0aWMoKQ0KDQpwYXJhbXNfb3B0aW0gPC0gcGFyYW1zX3YwNA0KDQp2YXJ5IDwtIGMobG9nMTAocGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJFJfbWF4KSwNCiAgICAgICAgICBwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkZXJlcHJvLA0KICAgICAgICAgIHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRpbnRlcmFjdGlvbl9yZXNvdXJjZSwNCiAgICAgICAgICBwYXJhbXNfb3B0aW1AcmVzb3VyY2VfcGFyYW1zJHdfcHBfY3V0b2ZmKQ0KDQpwYXJhbXNfb3B0aW08LXNldFBhcmFtcyhwYXJhbXNfb3B0aW0pDQoNCiMgc2V0IHVwIHdvcmtlcnMNCm5vQ29yZXMgPC0gcGFyYWxsZWw6OmRldGVjdENvcmVzKCkgLSAxICMga2VlcCBzb21lIHNwYXJlIGNvcmUNCmNsIDwtIHBhcmFsbGVsOjptYWtlQ2x1c3Rlcihub0NvcmVzLCBzZXR1cF90aW1lb3V0ID0gMC41KQ0Kc2V0RGVmYXVsdENsdXN0ZXIoY2wgPSBjbCkNCmNsdXN0ZXJFeHBvcnQoY2wsIHZhcmxpc3QgPSAiY2wiLGVudmlyPWVudmlyb25tZW50KCkpDQpjbHVzdGVyRXZhbFEoY2wsIHsNCiAgbGlicmFyeShtaXplckV4cGVyaW1lbnRhbCkNCiAgbGlicmFyeShvcHRpbVBhcmFsbGVsKQ0KfSkNCg0Kb3B0aW1fcmVzdWx0IDwtIG9wdGltUGFyYWxsZWw6Om9wdGltUGFyYWxsZWwocGFyPXZhcnksZ2V0RXJyb3JDdXN0b21fdjAyLHBhcmFtcz1wYXJhbXNfb3B0aW0sIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGF0ID0gcGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJGJpb21hc3Nfb2JzZXJ2ZWQsIG1ldGhvZCA9IkwtQkZHUy1CIiwgDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBsb3dlcj1jKHJlcCgtMTUsMTUpLHJlcCgxZS03LDE1KSxyZXAoLjEsMTUpLCAwLjEpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgdXBwZXI9IGMocmVwKDE1LDE1KSxyZXAoMC45OSwxNSkscmVwKC45OSwxNSksMTAwMCksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBwYXJhbGxlbD1saXN0KGxvZ2luZm89VFJVRSwgZm9yd2FyZD1UUlVFKSkNCnN0b3BDbHVzdGVyKGNsKQ0KDQp0b2MoKQ0KDQpzYXZlUkRTKG9wdGltX3Jlc3VsdCwgZmlsZT0ib3B0aW1fcmVzdWx0X3YwMi5SRFMiKQ0KIyBvcHRpbV9yZXN1bHQgPC0gcmVhZFJEUygib3B0aW1fcmVzdWx0X3YwMC5SRFMiKQ0KYGBgDQoNCg0KDQpgYGB7cn0NCiNwdXQgdGhlc2UgbmV3IHZhbHMgaW50b3NwZWNpZXNfcGFyYW1zIGFuZCBnbyBiYWNrIHRvIHRoZSB0b3Agb2YgdGhpcyBwYWdlIHRvIHJlLWNoZWNrIHRoZSBjYWxpYnJhdGlvbiANCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkUl9tYXg8LTEwXm9wdGltX3Jlc3VsdCRwYXJbMToxNV0NCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkZXJlcHJvPC1vcHRpbV9yZXN1bHQkcGFyWzE2OjMwXQ0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRpbnRlcmFjdGlvbl9yZXNvdXJjZSA8LW9wdGltX3Jlc3VsdCRwYXJbMzE6NDVdDQpyZXNvdXJjZV9wYXJhbXMocGFyYW1zX29wdGltKSR3X3BwX2N1dG9mZiA8LSBvcHRpbV9yZXN1bHQkcGFyWzQ2XQ0KDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJFJfbWF4DQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGVyZXBybw0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRpbnRlcmFjdGlvbl9yZXNvdXJjZQ0KYGBgDQoNCg0KYGBge3J9DQpzaW1fdjA0IDwtIHByb2plY3QocGFyYW1zX29wdGltLCB0X21heCA9IDEwMDApDQoNCnBsb3RCaW9tYXNzKHNpbV92MDQpDQpwbG90QmlvbWFzc09ic2VydmVkVnNNb2RlbEN1c3RvbShzaW1fdjA0KQ0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWwoc2ltX3YwNCkNCnBsb3REaWV0KHBhcmFtc19vcHRpbSkNCmBgYA0KDQpgYGB7cn0NCnBhcmFtc19zdGVhZHkgPC0gcGFyYW1zX29wdGltDQoNCmluaXRpYWxOKHBhcmFtc19zdGVhZHkpIDwtIHNpbV92MDRAbltkaW0oc2ltX3YwNEBuKVsxXSwsXQ0KDQpwbG90RGlldChwYXJhbXNfc3RlYWR5KQ0KDQpzaW1fdjA1IDwtIHByb2plY3QocGFyYW1zX3N0ZWFkeSwgdF9tYXggPSAxMDApDQoNCnBsb3RCaW9tYXNzKHNpbV92MDUpDQoNCnBsb3RHcm93dGhDdXJ2ZXMoc2ltX3YwNSwgc3BlY2llc19wYW5lbCA9IFQpDQpgYGANCg0KYGBge3J9DQpwYXJhbXNfdHVuZWRfdjAyIDwtIHR1bmVQYXJhbXMocGFyYW1zX3N0ZWFkeSkNCg0KcGFyYW1zX3R1bmVkX3YwMiA8LSBzdGVhZHkocGFyYW1zX3R1bmVkX3YwMikNCg0KcGFyYW1zX3R1bmVkX3YwMkBzcGVjaWVzX3BhcmFtcyRlcmVwcm8NCmBgYA0Kc3BlY2llcyAgICAgICAgd19taW4gICAgIHdfbWF0ICAgIHdfbWF4IGJldGEga192YiAgaCBtaW5fZGVwdGggbWF4X2RlcHRoIHdhdGVyLmNvbHVtbi51c2UgcF90aW1lX3ByeWR6IHBjX2FubnVhbF9vZmZzcHJpbmcNCjUgICBzYWxwcyAzLjE2MjI3OGUtMDUgMC4yNTExODg2IDI1LjExODg2IDEwMDAgICBOQSAzMyAgICAgICAgIDAgICAgICAgNTAwICAgICAgICAgIG5vbiBEVk0gICAgICAgICAgICAxICAgICAgICAgICAgICAgICAgTkENCiAgYmlvbWFzc19vYnNlcnZlZA0KNSAgICAgICAgICAgIDAuNjUyDQpgYGB7cn0NCnNwZWNpZXNfcGFyYW1zIDwtIGRhdGEuZnJhbWUoDQogICAgc3BlY2llcyA9ICJzYWxwcyIsDQogICAgd19taW4gPSAzLjE2MjI3OGUtMDUsDQogICAgd19tYXggPSAyNS4xMTg4NiwNCiAgICB3X21hdCA9IDAuMjUxMTg4NiwNCiAgICBiZXRhID0gMTAwMDAwMCwNCiAgICBzaWdtYSA9IDIsDQogICAgYmlvbWFzc19vYnNlcnZlZCA9IDAuNjUyLA0KICAgIGtfdmIgPSAwLjQsDQogICAgeWllbGRfb2JzZXJ2ZWQgPSAwLA0KICAgIGJpb21hc3NfY3V0b2ZmID0gMC4xLA0KICAgIGJpb21hc3NfY3V0b2ZmTG93ID0gMC4xLA0KICAgIGJpb21hc3NfY3V0b2ZmSGlnaCA9IDI1LjExODg2DQopDQoNCnBhcmFtc192MDUgPC0gYWRkU3BlY2llcyhwYXJhbXNfdHVuZWRfdjAyLCBzcGVjaWVzX3BhcmFtcykNCnBsb3RTcGVjdHJhKHBhcmFtc192MDUpDQoNCnBhcmFtc192MDVAc3BlY2llc19wYXJhbXMkZXJlcHJvDQoNCmBgYA0KDQoNCmBgYHtyfQ0KIyBmdW5jdGlvbiBydW5uaW5nIHR1bmVQYXJhbXMgZnVuY3Rpb24gaW4gYSByb3cgZm9yIGEgcXVpY2sgc3RhcnQgdG8gYSBjYWxpYnJhdGlvbg0KZmFzdENhbGliIDwtIGZ1bmN0aW9uKHBhcmFtcywgbWF0Y2ggPSBGKQ0Kew0KcGFyYW1zIDwtIGNhbGlicmF0ZUJpb21hc3MocGFyYW1zKSAjIGNoYW5nZXMga2FwcGEgYW5kIHJtYXgNCmlmKG1hdGNoKSBwYXJhbXMgPC0gbWF0Y2hCaW9tYXNzZXMocGFyYW1zKSAjIHNldCBybWF4IHRvIGluZiBhbmQgYWRqdXN0IGVyZXBybw0KcGFyYW1zIDwtIHN0ZWFkeShwYXJhbXMsIHRvbCA9IDAuMDAxKQ0Kc2ltIDwtIHByb2plY3QocGFyYW1zLCB0X21heCA9IDEwMDApDQpyZXR1cm4oc2ltKQ0KfQ0KYGBgDQoNCg0KYGBge3J9DQpzaW1fdjA2IDwtIGZhc3RDYWxpYihwYXJhbXNfdjA1KQ0KDQpwbG90QmlvbWFzcyhzaW1fdjA2KQ0KcGxvdERpZXQocGFyYW1zX3YwNSkNCmBgYA0KDQpgYGB7cn0NCnNvX3RoZXRhIDwtIGFzLm1hdHJpeChwYXJhbXNfdjA1QGludGVyYWN0aW9uKQ0KDQpzYWxwX3JvdyA8LSBjKDEsDQowLjEyNTAwMDAsDQowLjEyNTAwMDAsDQowLjEzNzM2MjYsDQowLjA1MDAwMDAsDQowLjEwMDAwMDAsDQowLjI1MDAwMDAsDQowLjI1MDAwMDAsDQowLjEwMDAwMDAsDQowLjE2NjY2NjcsDQowLjE2NjY2NjcsDQowLjA1MDAwMDAsDQowLjIwMDAwMDAsDQowLjEyNTAwMDAsDQowLjIwMDAwMDAsDQoxKQ0KYGBgDQoNCg0KYGBge3J9DQpzb190aGV0YVsxNixdICMgYXV0b21hdGVkIHNhbHBzIGludGVyYWN0aW9uIHdoZW4gYWRkZWQgdXNpbmcgYWRkU3BlY2llcygpDQoNCnNvX3RoZXRhWzE2LF0gPC0gc2FscF9yb3cNCg0Kc29fdGhldGFbLDE2XSAjIGF1dG9tYXRlZCBzYWxwcyBpbnRlcmFjdGlvbiB3aGVuIGFkZGVkIHVzaW5nIGFkZFNwZWNpZXMoKQ0KDQpzb190aGV0YVssMTZdIDwtIHNhbHBfcm93DQoNCnNvX3RoZXRhDQpgYGANCg0KDQpgYGB7cn0NCnBhcmFtc192MDYgPC0gc2V0UGFyYW1zKHBhcmFtc192MDUsIGludGVyYWN0aW9uID0gc29fdGhldGEpDQoNCnBhcmFtc192MDZAaW50ZXJhY3Rpb24NCmBgYA0KDQoNCmBgYHtyfQ0Kc2ltX3YwNyA8LSBmYXN0Q2FsaWIocGFyYW1zX3YwNikNCg0KcGxvdEJpb21hc3Moc2ltX3YwNykNCmBgYA0KDQpgYGB7cn0NCmluaXRpYWxOKHBhcmFtc192MDYpIDwtIHNpbV92MDdAbltkaW0oc2ltX3YwN0BuKVsxXSwsXQ0KYGBgDQoNCg0KYGBge3J9DQpzaW1fdjA4IDwtIHByb2plY3QocGFyYW1zX3YwNiwgdF9tYXggPSAxMDApDQpwbG90QmlvbWFzcyhzaW1fdjA4KQ0KcGxvdERpZXQocGFyYW1zX3YwNikNCmBgYA0KDQpGaWxsIGluIG1pc3NpbmcgYW5kIGFkanVzdCBiZXRhIHZhbHVlcyBmb3Igem9vcGxhbmt0b24gZ3JvdXBzIHVzaW5nIEhlbmVnaGFuIGV0IGFsLiAyMDIwICgxMC4xMDE2L2ouZWNvbG1vZGVsLjIwMjAuMTA5MjY1KQ0KDQpOZWVkIGEgdmFsdWUgZm9yIG1pY3Jvem9vcGxhbmt0b24NCm1pY3Jvem9vcGxhbmt0b24gKGluIE1jQ29ybWFjayBldCBhbC4gMjAyMCkgaXMgY29tcG9zZWQgb2YgSGV0ZXJvdHJvcGhpYyBkaW5vZmxhZ2VsbGF0ZXMsIHRpbnRpbm5pZHMsIGNpbGlhdGVzLCBjb3BlcG9kIG5hdXBsaWkNCkhlbmVnaGFuIGV0IGFsLiBsb2cxMFBQTVIgdmFsdWVzIGZvcjoNCkhldGVyby5GbGFnZWxsYXRlcyA9IDAuMuKAkzAuNzIgLT4gMC40Ng0KSGV0ZXJvLkNpbGlhdGVzID0gMi414oCTMi45IC0+IDIuNw0KTWVhbjogKDIuNyswLjQ2KS8yID0gMS41OA0KMTBeMS41OCA9IDM4LjAxODk0DQoNCk5lZWQgdG8gYWRqdXN0IHZhbHVlcyBmb3IgbWVzb3pvbywgb3RoZXIgbWFjcm96b28sIGV1cGhhdXNpaWRzLCBzYWxwcw0KDQpIZW5lZ2hhbiBldCBhbC4gbG9nMTBQUE1SIHZhbHVlcyAobWlkcG9pbnRzIGZvciByYW5nZSkgZm9yOg0Kc2FscHMgPSA2LjjigJMxMS43IC0+IDkuMjUNCjEwXjkuMjUgPSAxNzc4Mjc5NDEwDQoNCmV1cGhhdXNpaWRzID0gNi424oCTNy44IC0+IDcuMg0KMTBeNy4yID0gMTU4NDg5MzINCg0KbWVzb3pvbyAoY29wZXBvZHMpDQpPbW5pLkNvcC4gPSAzLjbigJM0LjYgLT4gNC4xDQpDYXJuLkNvcC4gPSAwLjjigJMxLjkgLT4gMS4zNQ0KTWVhbjogKDQuMSsxLjM1KS8yID0gMi43MjUNCjEwXjIuNzI1ID0gNTMwLjg4NDQNCg0Kb3RoZXIgbWFjcm96b28gKCkNCkNoYWV0b2duYXRocyA9IDEuOeKAkzMuNCAtPiAyLjY1DQoxMF4yLjY1ID0gNDQ2LjY4MzYNCg0KYGBge3J9DQpwYXJhbXNfdHVuZWRfdjAzIDwtIHR1bmVQYXJhbXMocGFyYW1zX3YwNikNCnBhcmFtc190dW5lZF92MDQgPC0gdHVuZVBhcmFtcyhwYXJhbXNfdHVuZWRfdjAzKSAjIHVwZGF0ZWQgc2FscHMgUFBNUiB0byBtYXRjaCBlbXBpcmljYWwgZXN0aW1hdGVzIGZyb20gSGVuZWdoYW4gZXQgYWwgMjAyMA0KYGBgDQoNCg0KYGBge3J9DQpnZXRFcnJvckN1c3RvbV92MyA8LSBmdW5jdGlvbih2YXJ5LCBwYXJhbXMsIGRhdCwgdG9sID0gMC4wMDEsIA0KICAgIHRpbWV0b3J1biA9IDEwKQ0Kew0KICBwYXJhbXNAc3BlY2llc19wYXJhbXMkUl9tYXhbMToxNl08LTEwXnZhcnlbMToxNl0gIyBSX21heCBmb3IgMTUgc3BlY2llcw0KICBwYXJhbXNAc3BlY2llc19wYXJhbXMkZXJlcHJvWzE6MTZdPC12YXJ5WzE3OjMyXSAjIGVyZXBybyBmb3IgMTUgc3BlY2llcw0KICBwYXJhbXNAc3BlY2llc19wYXJhbXMkaW50ZXJhY3Rpb25fcmVzb3VyY2VbMToxNl0gPC0gdmFyeVszMzo0OF0gIyBpbnRlcmFjdGlvbl9yZXNvdXJjZSBmb3IgMTUgc3BlY2llcw0KICBwYXJhbXMgPC0gc2V0UGFyYW1zKHBhcmFtcykNCiAgIyBpbnRlcmFjdGlvbiA8LSBwYXJhbXNAaW50ZXJhY3Rpb24NCiAgIyBpbnRlcmFjdGlvbltdIDwtIG1hdHJpeCh2YXJ5WzI4OjEwOF0sbnJvdyA9IDkpICMgc3RvcCBhdCA1NCBpZiBsb29raW5nIG9ubHkgYXQgMyBiaWdnZXN0IHNwZWNpZXMNCiAgDQogICMgcGFyYW1zIDwtIHNldEludGVyYWN0aW9uKHBhcmFtcyxpbnRlcmFjdGlvbikNCiAgICBwYXJhbXMgPC0gcHJvamVjdFRvU3RlYWR5KHBhcmFtcywgZGlzdGFuY2VfZnVuYyA9IGRpc3RhbmNlU1NMb2dOLCANCiAgICAgICAgdG9sID0gdG9sLCB0X21heCA9IDIwMCwgcmV0dXJuX3NpbSA9IEYpDQogICAgDQogICAgc2ltIDwtIHByb2plY3QocGFyYW1zLCB0X21heCA9IHRpbWV0b3J1biwgcHJvZ3Jlc3NfYmFyID0gRikNCiAgICANCiAgICBzaW1fYmlvbWFzcyA9IHJlcCgwLCBsZW5ndGgocGFyYW1zQHNwZWNpZXNfcGFyYW1zJHNwZWNpZXMpKQ0KICAgIA0KICAgICAgICBjdXRvZmZMb3cgPC0gcGFyYW1zQHNwZWNpZXNfcGFyYW1zJGJpb21hc3NfY3V0b2ZmTG93DQogICAgaWYgKGlzLm51bGwoY3V0b2ZmTG93KSkgDQogICAgICAgIGN1dG9mZkxvdyA8LSByZXAoMCwgbm9fc3ApDQogICAgY3V0b2ZmTG93W2lzLm5hKGN1dG9mZkxvdyldIDwtIDANCiAgICANCiAgICAgICAgY3V0b2ZmSGlnaCA8LSBwYXJhbXNAc3BlY2llc19wYXJhbXMkYmlvbWFzc19jdXRvZmZIaWdoDQogICAgaWYgKGlzLm51bGwoY3V0b2ZmSGlnaCkpIA0KICAgICAgICBjdXRvZmZIaWdoIDwtIHJlcCgwLCBub19zcCkNCiAgICBjdXRvZmZIaWdoW2lzLm5hKGN1dG9mZkhpZ2gpXSA8LSAwDQogICAgICAgIA0KICAgIGZvciAoaiBpbiAxOmxlbmd0aChzaW1fYmlvbWFzcykpIHsNCiAgICAgICAgc2ltX2Jpb21hc3Nbal0gPSBzdW0oKHNpbUBuW2RpbShzaW1AbilbMV0saixdICogcGFyYW1zQHcgKiANCiAgICAgICAgICAgIHBhcmFtc0BkdylbcGFyYW1zQHcgPj0gY3V0b2ZmTG93W2pdICYgY3V0b2ZmSGlnaFtqXSA+PSBwYXJhbXNAd10pDQogICAgfQ0KICAgIA0KICAgICBwcmVkIDwtIGxvZyhzaW1fYmlvbWFzcykNCiAgICBkYXQgPC0gbG9nKGRhdCkNCiAgICBkaXNjcmVwIDwtIHByZWQgLSBkYXQNCiAgICBkaXNjcmVwIDwtIChzdW0oZGlzY3JlcF4yKSkNCiAgICByZXR1cm4oZGlzY3JlcCkNCn0NCmBgYA0KDQoNCg0KDQpgYGB7cn0NCiMgY3JlYXRlIHNldCBvZiBwYXJhbXMgZm9yIHRoZSBvcHRpbWlzYXRpb24gcHJvY2Vzcw0KdGljKCkNCg0KcGFyYW1zX29wdGltIDwtIHBhcmFtc190dW5lZF92MDQNCg0KdmFyeSA8LSBjKGxvZzEwKHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRSX21heCksDQogICAgICAgICAgcGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJGVyZXBybywNCiAgICAgICAgICBwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkaW50ZXJhY3Rpb25fcmVzb3VyY2UpDQoNCnBhcmFtc19vcHRpbTwtc2V0UGFyYW1zKHBhcmFtc19vcHRpbSkNCg0KIyBzZXQgdXAgd29ya2Vycw0Kbm9Db3JlcyA8LSBwYXJhbGxlbDo6ZGV0ZWN0Q29yZXMoKSAtIDEgIyBrZWVwIHNvbWUgc3BhcmUgY29yZQ0KY2wgPC0gcGFyYWxsZWw6Om1ha2VDbHVzdGVyKG5vQ29yZXMsIHNldHVwX3RpbWVvdXQgPSAwLjUpDQpzZXREZWZhdWx0Q2x1c3RlcihjbCA9IGNsKQ0KY2x1c3RlckV4cG9ydChjbCwgdmFybGlzdCA9ICJjbCIsZW52aXI9ZW52aXJvbm1lbnQoKSkNCmNsdXN0ZXJFdmFsUShjbCwgew0KICBsaWJyYXJ5KG1pemVyRXhwZXJpbWVudGFsKQ0KICBsaWJyYXJ5KG9wdGltUGFyYWxsZWwpDQp9KQ0KDQpvcHRpbV9yZXN1bHQgPC0gb3B0aW1QYXJhbGxlbDo6b3B0aW1QYXJhbGxlbChwYXI9dmFyeSxnZXRFcnJvckN1c3RvbV92MyxwYXJhbXM9cGFyYW1zX29wdGltLCANCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdCA9IHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRiaW9tYXNzX29ic2VydmVkLCBtZXRob2QgPSJMLUJGR1MtQiIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbG93ZXI9YyhyZXAoLTE1LDE2KSxyZXAoMWUtNywxNikscmVwKC4xLDE2KSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB1cHBlcj0gYyhyZXAoMTUsMTYpLHJlcCgwLjk5LDE2KSxyZXAoLjk5LDE2KSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBwYXJhbGxlbD1saXN0KGxvZ2luZm89VFJVRSwgZm9yd2FyZD1UUlVFKSkNCnN0b3BDbHVzdGVyKGNsKQ0KDQp0b2MoKQ0KDQpzYXZlUkRTKG9wdGltX3Jlc3VsdCwgZmlsZT0ib3B0aW1fcmVzdWx0X3YwMy5SRFMiKQ0KIyBvcHRpbV9yZXN1bHQgPC0gcmVhZFJEUygib3B0aW1fcmVzdWx0X3YwMC5SRFMiKQ0KYGBgDQoNCg0KYGBge3J9DQojcHV0IHRoZXNlIG5ldyB2YWxzIGludG9zcGVjaWVzX3BhcmFtcyBhbmQgZ28gYmFjayB0byB0aGUgdG9wIG9mIHRoaXMgcGFnZSB0byByZS1jaGVjayB0aGUgY2FsaWJyYXRpb24gDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJFJfbWF4PC0xMF5vcHRpbV9yZXN1bHQkcGFyWzE6MTZdDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGVyZXBybzwtb3B0aW1fcmVzdWx0JHBhclsxNzozMl0NCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkaW50ZXJhY3Rpb25fcmVzb3VyY2UgPC1vcHRpbV9yZXN1bHQkcGFyWzMzOjQ4XQ0KDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJFJfbWF4DQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGVyZXBybw0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRpbnRlcmFjdGlvbl9yZXNvdXJjZQ0KYGBgDQoNCmBgYHtyfQ0Kc2ltX3YwOSA8LSBwcm9qZWN0KHBhcmFtc19vcHRpbSwgdF9tYXggPSAyMDApDQpwbG90QmlvbWFzcyhzaW1fdjA5KQ0KYGBgDQoNCg0KYGBge3J9DQpwYXJhbXNfdjA3IDwtIHBhcmFtc19vcHRpbQ0KaW5pdGlhbE4ocGFyYW1zX3YwNykgPC0gc2ltX3YwOUBuW2RpbShzaW1fdjA5QG4pWzFdLCxdDQpgYGANCg0KDQpgYGB7cn0NCnNpbV92MTAgPC0gcHJvamVjdChwYXJhbXNfdjA3LCB0X21heCA9IDIwMCkNCnBsb3RCaW9tYXNzKHNpbV92MTApDQpwbG90RGlldChwYXJhbXNfdjA3KQ0KYGBgDQoNCmBgYHtyfQ0KaW5pdGlhbE4ocGFyYW1zX3YwNykgPC0gc2ltX3YxMEBuW2RpbShzaW1fdjEwQG4pWzFdLCxdDQpgYGANCg0KDQoNCmBgYHtyfQ0Kc2ltX3YxMSA8LSBwcm9qZWN0KHBhcmFtc192MDcsIHRfbWF4ID0gMjAwKQ0KcGxvdEJpb21hc3Moc2ltX3YxMSkNCmBgYA0KDQoNCmBgYHtyfQ0KcGFyYW1zX3R1bmVkX3YwNSA8LSB0dW5lUGFyYW1zKHBhcmFtc192MDcpDQoNCnBhcmFtc190dW5lZF92MDUgPC0gc3RlYWR5KHBhcmFtc190dW5lZF92MDUpDQoNCnBhcmFtc190dW5lZF92MDVAc3BlY2llc19wYXJhbXMkZXJlcHJvDQpgYGANCg0KDQpgYGB7cn0NCiMgc2F2ZVJEUyhwYXJhbXNfdHVuZWRfdjA1LCAicGFyYW1zX29wdGltX3YwMi5yZHMiKQ0KDQpwYXJhbXNfdHVuZWRfdjA1IDwtIHJlYWRSRFMoInBhcmFtc19vcHRpbV92MDIucmRzIikNCmBgYA0KDQoNCg0KYGBge3J9DQpzaW1fdjEyIDwtIHByb2plY3QocGFyYW1zX3R1bmVkX3YwNSwgdF9tYXggPSA1MDApDQoNCmluaXRpYWxOKHBhcmFtc190dW5lZF92MDUpIDwtIHNpbV92MTJAbltkaW0oc2ltX3YxMkBuKVsxXSwsXQ0KDQpzaW1fdjEzIDwtIHByb2plY3QocGFyYW1zX3R1bmVkX3YwNSwgdF9tYXggPSA1MDApDQpgYGANCg0KYGBge3J9DQpwbG90bHlCaW9tYXNzKHNpbV92MTMpDQpwbG90QmlvbWFzc09ic2VydmVkVnNNb2RlbChzaW1fdjEzKQ0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWxDdXN0b20oc2ltX3YxMykNCmBgYA0KDQpHcmFkdWFsbHkgcmVkdWNpbmcgcmVzb3VyY2UgbWF4aW11bSBzaXplIGluIHR1bmVQYXJhbXMoKSwgYXMgZG9pbmcgaXQgZGlyZWN0bHkgaW4gdGhlIHNldFBhcmFtcygpIHJvdXRlIHdpbGwgdGVsbCB5b3UgdGhlIHdfcHBfY3V0b2ZmIGhhcyBjaGFuZ2VkLCBidXQgaXQgd2lsbCBzdGlsbCBpbmNvcnBvcmF0ZSBhIGJhY2tncm91bmQgcmVzb3VyY2UgdXAgdG8gdGhlIHdfcHBfY3V0b2ZmIHRoYXQgd2FzIG9yaWdpbmFsbHkgdXNlZCBpbiBuZXdNdWx0aXNwZWNpZXNQYXJhbXMoKQ0KDQpgYGB7cn0NCnBhcmFtc190dW5lZF92MDYgPC0gdHVuZVBhcmFtcyhwYXJhbXNfdHVuZWRfdjA1KQ0KDQpwYXJhbXNfdHVuZWRfdjA2IDwtIHN0ZWFkeShwYXJhbXNfdHVuZWRfdjA2KQ0KDQpwYXJhbXNfdHVuZWRfdjA2QHNwZWNpZXNfcGFyYW1zJGVyZXBybw0KcGFyYW1zX3R1bmVkX3YwNkBzcGVjaWVzX3BhcmFtcyRSX21heA0KDQpgYGANCmBgYHtyfQ0KcGFyYW1zX3R1bmVkX3YxMCA8LSB0dW5lUGFyYW1zKHBhcmFtc190dW5lZF92MDYpDQoNCnBhcmFtc190dW5lZF92MTAgPC0gc3RlYWR5KHBhcmFtc190dW5lZF92MTApDQoNCnBhcmFtc190dW5lZF92MTBAc3BlY2llc19wYXJhbXMkZXJlcHJvDQpwYXJhbXNfdHVuZWRfdjEwQHNwZWNpZXNfcGFyYW1zJFJfbWF4DQoNCnBhcmFtc190dW5lZF92MTBAcmVzb3VyY2VfcGFyYW1zJHdfcHBfY3V0b2ZmDQoNCmBgYA0KDQoNCmBgYHtyfQ0KIyBzYXZlUkRTKHBhcmFtc190dW5lZF92MDYsICJwYXJhbXNfb3B0aW1fdjAyX3dfcHBfMTAwLnJkcyIpDQojIHNhdmVSRFMocGFyYW1zX3R1bmVkX3YxMCwgInBhcmFtc19vcHRpbV92MDNfd19wcF8xMDAucmRzIikNCg0KYGBgDQoNCg0KYGBge3J9DQojIGNyZWF0ZSBzZXQgb2YgcGFyYW1zIGZvciB0aGUgb3B0aW1pc2F0aW9uIHByb2Nlc3MNCnRpYygpDQoNCnBhcmFtc19vcHRpbSA8LSBwYXJhbXNfdHVuZWRfdjEwDQoNCnZhcnkgPC0gYyhsb2cxMChwYXJhbXNfb3B0aW1Ac3BlY2llc19wYXJhbXMkUl9tYXgpLA0KICAgICAgICAgIHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRlcmVwcm8sDQogICAgICAgICAgcGFyYW1zX29wdGltQHNwZWNpZXNfcGFyYW1zJGludGVyYWN0aW9uX3Jlc291cmNlKQ0KDQpwYXJhbXNfb3B0aW08LXNldFBhcmFtcyhwYXJhbXNfb3B0aW0pDQoNCiMgc2V0IHVwIHdvcmtlcnMNCm5vQ29yZXMgPC0gcGFyYWxsZWw6OmRldGVjdENvcmVzKCkgLSAxICMga2VlcCBzb21lIHNwYXJlIGNvcmUNCmNsIDwtIHBhcmFsbGVsOjptYWtlQ2x1c3Rlcihub0NvcmVzLCBzZXR1cF90aW1lb3V0ID0gMC41KQ0Kc2V0RGVmYXVsdENsdXN0ZXIoY2wgPSBjbCkNCmNsdXN0ZXJFeHBvcnQoY2wsIHZhcmxpc3QgPSAiY2wiLGVudmlyPWVudmlyb25tZW50KCkpDQpjbHVzdGVyRXZhbFEoY2wsIHsNCiAgbGlicmFyeShtaXplckV4cGVyaW1lbnRhbCkNCiAgbGlicmFyeShvcHRpbVBhcmFsbGVsKQ0KfSkNCg0Kb3B0aW1fcmVzdWx0IDwtIG9wdGltUGFyYWxsZWw6Om9wdGltUGFyYWxsZWwocGFyPXZhcnksZ2V0RXJyb3JDdXN0b21fdjMscGFyYW1zPXBhcmFtc19vcHRpbSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdCA9IHBhcmFtc19vcHRpbUBzcGVjaWVzX3BhcmFtcyRiaW9tYXNzX29ic2VydmVkLCBtZXRob2QgPSJMLUJGR1MtQiIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBsb3dlcj1jKHJlcCgtMTUsMTYpLHJlcCgxZS03LDE2KSxyZXAoLjEsMTYpKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVwcGVyPSBjKHJlcCgxNSwxNikscmVwKDAuOTksMTYpLHJlcCguOTksMTYpKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHBhcmFsbGVsPWxpc3QobG9naW5mbz1UUlVFLCBmb3J3YXJkPVRSVUUpKQ0Kc3RvcENsdXN0ZXIoY2wpDQoNCnRvYygpDQoNCnNhdmVSRFMob3B0aW1fcmVzdWx0LCBmaWxlPSJvcHRpbV9yZXN1bHRfdjA0LlJEUyIpDQoNCiMgb3B0aW1fcmVzdWx0IDwtIHJlYWRSRFMoIm9wdGltX3Jlc3VsdF92MDQuUkRTIikNCmBgYA0KDQoNCmBgYHtyfQ0KI3B1dCB0aGVzZSBuZXcgdmFscyBpbnRvc3BlY2llc19wYXJhbXMgYW5kIGdvIGJhY2sgdG8gdGhlIHRvcCBvZiB0aGlzIHBhZ2UgdG8gcmUtY2hlY2sgdGhlIGNhbGlicmF0aW9uIA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heDwtMTBeb3B0aW1fcmVzdWx0JHBhclsxOjE2XQ0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm88LW9wdGltX3Jlc3VsdCRwYXJbMTc6MzJdDQpzcGVjaWVzX3BhcmFtcyhwYXJhbXNfb3B0aW0pJGludGVyYWN0aW9uX3Jlc291cmNlIDwtb3B0aW1fcmVzdWx0JHBhclszMzo0OF0NCg0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRSX21heA0Kc3BlY2llc19wYXJhbXMocGFyYW1zX29wdGltKSRlcmVwcm8NCnNwZWNpZXNfcGFyYW1zKHBhcmFtc19vcHRpbSkkaW50ZXJhY3Rpb25fcmVzb3VyY2UNCmBgYA0KDQoNCg0KYGBge3J9DQpzaW1fdjE0IDwtIHByb2plY3QocGFyYW1zX29wdGltLCB0X21heCA9IDIwMDApDQpwbG90QmlvbWFzcyhzaW1fdjE0KQ0KYGBgDQoNCg0KYGBge3J9DQpwYXJhbXNfdjA4IDwtIHBhcmFtc19vcHRpbQ0KaW5pdGlhbE4ocGFyYW1zX3YwOCkgPC0gc2ltX3YxNEBuW2RpbShzaW1fdjE0QG4pWzFdLCxdDQpgYGANCg0KDQpgYGB7cn0NCnBhcmFtc190dW5lZF92MDcgPC0gdHVuZVBhcmFtcyhwYXJhbXNfdjA4KQ0KDQpwYXJhbXNfdHVuZWRfdjA3IDwtIHN0ZWFkeShwYXJhbXNfdHVuZWRfdjA3KQ0KDQpwYXJhbXNfdHVuZWRfdjA3QHNwZWNpZXNfcGFyYW1zJGVyZXBybw0KcGFyYW1zX3R1bmVkX3YwN0BzcGVjaWVzX3BhcmFtcyRSX21heA0KDQojIHNhdmVSRFMocGFyYW1zX3R1bmVkX3YwNywgZmlsZT0icGFyYW1zX29wdGltX3YwNF93X3BwXzEwMC5yZHMiKQ0KDQpgYGANCg0KDQpgYGB7cn0NCnNpbV92MTUgPC0gcHJvamVjdChwYXJhbXNfdHVuZWRfdjA3LCB0X21heCA9IDIwMDApDQpwbG90QmlvbWFzcyhzaW1fdjE1KQ0KcGxvdERpZXQocGFyYW1zX3R1bmVkX3YwNykNCmBgYA0KDQoNCmBgYHtyfQ0KaW5pdGlhbE4ocGFyYW1zX3R1bmVkX3YwNykgPC0gc2ltX3YxNUBuW2RpbShzaW1fdjE1QG4pWzFdLCxdDQoNCnNpbV92MTYgPC0gcHJvamVjdChwYXJhbXNfdHVuZWRfdjA3LCB0X21heCA9IDUwMCkNCnBsb3RseUJpb21hc3Moc2ltX3YxNikNCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsKHNpbV92MTYpDQpgYGANCg0KYGBge3J9DQpwYXJhbXNfbG9vcCA8LSBwYXJhbXNfdHVuZWRfdjA3IHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpIHw+DQogICAgbWF0Y2hCaW9tYXNzZXMoKSB8PiBzdGVhZHkoKSB8PiBtYXRjaEJpb21hc3NlcygpIHw+IHN0ZWFkeSgpDQoNCnBhcmFtc19sb29wQHNwZWNpZXNfcGFyYW1zJGVyZXBybw0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWwocGFyYW1zX2xvb3ApDQoNCmBgYA0KDQoNCmBgYHtyfQ0Kc2ltX3YxNyA8LSBwcm9qZWN0KHBhcmFtc19sb29wLCB0X21heCA9IDUwMCkNCnBsb3RseUJpb21hc3Moc2ltX3YxNykNCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsKHNpbV92MTcpDQpwbG90QmlvbWFzc09ic2VydmVkVnNNb2RlbEN1c3RvbShzaW1fdjE3KQ0KcGxvdEJpb21hc3NSZWxhdGl2ZShzaW1fdjE3KQ0KcGxvdERpZXQocGFyYW1zX2xvb3ApDQpgYGANCg0KYGBge3J9DQppbml0aWFsTihwYXJhbXNfbG9vcCkgPC0gc2ltX3YxN0BuW2RpbShzaW1fdjE3QG4pWzFdLCxdDQpgYGANCg0KDQpgYGB7cn0NCnNpbV92MTggPC0gcHJvamVjdChwYXJhbXNfbG9vcCwgdF9tYXggPSAxMDAwKQ0KcGxvdGx5QmlvbWFzcyhzaW1fdjE4KQ0KcGxvdEJpb21hc3NPYnNlcnZlZFZzTW9kZWwoc2ltX3YxOCkNCnBsb3RCaW9tYXNzT2JzZXJ2ZWRWc01vZGVsQ3VzdG9tKHNpbV92MTgpDQpwbG90QmlvbWFzc1JlbGF0aXZlKHNpbV92MTgpDQpwbG90RGlldChwYXJhbXNfbG9vcCkNCmBgYA0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg==